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ABSTRACT 

This report addresses research into active (feedback) control of 
the Satellite Power Astern (SPS), a large oechanically flexible aerospace 
structure. The basic objective of this analysis is to develop an overall 
assessment of the interaction of the SPS structure and its active 
control systems. 
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Executive Si:.wiiB»T»y! NAS-9-l6053 


In order to take maxlsnm advantage of the apace environoent, especially 
the available solar energy, future spacecraft and satellites viU tend to 
be large (0.1 - 100 kilometers in their largest dimensions) and very mech- 
anically flexible structures. It is likely that such large space structures 
(LSS) will exhibit the foUoving characteristics: 

(a) high total mess but lov mass density (mass distributed throughout 
the structure over a large area) ; 

(b) a very large number of low frequency resonances (many below 
1 hz); 

(c) very low natural damping (hence, easily excited resonances) 

Because of their high level of flexibility, they will tend to react to 
forces, such as thermal and gravity gradients, solar press'ore, and small 
impulse loadings, in ways which were heretofore ignorable in smaller struc- 
tures. In partici;lar, a concept like the Satellite Power System (SPS) 
will be especially vulnerable to these distrubance forces end yet, simiU.- 
taneousiy, be expected tc meet stringent performance requirements for 
shape fidelity, orientation alignment, vibration suppression, and pointing 
accuracy; however, all LSS will face similar difficulties. 

Active control must be used to augment such structures in order to 
meet these specifications. The basic element of active control is feedback ; 
sensors at various critical locations on the structure produce electronic 
signals proportional to local disturbance-induced errors in position, 
velocity, or acceleration; these signals are processed by control algor- 
ithms implemented in an on-line digital computer; the computer generates 
new electronic signals which are control commands applied to control 
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ACtuAtors such as thrust engines or toraulng devices: the control 
i Actustors, plAced At vArious locAtions on the structure, exert forces 

to bring the structure hAck into line And reduce the effects of Any dis- 
turhAnces. Althou^ the plAcement of AppropriAte ActuAtors end sensors 
^ is An iiBportAnt itea, the control Algorithm , isqpleaented in reAl>tia» by 

the on-line computer, is the principle component in this feedbeck link 
from sensors to actuators; it must be carefully designed to carry out its 
‘ complex tasks vithin the physical limitations of the actuators and sensors 

and the canputational limitations of the computer. 

The essential element of the control algorithm is a model of the j 

j 

structural dynamics . This model is used to predict the behavior of the J 

structure and compare it vith the actual sensor measurements. It uses 

any errors it detects to improve its ;>ercepticn of the instantaneous 

structiiral behavior and generate control commands to cause this behavior 

to meet the system specifications, e.g. to keep the SPS microwave system 

accurately pointed at the Earth. For the rather rigid structures which 

comprise current spacecraft ana satellites, such models are relatively 

simple; however, LSS require very complicated sets of partial differential 

equations and corresponding boundary conditions to predict the behavior \ 

at every location on the structure. A LSS is an example of -v^at is called 1 

a distributed narameter system (DPS) which in theory must be described by 

an infinite number of modes of vibration, but even in practice a very 

large number of these modes is required to adequately approximate the 

djTiamical behavior of the structure. 

The number of modes included in the structure model is related to the l 

computer memory capacity and its access time, i.e. how fast the computer j 




CUB cftleulatt. If only a few modes are included in the controller's 
a»del, the ccoputer can do the calculations very fast and there is 
ia^rceptihle delay in the feedback system. However, if a large number 
of modes is included in order to accurately model the structural behavior, 
the computer must respond more slowly and the actuator feedback commends 
are delayed and arrive too late to correctly reduce the distrubance effects. 
In some cases, they arrive so late that they increase the disturbance 
effects and cause the system to be unstable. This dilemma of how many 
nodes to keep (and which ones to throw out) and still be able to compute 
control c omm a n ds rapidly is called model reduction ; it is the crux of the 
LSS active control problem. Even if all the low frequency nodes of a 
structure like SPS were retained in the controller model, they would 
more than exhaust the a'\niilable computer capacity. Consequently, a model 
reduction must be performed and a trade-off made: Some modes are retained 

and others (thought to be less critical for performance) are left out of 
the controller's model. 

liow we cone to the controller-structure interaction problem. By 
necessity, we are forced to design a controller which actively controls 
the structure using a reduced model. The controller computer is only 
avsire of the modes included in its model of the structure. The residual 
(i.e. unmodeled) modes are an intrinsic part of the actual structure, but 
the controller has no internal knowledge of them. When a particular 
sensor detects a local displacement error, this error is made up of contri- 
butions from all modes - modeled and residual; the part of the sensor out- 
put caused by residual modes is called observation snillover . The control 
computer accepts the sensor signal and compares it with its reduced internal 
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nod«l to produce control coo&azids. Such a command causes a force to he 
applied at an actuator; this force excites all modes • modeled or residual. 
The excitation of the resld'ial nodes is called control spillover ( control 
energy "spills over" into the residual modes). 

Consequently, through the structure controller ve have feedback 
around the irodeled s»des of the structure (vhich ve want) and feedback 
around the unmodeled modes through the control and ohservation spillover 
terms. This latter Is the controller-structure Interaction and it must 
not he ignored in the design and evaluation of active large space structure 
control. It may sound foolishly ohvious hut an actively controlled struc- 
ture is quite different from a passive structure; hy electronic feedback, 
the active controller changes the structure's characteristics, and its 
response to disturbances. It is intended that this change vill improve the 
performance of the structure; however, if the controller-structure inter- 
action assessment is not considered as a fundamental part of any active 
control design for LSS such as S?S, the sought after performance improve- 
ments vill be lost. In fact, in some cases, the controller interaction 
with the unmodeled modes can cause them to be unstable which would mean an 
eventual loss of control of the struct'ore. 

This report presents a "top-down" approach to the design and analysis 
of active control systems for large aerospace structures, like the SPS. 

The introduction in Section I describes the LSS control issues smd the 
objectives of this study in greater depth and gives some current references 
on the general problem. Sections II and III describe the main research 
rerfcrmed: 
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(a) development of a mathematical framework for the LSS control 
problem; 

(b) development of iMthods to synthesize active structure con- 
trollers from distributed parameter or finite element models 
of the structure; 

(c) development of methods to analyze and assess the controller- 
structure interaction and its impact on the successful 
operation of actively controlled structures. 

Section IV describes the application of this approach to a generic model 
of the SPS which emphasizes the most stringent SPS control requirement: 
accurately pointing the microwave power transmission sub-system. Section V 
proposes and assesses methods to compensate and reduce the deleterious 
effects of the controller-structure interaction; this discussion includes 
structure, as well as controller, modifications. Section VI discusses 
an important issue for actively controlled structures: the control is only 

as good as the model. We have already seen the effects of the unmodeled 
modes; however, even the modeled modes may be ratner poorly described due 
to our lack of knowledge of the zero gravity behavior cf large structures. 
Consequently, some of the parameters of the structure may be poorly modeled; 
hence, there is a need for cn-board estimation cf structural parameters and 
the ability of the control algorithm to make use of these parameter esti- 
mates to adapt (i.e. t\ane) itself to the structure it is controlling. 
Sections VII and VIII present discussions and conclusions of the results 
of this study and recommendations for future research directions both in 


systems theory* and experiment design. The papers published 'under this 


contract appear as Appendices A-S and G; Appendix F contains a paper which 
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I. imODUCTIOR TO THE BPS kCSTTS COSTROL 


1.0 Introduction 

Th« Sntcllito Povtr System (SPS), as described in [l], is an extremely 
large aerospace structure « its aeration will require integration of struc- 
tural design and active control systems to generate and transmit, via 
microfv^ves, large quantities of eleetzlcal pover from geosynchronous orbit 
to receiving antennae on the earth's surface. The solar collector portim 
cf the satellite is the largest piece of the structure, see [l] Figure 
it has a lav structural mass (in proportira tc the total satellite mass) and 
is mechanically flexible vith many low frequency resoiances. In addition, 
the sicrovave transmitting antenna, although much smaller in size, is a 
flexible substructure uhich must be fine-pointed at the earth with minimal 


structural dynamic interaction with the solar arra;'. 

The following areas are most likely to i*equire sob» active control tc 

a\:gment the structural design during in-space operation; 

(a) attitude control for solar orientatior cf the array structure; 

!b) shape or figure control cf the array structure to maintain the 
nominal form; 

(c) fine earth pointing control cf the microwave transmission 
antenna- 


Active (or feedback) control would be iarlenented by on-line digital pro- 
cessing ^f signals, from a variety of sensors located or. the structure, tc 
produce control actuator commands, for thrust and momentum devices, to make 
the necessary corrections to achieve the above control objectives. These 
control tasks might be achieved with separate, or ''decentralized", control- 
lers or by a central control processor or some combination cf these; the 
tasks would have tc be accomplished simultanec'isly or. the SPS. 


2.0 f«ehaie>l Dltguiiloa 
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2.1 Protelta Stattaint 

Sine* tht S?S it t btebtnietUy fltxiblt itructurt, Itt ^Tutaic bt- 
htvlour it aodtltd by ^ eontlnvtua dttcriptioB, i.t. t ttt of coupled par- 
tial difftrential tquttioaa which aodcl the flesdblt tv^ttructurtt . *.g. 
btaat t plattt » ati^s^tt , and a ttt of boundary conditiont which aodtl 
the inttreozmtetion of thttt tubitructurti to cash ether and to rigid aaaaea 
in the total STO ttructure. Ihit scans that SPS it in the category of 
distribute d narsaeter tvatena (DPS), i.e. those described by partial dif- 
ferential equations as distinguished froa luased oartaeter systeas which 
are described by ordinazy differential equations. 

Such a DPS description of SPS would have the fora: 


where u(x,t) represents the (vector of) displaceaeats (rotational and 
trtaslational) froa the equilibriuB position of the structure H at the 
point X aiid tiae t due to the external force disbributicn ?(x,t). The struc- 
tural intemal restoring forces A u sad dsspiag D^u* ore generated by partial 

O w w 

differential operators I^with appropriate boundary' conditions on u, e.g. 


A jU ■+■ AgU* ■ 0, on bo'iaidary of H 

The aass distribution of the structure is given by a(x). The force 
bunion is coa^osed of disturbance forces and control forces 
F(x,t) ■ Fjj(x,t) ♦ Fj,(x,t) 


r * <5 'i 


distri- 


(2.3) 


Control aay be introduced via !•'. force or acaentua actuators with 
aaplitudes f, (t): 

A 


ft « 



f A 
V ^ < 


where b,(x) are the inf?.uence f-jsetiens of 


the indivi;. 


actuators and 
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y(t) ■ CgU^ 


f(t) It tht vtctor of actuAtor conaan dt Observation of the tyataa 

it dona via P tentort: 

(2.5) 

vhert dia y(t) ■ P and (2.5) allcvt fc position, velocity, acceleration, or 
CM^inationt of tbete aetsureBentt . There would alto be tensor and actuator 
tl^aaict associated with (2.1i}>(2.5) but these dyntaics have net been shown; 
they would be part of the assesssant. 

^e perforaance desired of the SPS in tcras of pointing accuracy, shape 
control, and solar orientation aay be stated in teras of clssed*locp tran- 
sient and stsa^-state behavior or as a perforaance function, such as modal 
energy' or pointing error, to be minimized. 3cth t:*pes of perforaance speci- 
fication will be taken into account in this assessaent. 

From DPS theory, the ^naaical behavior of the uncontrolled SPS stric- 
ture is described by an infinite nuaflser of vibration modes - the rigid 
bod” (zero frequency) aedes plus all the independent (resonant frequency) 
flex *odes. In practice, finite eleaent struct -ural anal;.'sis techniques 
are used to produce approxiaate aodal data (node frequencies and aode shapes) 
which yields a large-scale, luaped pareaeter rodel of the £?£ (i.e. a l^gc 
but net infinite r.'uaber cf modes'. Tnis large sctle model suffers from 
Bodel error (i.e. error in the aodal data) and uruce deled residual aodes cf 
the SPS. 

In addition, although the large scale acdel ctn be used for analysis 
and ccaputer simulation purposes , it is usually too large to use for an 
on-board, real-tia^ controller. 7ne limited capacity :f such controllers 
will req*uire that a reduced- o rder mode- cf the SPS be produced; this en- 
tails selection of a saallcr subset of modes (or an aggregation of aodes) 
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t^iat gives a reasonably faithful representation of the total system be> 
havior. The active controller design would be based on this reduced-order 
model. Such a controller would take the form: 
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{ 


f(t) ■ y(t) + z(t) 


z(t) ■ y(t) + Hgg z(t) 


( 2 . 6 ) 


where din z(t) ■ S < ». The controller dimension S and the parameters H.,, , 

•kX 

^12’ ^21’ ^22 determined directly from the reduced-order model 

of the SPS. The size of S and the operations involved in (2.6) are directly 
related to the capacity of the on-board control computer (or subcomputer) . 


2^2 Fundamental Research Problesis for Active Control of SPS 

The msst fundamental control research problem for the SPS is the 
interaction of any on-board active controller (2.6), based on a reduced- 
order model, with the actual SPS structure (2.1)-(2.5). This structure- 
controller interaction occurs becaiise of modeling error, control-device 
dj'naaics, and spillove.- due to modes not included in the reduced-order 
model (for "spillover", see [2]-[3]); this interaction is the most crucial 
research issue to be considered in this analysis. Other important issues to 
be assessed are the following: 

(1) selection of appropriate structure vibration nodes to control to 
achieve desired system performance; 

(2) the effect of residual (uncontrolled) nodes on system performance; 

(3) the effect of modeling error (i.e. modal approximation) on system 
performance; 

(1) the use of appropriate compensation, if needed, to counteract 
residual mode and model error effects 
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(3) the affect of ectiMitor end sensor dyneoics on controller design 
end systea perforaeneei 

(6) controller design with multiple distributed ectuetors end sensors; 

(T) control device (ectuetor or sensor) location, quantity, and char- 
acteristics necessary to achieve the basic control objectives ; 

(8) examination of the trade-offs available in the basic control 
approaches, including stractural design, passive dampers, direct output 
feedback, and modal control with condensation feedback and 

(9) separation control of the SPS array structure from the microwave 
antenna to reduce the effect of vibrations in the main structure on the 
beam-pointing accuracy. 

2.3 Aonroach 

This analysis will use the flexible structure control theory of 
[2]-[6]to develop a genereil mathematicad framework for the assessment of 
basic active control systems procedures for the SPS. The development of 
this framework will be the principle element of this assessment. It will 
be ible to begin with a continuum description of SPS, or a finite element 
description, or even a simplified continuum model, such as a plate model, 
and fellow through, in a "tep-dov-" manner, to the basic controller logic 
and provide an assessment of the fundamental issues raised in Section 2.2. 
Particular erphasis will be given to making the framework broad enough to 
include new control trends, e.g. [T]-[9], and, yet, capable of yielding 
answers to the above issues and to the assessment of ether approaches to 
(reduced-order) modeling, e.g. traveling-wave vs. modal descriptions of 
SPS. It will not be limited to the current SPS reference system [l]. 
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This approach will emphasise: 

(a) generic problems of SPS control; 

(b) the basic design logic for controller sj*nthesis; 

(c) the trade-offs inherent in the basic control approaches for the 
attitude, figure^and pointing control tasks* 

(d) the interrelationship and interaction of the overall structural 
aesign and its models with the active controller design end operation. 

(e) the distributed parameter nature of the modeling, control and 
interaction problems. 
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Our tasic objectives in this study are to develop fundsaental con- 
troller synthesis approaches for SPS structure control and to assess the 
interaction of the SPS structure vith its active controller(s). An innle- 
nentahle active structure controller is one which can make use of a vide 
variety of possible sensors and actuators to synthesize on-line . via 
digital computer (s), the feedback control commands necessary to maintain 
the structural shape, attitude, and/or pointing req.uirements in the pres- 
ence of disturbances. Such a controller must be finite-dimensional. The 
controller dimension must not be too high because this dimension corre- 
sponds to the on-line computer capacity (i.e., the size of the memory and 
the speed of memory access). 

The fundamental issue is the finite-dimensional . ianlementable con- 
trol of an infinite-dimensional , distributed parameter structure. Any 
controller synthesis technirue must produce finits-dime^sionaJ controllers 
for the distributed parameter structure; this will always involve some 
kind of model reduction of the infinite dimensional system. However, due 
tc this model reduction, there will always be a ccntrcller-structure in- 
teraction with residual (unmodeled) structure dynamics. In general, the 
sources of controller-structure interaction are the following: 

(a) structure modeling error (i.e., poorly known model parameters) 

(b) unmodeled control device dj'namics 

(c) spillover due to structure dynamics unknown to the controller. 

Therefore, it is not enough to produce active controller synthesis 

procedures for structural control; the controller-structure interaction 
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BUft be analyzed for any auch procedure. This analvala la an Intrlnilc 
part of the uroblea and abould be done with a view toward controller re- 
design for reduced interactions. Without a general aathematical fraae- 
work which takes the above issues into account, it is iapossible to 
handle the active structure control problem in any but a trial and error 
fashion. 

In Paper A: Feedback Control of Dissipative Hyperbolic Distributed 

Parameter Systems with Finite Dimensional Controllers (presented at the 
19T9 USE Control and Decision Conference; revised and expanded version 
submitted to SIAM Journal of Control and Optimization) we have presented 
such a general mathematical framework for model reduction, controller 
synthesis, and interaction analysis. This framework is, used to assess 
modal control approaches to active structures. In modal control , it is 
assumed that the vibration mode frequencies and mode shapes of the struc- 
ture are known. A group of critical modes are chosen, e.g. , by their 
excitability by the disturbance frequency spectrum cr by their contribu- 
tion to a performance index such as the microwave antenna pointing 
accuracy; the choice of critical modes is not a trivial issue by any 
means and merits further investigation. These critical modes form the 
reduced-order model of the structure. The nodal controller is synthe- 
sized with only the critical mode data. In Paper A, the most general 
linear controller is considered; this allows trade-offs in the controller 
design, the number and location of ccr.trcl devices (actuators and 
sensors), and the on-line computer capacity. Very simple modal criteria 
are developed to assist the designer in assessing the available design 
trade-offs. The controller- structure interaction is analyzed, and 
bounds on spillover and mod*; error are produced; these bounds indicate 






Q 


the tolerable levels of these interaction tern, for stable operation 


of the active structure. When the bounds are satisfied it makes sense 
to use the modal controller as designed; vhen they are not satisfied, 
the designed controller should be used with caution and, more reasonably, 
it should be redesigned or spillover compensation should be added. We 
viU say more later about spillover ccopensation. The mathematical 
fraaevork presented in Paper A is very general and can handle many other 
types of active structure control; it is not limited tc modal control . 

The most ccaoaon tool a'^^lable to structural dynamicists is NASTRN 
which is a large-scale computer code based on a finite element analysis 
of a structure. Approximate node shapes and frequencies are available 
from NASTPN auid can be used in active controller synthesis. The fact 
that the modal data is approximate adds the effect of model error to that 
of spillover to further complicate the controller- structure interaction. 
In Paper E: Finite Element Models and Feedback Control of Flexible Aero- 

spaoe Structures, the structure con’rol synthesis and analysis problem is 
considered in detail for finite element models of :he structure . Bound.- 
on tolerable model error and spillover are produced to indicate when an 
implemeniable active controller based on a consistent finite element 


arrroximation of a struct-ure like S?S makes sense and can be exoected 


to operate stably . 

Among the practical constraints on structure control, we are inclined 
to assume that point actuators and sensors will be used; these may be 
placed at many locations on the structure. In some special situations 
a distributed actuator, e.g., use of electromagnetic field interaction 
of the SPS and the earth, may be available. Such distributed actuators 
and sensors are not excluded from the framework established in Paper A 


and may be used to advantage when they exist. 
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It should h« mentioned that In all of our work on active structure 
control we are talking shout feedback or closed- loon control of the 
structure. In "Control-Structxare Interaction in a Free Beam” by S. 
Shrivastava and R. Ried (NASA Report JSC 16699, May 1980), an energy 
approach is used to study the controller-structure interaction for a 
free beam with on-off thruster control. The paper illustrates the 
excitability of structure modes to such nonlinear control forces. This 
is certainly an important consideration; however, one should be careful 
not to draw many detailed conclusions from such an analysis because the 
control is ooen-loox> . Many of the controller-struccure interaction 
effects appear because of the feedback nature cf active structure control. 
The analysis of Shrivastava and Ried aids in understanding the complex 
way in which a flexible structure responds to open-loop thruster loads; 
it is quite complementary to our analysis which considers linear forces 
generated by feedback control systems. 

The above mentioned papers A and E are A.ppendices A and B cf this 
report . 
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III. ACTIVE STRUCTURE CONTROL USING SINGULAR PERTURBATIONS 



In the previoue section » ve presented s genersl frsa»vork for setive 
structure control synthesis end snslysis. It ves based on a regular 
■perturbation vievpoint. i.e., the spillo'ver and model error terms vere 
taken as perturbations of the desired system; this is hov ve arri-ved 
at our spillover and model error bounds for closed-loo? stability. A 
different vievpoint vhich uses singular (as opposed to regular) perturba- 
tions is described and analysed in Paper C: Reduced-Order Feedback 

Ccntrol of Distributed Parameter Systems Via Singular Perturbations 
(presented at the 198O Princeton Conference on Information Sciences 
and Systems; revised version submitted to Journal of Mathematical Analy- 
sis and Applications) amd Paper D: Stability of Distributed Parameter 

Systems with Finite Dimensional Compensators Via Singular Pert'urbation 
Methods (presented at 1980 International Congress on Applied Systems 
Research and Cybernetics). Applications cf the general approaches to 
aerospace structures are discussed in Paper I: Clcsed-loop Stability 

cf Large Space Structures Via Singular and Regular Pert’orbation Tech- 
niques (presented at 1980 IEEE Control and Decision Conference). 

-"ic principal difficulty in applying the sing-ular perturbation 
approach is to arrive at a physically meaningful singular fona'ulaticn. 
When this can be done the system separates into a critical subsystem 
that must be controlled and a residual subsystem whose overall effect 
can be approximated statically. This makes the controller synthesis 
and air.alysis very different frcm that of Section II and we feel this 
Justifies a separate section devoted tc this approach- • ... 


‘•v. 


No single approach to the active structure control problem can be 
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expected to yield an the Magyrs . Each one has its a^cial advantages 
and disadvantages in a given a]^lieation. In this and the previous 
section, we have presented two very general viewpoints which incorpor- 
ate the controller synthesis and the interaction analysis which we feel 
is crucial to the solution of this complex problem. 

Papers C, D, and £ are Appendices C, D, and £ of this report. 


IV. CQSTRQLLER-STRUCTURE ICTmCTlOS FOB S?8: 
APPLICATIOlf TO A GEHERIC riODg. 


Although the ueeenient a^roaehtt deicrihed in the previous see* 
tions msy be applied to a vide variety of flexible struet.ure ludels, 
including finite eleaent aodels, ve have chosen to look carefully at 
their application to a sinple aodel of SPS and the (potentially) aost 
difficult control problea: accurately pointing the aicrovave antennae 

on the flexible solar array structure. Of course, mart detailed acdels 
will exhibit Bore aspects of this coaplex problea but we are operating 
on the trite-but-true principle: You have to start somewhere I 

The generic aodel we have used for SPS is a flexible bean with 
luaped aasses at both ends; this aodel was also considered in 
"Projected Models of Solar Power Systea" by J. Juang (Jet Propulsion 
Laboratory Engineering Meac 3^7-58, Feb. I960). 

The dynanics are given by: 


{ 


W w w w 

y - u ♦ Eq 


u, + A. u - B- f 


:ere 


(4.1) 


a u 


xxxx 


- S u. 


XX 


(h.2) 


with appropriate boundary conditions and u(x, t) is the transverse dis 
placement off equilibriua of the beam of length t . The dasping oper- 


ator Tp. is chosen later; its fcm is : 

V 

only adds a saall arcunt of dissipation 


t of great iarcrtar.ce since it 
0 the structure. We consider 


several actuator-sensor arrar.geaents. 


Tht PTfenswcc indtx !• given ty: 

J ■ I [E^(0, t) + E*(i?, t) ->• f* R f] dt (i*.3) 

■'0 

where the i>clntin£ error E(x, t) it given by; 

E(x, t) ■ e(x, t) - e^(x) 

where Sj.(x) it the desired pointine direction at the location x and 


9 (x. 


3v‘ x, t! 


3x 


is the angle of the ncraal 'or pointing) vector at location x at tiae 
t . The control penalty tera has a positive definite weighting aatrix 
B ; this penalty tera can be selected to heep the acxnator c onaan ds 
within practical liisits. 

Tne systea (-.1) together with the perforaance index (u.3) fonts 
a distributed ’^araaeter crtiaal control rro'cl^a ; this can be brought 
into the fora of an optiaal regulator by considering: 


e(x, t; 


u(x, t) - u (x) 
r 


( L 


where u^(x) is chosen so that' 


3u. 


3x 


It is easily shown that 


{ 


•tt * "0 * *‘0 ' 


• ^ 




y ■ =0 ' * ‘0 *t * =c > 
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vh«a (m it tht ettt by proptr choiet of ): 

*0“r’° 

vt ebtU ehoott to mctaurt 

y ■ y - Cq Uj, (1*.9) 

intttad of y ; this c&£ alviyi I > dont btctuit it t k&ovn tint- 

Invtritnt ttra. Alto, not« that tbt pointing trror btcoatt: 

S(x, t)-|| 

Combining (^.5), (**.7), (^*.9), (*..10), and (1*.3). w« obtain a 
diitributtd oaramater ontimal regulator problem , at dttircd. 7nit prob- 
lem may bt analyzed using tbe methods of Paper A. We have developed 
active n»dai controllers by tve straightforward methods: 

(a) produce the infinite dimensional controller (which cannot be 
implemented) and truncate it to obtain the finite dimensional ccntrolle: 

(b) truncate the distributed parameter system and performance 
index, and generate the finite dimensional controller from this reduced* 
order model. 

The latter method seems the nest natural from an engineering viei.*- 
point. When the exnct modes are known (at they are in this model), the 
two methods yield essentially the same controller. 

This application and the study of its properties form the basis fo; 
the masters thesis of my research assistant, K. l-!ashayekhi ; this thesis 
will be completed and available by 1, 1 


, 931 . 





SoM aUcviatlM of the cmtroUor-ttructuro inttroetion may ba 

* 

obtaiaad by puraly itruetural modification or radatign. As optimal ovar 
all structural dasign may ba dona to radistributa tb'^ mast and iscraasa 


tba stiflbats. From our control point of viav, wa naad incraasad daan?» 
ing vithout tarious loss of stiffnass ; ainca tha stiffnasi is m a in l y a 
function of configuration and tha damping is a function of materials and 
construction, this miy ba attainabla vith structural reconfiguration and 
use of ccaposita materials. Also, wa naad to know tha system parameters 
as well as possible; perhaps, an attempt should be made tc structiirally 


isolate subsystems or tc yield a structxira composed of wall-approximated 
subsystems. One thing is vary clear, the best contromise of structural 
stiffness vs. dai^ing would be made in the treliainar^* desien phase ; such 
decisions require a basic •understanding cf the effects cf strucforal con- 
figuration and materials tc meet the mission requirements. 

Structural add-ons include coatings to increase high frequency damp- 
ing and (passiTe or active; member da.tters to increase lov frequency 
damping. The latter are especially desirable because it is the low fre- 
quency damping that most needs augmentation; perhaps, some combination of 
dampers plus a more sophisticated active control can be used to achieve 
the structural control requirements. 

Of course, much of the structural design is constrained by factors 




ultimate functic 


;f the overall struofure, e.g.. trans- 


portability and ccnstracticn in orbit, nonetheless, in the more distant 
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future ve would expect to see control-configured structures » i.e., 
structure and active control designed and optimized as one integrated 
package. Such an Integrated optimization is not simply a matter of 
producing a "stper coo^uter program"; a much deeper understanding of 
structural behavior in the presence of active control will be rei],uired 
to generate control-configured aerospace Ftructures. 


e 


V 


In general) ve feel that structural oodification alone is quite 
limited as a means to reduce the controller-structure interaction. The 
interaction prohlos occurs because of the combination of mechanical 
structure and control electronics i this suggests further use of elec- 
tronic cos^ensation as an alternative. 

Actually, there are many vays to introduce this stillover compensa- 
tion . The basic requirements are the folloving: 

(a) c<m^ensation should be added-on to the existing controller 

(b) compensation shovild not substantially increase the ntoiber of 
control devices or the dimension of the controller (i.e., the on-line 
computer capacity). 

The methods we have investigated are prefiltering, cancellation, 
residual aggregation, and enhancement. Each of these has its advantages 
and disadvantages ; none is a panacea. The most natural seems to be 
■crefjlterin/?. i.e., filtering the sensor outputs tc. reduce the spillover 
frequency components; however, phase distortion due to the filter can 
have as deleterious an effect on the stability as the spillover it 
removes . 

We feel this area has not received sufficient attention; there are 
many promising directions for electronic compensation methods which 
need further investigation. 

3 . C Stecial Purocse flocalized) Controllers 

Decentralized control in general sounds like a good bet for active 
structure control. 'It has been used successfully in the area, of electric 




power generation and distribution; this application area has so&e sur- 
prising theoretical similarities to flexible aerospace structures. 
Therefore, ve do not want to discount the possible benefits of a de- 
centralized approach to structure control; however, ve also do not want 
to forget the extremely difficult technical problems that* continue to 
impede the development of a general theory of decentralized control. 

Still, a localized approach has a great deal of serit if ve exploit 
properties of the particular application. Specifically, the SPS micro- 
wave pointing problem mi^t be attacked in the following localized way: 
treat the flexible antenna as the precision system, i.e., the one in need 
of accurate control, and the solar array as a disturbance system acting 
upon the antenna; try to design a disturbance counteracting control to 
point the antenna and simultaneously reduce disturbances from the solar 
array structure. Such an approach was first considered in "An Active 
Modal Control System Philosophy for a Class of Large Space Structures,” 
by K. Balas and J. Canavin (presented at 19TT Blacksburg Symposium). 

Of course, a centralized controller with a complete dj'r.amic model 
of the array structure with antennae attached car. do a much better Job 
than the above-described localized control. However, the centralized 
controller is much more ccmplicaned and nay prove impractical to imple- 
ment while the more simple localized controller could be implemented 
and may prove to be nearly as effective as the centralized scheme. The 
trade-offs in such localized structure control schemes remain to be more 
carefully examined. 
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Active control of structures is based on a complete knowledge of 
certain system parameters i for exaiqpley in modal control, the mode fre- 
quencies and mode shapes must be known (at least, the mode shapes must 
be known at the actuator-sensor locations). In practice, these param- 
eters are poorly known and may vary with time. The control designer has 
two choices: 

(a) make the design robust, i.e., capable of carrying out the 
control task even though the parameters have substantial error; 

(b) make the design adaptive , i.e., use the sensor information to 
estimate poorly known parameters on-line and have the controller tune 
itself up. 

It is difficult to produce robust controllers for highly accurate 
tasks without introducing unacceptably high gains; therefore, we have 
concentrated on the adaptive approach. However, one should certainly 
check the sensitivity of a given structure controller to parameter inac- 
c\u:acies before proceeding with parameter estimation and/or adaptive 
control. 

A general approach to parameter estimation and adaptive control for 
aerospace structiares is presented in Paper F: Towacrd Adaptive Control 

of Large Structures in Space by M. Balas and C. R. Johnson, Jr. (from 
Applications of Adaptive Control , K. S. Rarendra and R. Monopoli, eds.. 
Academic Press, 19?0). This paper was not written under this contract; 
however, we include it because it discusses the special problems of this 
field within the context of our general mathematical framework. In 
addition, seme preliminary investigation (accomplished with RASA support 
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this year) on adaptive control cf reduced-order systems is presented in 
Paper G: Reduced-Order Adaptive Controller Studies by C. R. Johnson, 

Jr., and M. Balas (presented at I960 Joint Automatic Control Confer- 
ence ) . 

Papers F and G are Appendices F and G of this report. 
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VII. DISa^SSION AKS CONCLUSIONS 

No one who has vorked in the field of aerospace struct\u:es vould 
argue that active control of a structure as large and unruly as the 
proposed Satellite Power System is anything but a monumentally diffi- 
c\ilt tasx. Furthermore, no one but a fool or a madman vould undertake 
to solve the basic SPS control problems in a one-year study. Hoping 
not to be judged the former and unsure about the latter, ve have never- 
theless undertaken the study of some of these problems under this con- 
tract. 

Throughout this work, our emphasis has been on the generic problems 
and the distributed parameter nat/ore of SPS active structure control 
with the goal of obtaining a deeper understanding cf structure control 
technology by focusing on an extremely complicated and difficult structure. 
We have gone in many directions during this study; cfoen, in order to 
cover as much ground as possible, we have developed a given direction 
only to the point where it began to reveal ixs potential. Still, this 
work is not without substance. We have developed a mathematical frame- 
work for the assessment of the active control of aerospace structTires 
like ( although net restricted to ) the SPS. Prom this we have generated 
quite general techniques to: 

(a) synthesize structure controllers from distributed parameter 
or finite element models of the structure; 

(b) analyze the controller- struct'ore interaction on the closed-loop 
stability; 

(c) compensate . by struct’jre or controller mcdificaticn, to reduce 


the interaction effects. 


This has been done vithln the constraints of aaxlaug use of exist- 
ing control and structure knovledge and tools, allowance for a variety 
of SPS models and configurations, and ability to include new develop- 
B»nts in control and/or structure technology. 

We have accomplished our principal objectives in this study: 

(1) to develop basic active structure control approaches for SPS 
and assess the design trade-offs 

(2) to develop mathematical methods to assess the effects of con- 
troller-structure interaction on the closed-loop stability 

(3) to develop numerical methods to compute the stability bounds 
of (2). 

This work has produced seven published papers and a masters thesis. 
We have presented various aspects of the work in several national con- 
ferences and in seminars at L. B. Johnson Space Center and the Jet 
Propulsion Laboratory throughout the year. 

We hope our contributions and the perspective presented here will 
further the understanding of Satellite Power System control 
control of aerospace structures in general. 


and active 
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VIII. RDCOMMENPATIOKS FOR FUTURE RESEARCH IK LARGS AEROSPACE STBUCTURES 


The research reported here has been on the development of systems 
theory. It is In this area that ve shall make our principal recommenda- 
tions. 

Two areas, out of many which we have touched on in this report, 
stand out as major issues in the successful development of active struc- 
ture control: 

(1) Structural modeling for active control , i.e., modeling of the 
structure explicitly for feedback control purposes as opposed to assess- 
ment of dynamical behavior under external, loads; 

(2) On-line parameter estimation and adaptive control . 

The first of these really is an indication of the different view- 
points of the structural dynamicist and the structural control designer; 
emphasis on this problem would go a long way toward a much-needed 
rapprochement of these two disciplines. One cannot realistically expect 
tc operate large aerospace structures without some resolution of the 
second issue; yet, on-line parameter estimation and adaptive control 
techniq,ues lead to some of the most complex theoretical issues in dis- 
tributed parameter systems. 

There is a clear need to look beyond the systems theory aspects of 
large aerospace structures and to begin to consider the development of 
gr' und-based, and eventually in-space, experiments in actively controlled 
structures. Of course, no successful ground-based experiment can begin 
to approach the dimensions of something like the S?S; however, this is 
not necessary to gain the basic understanding that is needed tc cope with 
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this nev area of technology. 

Instead, we say begin to investigate particular structural behavior. 
For exBsple, ground-based experiments to determine the dynamical behavior 
of flexible membranes attached to s«&i-rigid structures would help us to 
understand a major coiQ>onent of the SPS: the interaction of the solar- 

cell blankets with their supporting structure. Furthermore, active 
suspension control could be introduced into such an experiment to inves- 
tigate its capacity to suppress the blanket vibrations without seriously 
increasing vibrations in the supporting structure - a clear application 
of the systems concepts developed in this report to design and assess 
the controller-structure interactions. 

Other experiments in active control of beams and plates could be 
designed and carried out to assess the potential of active structures 
and to illuminate the problems. Although computer simulation is an 
excellent tool for analysis of struct;iral behavior, it is only as good 
as the models of the phenomena to be studied. 7nese models require sub- 
stantial improvement in order to successfully predict the complex 
behavior of large aerospace structures. Such improvements can only 
grow out of a better understanding of structures and their interaction 
with control - an understanding which car. be greatly aided by laboratory, 
and space-shuttle-based, experiments with generic structures. 
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A feedback control of dissipative hyperbolic distributed parameter systems 

WITH finite dimensional CONTROLLERS 
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ABSTRACT 


The class of dissipative hyperbolic distributee parameter systens 
Is defined here; It includes many physical phenomena such as mechanically 
flexible structures and certain linear wave propagation problems. Feedback 
control theory for this class of systems Is investigated. Although some re- 
sults on infinite dimensional controllers are presented, the enphasis is on 
practical (I.e., finite d1n«ns1onal) controllers based on reduced-order 
models of the distributed parameter system and the closed-loop operation of 
such controllers in the actual system. Sufficient conditions, in terns of 
"spillover'' bounds for the unirodeled residuals, are oresenteo “or exponen- 
tially stable closed-loop operation. 

1.0 INTRCDUCTION 

Many physical phenomena are best described by partial dif-ierential 
epuations, e.g., [1], Chapt. 1, and hence, are distributed parameter systems 
(OPS). Several books have been -written on DPS and control, e.g., IV-'J], as 
■well as, a large number of papers, e.g., the sur/eys [5l-[5]; aoolications p* 
DPS control range from ir^c.hanically 'lexiple strjcfures to wilcfire s«pcres- 
sion, e.g., the survey [10]. Our special area of interest has been apolica- 
tion to large flexible aerospace structures, e.g., tne sur/eys Iil]-[12]; 


/ 
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although iht rtsults obtatntd In this paptr may bt appllad to such struc- 
turts, thty can ba applltd to a much largtr class of probltms, namtly, 
dissipative hyperbolic DPS. 


2.0 DISTRIBUTED PARAMETER SYSTEM DESCRIPTION; DISSIPATIVE HYPERBOLIC SYST!?<S 


We consider OPS describedby the following: 


■ A v(t) * B f(t), v(0) • Vq 

' ' (2.1) 

y(t) • C v(t) 

where the systen state v(t) is in a Hilbert space H with inner product (•,•) 
and associated norm H*M* "Hie differential operator A is linear, time- 
invariant, has domain 0(A) dense in H, and generates a Cq semigroup U(t) 
on H. The OPS is called dissipative hyperbolic here when the semigroup 
U(t) has the following "growth'' property: 




i -0 • 




t > 3 


( 2 , 2 ) 


.^rere c > 0 and 1 1 ; when e > 0 , -we say the 3PS 1 s strictly dissisatl'/e . 
I* • C, '^e sss'wme ■ 1, For mere details on semigroucs, see [5j-[7j, 
[13]-[Uj. 

The control is introduced via M inputs: 


•4 

3 f(t) • ^ b, f^(t) 
i-1 
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wnere the sct'jator influence functions are in H. ’he system is ooserveo 
via P sensors whose outputs yj(t- form tne vector y{t; • Ly.j!t) yp(t)]‘ 
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and art givtn by 

yj(t) ■ (Cj, v(t)) (2.4) 

wiMrt thi stnsor Infiuenct functions Cj are In H. The rank of the Input oper> 
ator E Is M and that of the output operator C Is P. In most UPS applications 
the control and observation must be done with a finite number of devices; 
hence, these restrictions on S and C are natural. Usually, the control de- 
vices will be localized, I.e., the influence functions b^ and will be 
nearly Oirac delta functions. 

The class of dissipative hyperbolic systems Includes mechanically 
flexible structures, where t 1 $ related to the natural or materials damping 



In the structure, and certain wave propagation problems which are described 
by symmetric hyperbolic systems of partial differentlai equations, e.g., [15]; 
In fact. It Includes most oscillatory DPS. It does not Include (so called) 


parabolic problems like the heat equation; for these systems, the semigroup 
U(t) Is holomorphic [16] and the control oroblems are somewnat different. 

Control of parabolic DPS nas been considered In, for example, [5^, 

Our control aooroacn to dissipative hyperbolic DPS emphas-'ies 
plementabll 1 tv . I.e., from tne sensor cutouts y(t) i" [Z.~] , the controller 
must synthesize control commands f(t) In (2.3) to adequately staoilice ana 
Improve the performance of the OPS, and these commancs must be synthesized 
with an on-line conmuter. Consequently, the control synthesizer must be a 


finite dimensional dynamical system (usually cf low dimension). *his o’aces 
reasonable, practical constraints on tne DPS control crobltm [2C3. 



duced prcer modeling and implementable controllers are c'scussed in Secs. ~ 
and 3. Dur main results on closed-loco stacility '*-e presented in Sec. 3 wi t.n 


c 


application of these concepts to control of large flexible structures in 
Sec. 7 and general conclusions in Sec. 8. 
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3.0 INFINITE DIMENSIONAL CONTROLLER 


Although the results presented in this section may not be implem- 
table (and, in many applications, they are not), they do give some insight 
into the OPS control problem. 

Control laws for OPS, as in the case for limf)ed parameter systems, 
often take the form of gains multiplying full state feedback: 

f*(t) ■ G v(t) (3.1) 


where Q is a bounded control gain operator; this is especially true Wien a 
performance index is optimized, e.g., [6], Chapt. 14. However the states 
v(t) for a DPS are in an infinite dimensional space and, in most cases, can 
not be recovered with a finite number of sensors as in (2.4). Consequently, 
one is led to a state estimator (again, as in the lunged parameter case) of 
the form: 


» A y(t) + B f(t) + < y(t) 
y(t) » C v(t) 


y(t)), 0(G) « 0 


(3.2) 


wnere < is a bounded estimator gain ooerator constructed so that v(t) remains 
in 0(A). The estimator error, e(t) 2 v(t) - v(t), satisfies 


iSM . (A - K C) e(t), e(0) • -Vjj (3.3) 

It is natural to take the control law 

f(t) * 3 v(t) » G y(t) * 3e(t) 

in olace of (3.1) and hooe for tne best. 

The following theorem gives sufficient conditions ■*or the a.xpcnen- 
tially stable operation of the infinite dimensional controller (3.2) and 
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i c 


(3.4) in closed-loop with the OPS (2.1), (2.3), (2.4): 


■ f ( 

% 

\ ! 


* C 


THEOREM 3.1 

If (A, B) and (A*. C*) are exponentially stabilizable, then there 
exist bounded operators G and K such that the infinite dimensional 
controller produces an exponentially stable closed- loop system 
and v(t) converges exponentially to the OPS state v(t). 

The proof of this theorem uses results from [21] and appears in Appendix I; 
a DPS (A, 3) of the form (2.1) is exponentially stabilizable if a bounded 
gain operator 6 exists which causes the semigroup U.j(t) generated by A + B G 
to satisfy 

t 

! |Ui(t) I ! £ M.J e , t > 0 (3.5) 


where M.j £ 1 and e.j > 0. 

It is known from [22] that, if s » 0 in (2.2), then neither (A, B) 

★ * ★ 

nor (A , C ) can be exponentially stabilizable because 3 and C are finite 

rank (and, hence, compact). This means tnat, if the open-loop DPS is not 
strictly dissipative (no matter how small z is), then no amount of controller 
feedback will make the closed-loop system exponential ly stable. 

The infinite diirensional controller may also be written in convolution 

form: 


f(t) » G U (t - t) K y(-) dT 
0 


( ^ 6 ) 


where U(t' is the semigroup generated by A * 3 G - < C. However, it may 

not be possiole to implement either (3.6) or (3.2) and (3.-) with -'nite 

dimensional devices. Of course, either version of the infinite ci mens i oral 

controller may be approximately implemented and this is the subject of the 

next section. „ 
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4.0 REDUCED-ORDER MODELS QF DPS 


Even though the underlying structure of a OPS nay be described by 
(2.1)-(2.4) and, hence, is infinite dimensional, the controller may have to 
be based on a finite dimensional approximation of the original OPS, in order 
to be implementable in the sense of Section 2.0. Such an approximation or 
reduced-order model (ROM) is a projection (not necessarily orthogonal) of 
the OPS on H onto an appropriate finite dimensional subspace of H. The 
ROM subs pace has dimension N and its projection is denoted by Pj^; the 

residual subspace associated w'th completes the decomposition 
H » projection is denoted by Pj^. The total OPS state v 

is the sum of the ROM state v^^ = Pj^ v and the residual state Vj^ = Pj^ v: 

v(t) « Vj^(t) -r Vj^(t) (4.1) 


Tne choice of subspace and the projection Pj^ completely specifies 
the residual subspace and its corresponding projection Pj^; this choice 
is usually cictated Py the physical aoolication and knowledge (or lack of 
knowledge) about the scecific DPS involved. In many situations, the partial 
differential eouations or the corresponding boundary conditions are too com- 
plicated to permit simple closed-form solutions. Conseouently , an acoroxi- 
mation of DPS, such as the Rayleigh-Ritz-Galerkin or finite element method, 
C23]-C24], must be used to deal with the problem numerical ly; this is one way 
to produce a ROM of the OPS but it is by no means the only way to generate 
such approximations. 

We will restrict our attention in this paper co consideration of 
ROM's obtained from orojection (not necessarily orthogonal) onto a pair of 
reducino subsoaces H„, H, for the differential operator A, i.a., the following 

■ ii K 
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conditions hold: 


H . H„©Hp 

(4.2a) 

Pn (o(p)) 

(4.2b) 

A P„ V . 0 

(4.2c) 

P„ A Pn V . 0 

(4. 2d) 


From [25], p. 172, these conditions are equivalent to the condition that 


commutes with A; 


V e D(A) P|^ V e 0(A) 


A ?n, V » A V 


(4.3a) 

(4.3b) 


Of course, Pj^ also conmutes with A. 

Throughout the rest of this paper, any ROM considered for (2.1) 
will be based on some pair of reducing subsoaces an^ for A . Tne 
projection onto these subspaces may or may not be orthogonal; -^en it is 
orthogonal, 'mq will have and 




I 2 

V ' i + 


. ;2 




ana 


I I p 

i 


! ! a I ! 

I.Pr;! 


1 (4.5) 

An example of reducing subspaces occurs when is a modal subsoace . 


1 .e. , 


» span {3,, . . . , 


wnere A o, = o. for eigenvalues X.. with correspcnaing eigenfunctions ; 
these eigenvalues must be separated from the rest of the scectrum of A 
(i.e., a rectifiable, simple closed curve can be drawn in the comolex plane 
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SO as to enclose an open set containing . . . , Xj,j in its interior and 
the rest of the spectrum of A in its exterior). From [25], p. 178, and 
form a pair of reducing subspaces for A and s ^ ii a finite 

dimensional operator whose spectrum is exactly {X.| X|^} while the 

spectrum of Aj^ s Pj^ A is the rest of the spectrum of A. Note that the 
projections are orthogonal when A is normal [25], p. 277. 

We emphasize that other reducing subspaces may exist for A but 
the modal one is useful when a finite set of separable eigenvalues, and their 
corresponding eigenvectors, can be found; this is the case, for example, when 
A has contact resolvent ([25], p. 187). Also, the modal subspace is physi- 
cally reasonable for producing ROMs of many DPS such as mec.hanically flexible 
structures. It should be noted that, when the exact modes are not immediately 
available (even though they may be known to exist) and approximate modes are 
used for Hj^, the subspaces and are not necessarily reducing subspaces. 

We have the following result for dissipative hyoerbolic OPS: 

ORIGINAL PAGE IS 

THEOREM ^.1 OF POOR QUALITY 


If and are reducing subspaces for A in (2.1), then this 
OPS decomposes into 




(t) , 

Vn( 0) - Vg 

(A. 6) 

3v«(t) 





5t 

' " ’'a * ’ 

, Vj^(O) 

' '’a ■'0 

(4.7) 

■ y(t) » 

^R 

(t) 


(4.3) 

where = 

Pn '^R “ ' R 


^ ^N’ ’ °N '"N 

= C ?,,, 9t 

In addition 

, the projections 

and 

commute with the 

semigroup 


U(t) generated by A and Aj^ generates the Cq semigroup U(t) P^ 
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with the growth property; 

-et 

\\%it)\ \ < e , t > 0 (4.9) 

where = nPplI Mq w1th(MQ, s^as In (2.2). 

The proof of this theorem uses results in [25] and appears in Appendix II. 

The ROM corresponding to the projection of (2.1) onto the reducing 
subspaces and is defined by (4.6) and (4.3) with Cj^ assumed zero: 

Vn,(t) « Aj^, Vj^(t) + Bj, f(t) . v^(0) » Vq 

y(t) » \{t) (4.10) 

or (4.10) may be abbreviated as (Aj^, C^). In particular, when is the 

modal subspaca for N (separated) eigenvalues of A, we call the corresponding 
ROM - the modal ROM . 

The interaction terms Bj^ f and ^ Vj^ are called control and observa- 
tion spillover , respecti vely. It is through these terms that a feedback con- 
troller can interact with the residual subsystem (^.7). This would not be the 
interaction if the ROM had been based on nonreducing subspaces for A; in that 
case model error terms “ 4 v^ and A - wou^ld have appeared in 

(A. 6) and (4.7), respectively . However, these terms are zero when reducing 
subspaces are used for the ROM (as pointed out in Theorem A.l). 

Henceforth, the ROM (Aj^j, Bj^, Cj^^) will be assumed to be controllable 
and observable . It is not difficult to verify whether a given ROM satisfies 
this assumption via the usual finite dimensional space technicues; when the 
modal ROM is used, it is especially easy to verify, e.g., [12]. 

5.0 :MPUg^ENTABL£ Fc-uBAC:< CONTROLLERS "OR DPS 


The performance of the controlled DPS will usually be expressed in 
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terms that can be accomplished with ROM state feedback, e.g., pole placement 
or optimization of a performance Index for the ROM. For example, when the 
modal ROM is used, the modes chosen would be the most cri ti cal (I.e., easily 
excited by whatever disturbance sources are present) modes : increased 
damping and/or stiffening of these inodes would be the main goal of feedback 
control. We will suppose that the desired control law for the OPS (2.1) 

Is given (or adecuately approximated) by 

f*(t) « 6jj Vj^(t) (5.1) 

where the gain operator Gj^ can be selected via pole placen«nt or optimal 
regulator methods, [26], applied to the ROM (4.10). 

If this control law (5.1) could be implemented, the closed-loop 
system (2.1) and (5.1) would have the form (using (4.6)-(4.7) ) : 


VN(t) » (A^ + Gj^) V|^(t) 


(5.2) 


3vp(t) 

“TT-* S ^ ^ 


(5.3) 


where v(t) » 

fcllcw from the result: 


Tne closed-loop stability of (5.2) 


(5.3) would 


;HE0R£M 5.' 


Let the control law be given by (5.1). If B^^j) is controllable, 
then G^ can be chosen so that - e > 0 and 


!v(t)|| < Mq M^ M e 


'o' ■ 


s n 

I w > vj 


(5.4) 


where 

( . 


R ^ 




D I 

■ N' 


if nonorthogcnal projections and are used 

1/2 


’ '‘m > 

if orthogonal projections are used 
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and cu^ 5 S s ||B|^ G^^ll (the control spillover coefficient ), 

^ut 

s G^i generates e with 


'N ''N 


.^N*, 


<Hn 


.-«n‘ 


t > 0 


» *• — 


(5.5) 


and (^^, e) are as In (2.2). 

The proof of this theorem appears In Appendix III. Note that only control 
spillover 3 Is present In (5.4); this Is because direct access to v„(t) has 

n 

been assumed. 

Of course, one rarely has direct access to the ROM states Vj^j(t); at- 
best, one could try to recover them from the outputs of P * N sensors but 
observation spillover would Interfere with this approach. Therefore, the 
control law (5.1) can not be exactly implemented (except in the special case 
that the observation spillover In (4.8) Is zero). Instead, estimates of 
Gj^j Vnj(t) must be generated from the sensor outputs y(t). This can be done 
via the following imolementable, finite dimensional controT^er : 

y(t) H.J 2 2 (t) ;:.£a) 

z(t) « F z(t) * H y(t) + £ f(t), :(0) « 0 (5.5b) 

where dim z ■ S _< N « dim Vj^. 

Let ej^(t) 5 z(t) - T|^ v^(t) and, from (4.6), (4.3) and (5.5b), we 

obtain: 

e^(t) • F e^(t) + H v^(t) 

' (5 . ^) 

e^(0) - Vq 

when the following assumctizns a’^e made about 'sj £. 

('- 3 ) 
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f' 


i * 

f 

I 

t 

t 



I 

I 


i 

ic 


f 



\k win assume that the matrix F can have anv desired stability margin 


ofp» 1 ,e. » 


|e^^l I < Mp e ^ , t ^ 0 


(5.10) 


This (Mans that P can be chosen so that It shares no conmon eigenvalues with 
Ajjj; therefore, by [27], Oiapt. 8, there will exist a unique solution Tj^ of 
(5.8) and hence. E can be chosen to satisfy (5.9). 

From (5.6a) and the definition of ej^(t), we obtain: 


f(t) • Gjj M^it) + H^2 '^R 

when we make the assumption that and H ^2 chosen to satisfy: 

^n.*'N * ^12 '''n ■ ( = -12) 

This last assumption, plus (5.10), can be satisfied when - P<_S <.N, e.g., 
[28], but It can often be satisfied for S < N - P, e.g., [29]; the special 
case S ■ N Is the usual full state estimator or Kalman filter based on the 
ROM, e.g., [30]. Therefore, the assumotions (5. 10) and (5.12) can always be 
satisfied by some implementable controller (5.6) whose dimension S is less 
than N; the actual design of such a controller is a well-studied problem in 
finite dimensional space and will not concern us further; such a controller 
may be synthesized directly with a digital comcuter, [25], Chaot. 6. However, 
please note that the design of the controller (5.6), and all necessary assump- 
tions, can be done entirely In terms of the ROM. This is a reasonable aporoach 
since the residual subsystem parameters (A^, 3;^, C«) may be much less accurately 
known than those for the ROM; in TOSt cases only bounds on these oarameta**s 
would be available in pracf'ce. 

The question remains: How well can this i mo ''emen table controller 

be made to perform? This is the subject of the next section. 


6.0 CLOSED-LOOP OPERATION OF THE IMPLEMENTABLE ROM-BASED CQNTRQLLE.t^ 


Even though the design of the controller (5.6) Is based entirely 
on the ROM (4.10), this controller must operate In closed-loop with the 
actual OPS (4.6)-{4.8). Consequently, the residual subsystem (4.7) Is In- 
volved In this operation via the control and observation spillover terms. 

It Is a rare situation when these spillover terms are zero; this occurs If 
the actuator and sensor Influence functions b^ and c^ lie entirely In the 
ROM subspace Hj^. Most practical DPS problems will Involve a certain amount 
of spillover even though It may be small; this can lead to degraded perform- 
ance and, In some Instances, Instability [30]. 

Let q(t) = Cv]j{t) e^(t)f . From (4.6)-(4.S), (5.7) and (5.11), we 
obtain the closed- loop system ; 


'where 


q(t) • q(t) + 1^2 


3Vj^(t) 

st 

\ 


q('.) ^22 





r. 


-N 



n 


w 


r 


(6.1) 


'12 


®N ^11 **R 
H C, 


^•rj " [3^ ^2'? * ^R * ^R 


H, 


il "R 


with the Initial ccnCiticns -’.j ana 


« -e* 


!'q(t)ll • (i’v^(t)! 


2 


:2sl/2 


Mit have the following useful result: 


f 


If no obsorvation spillover Is present (Cj^ ■ 0) in (6.1), then F 
can be chosen so that > Q and 


1 '»(‘)ll*(l|v^(t)||‘+l|tN(t)lf)’‘< t*’"' n q(0) M . t > 0 (6.3) 


virfjere 


Mp (1 


!VV 


(6.4) 


and a^j), (Mp, Op) are as in (5.5) and (5.10), respectively. 

The proof of this result appears in Appendix IV. Note that (6.4) is not the 
only bound for in fact, the "sharper" we can make the bound on the 
better our later results bec^. 

The next two results are the principal ones in this paper; they give 
sufficient conditions for exponential stability of the closed-loop system in 
term of spillover bounds: 


THEOREM 6.2 


Asscane that (A^, C^) 1s controllable anc observable, that 1,^ 

is chosen in (5.12) so that • c > 0, and that r is chosen 

in (5.6b) so that Jp > If no observation soi Hover is cresert 
(i.e., ■ 0), then, the closed-loop system (6.1) satisfies: 


• -«• 


>(t)’ ' < M, M. M e ' 


• ' 'w I 


, t > C 


15 . 5 / 


where 


I 


IS 


»(t) i Cv^'^U) . 

v(t) • Vj,j(t) -► v,^(t) , 

||w(t)|l • (||v(t)||2 + Wl 

M* Is g1vtn by (5.4) in Thtortm 5.1, and 1$ boundtd as in (6.4) 

in Theoffin 6.1. 

The proof of this thtortm apptars in Apptndix V. 


THEOREM 6.3 


Ltt tht optn-loop systtm (2.1) bt a strictly dissioativt hyperbolic 
OPS (i.t., c > 0). Undtr tht samt assumptions as Thtortm 6.2 txctpt 
that obstryation spiUoytr is prtstnt . tht clostd-lcop systtm (5.1) 
stability margin is at Itast c ^whtrt 


M- M* r 


c "0 


( 6 . 6 ) 


whtn r s M3'' CqM Cj^M 'S su'^ficiently small £ > 0. 

In that case, 


I ' 


w(t)M < M* t‘” ' . 


t > 


(6.7) 


'C "0 ■ ■ ' 'o' 

Tht proof of this theorem scpea.'‘S in Apoendix 71. There is no ooint In ccn- 


sidering such a result whtn (2.1) is not strictly dissioativt hyperbolic^bt- 
cause t could not eytr by positiyt in (6.6). 


7.3 APPLICATION: CQNTPQl 0^ ~l£XI3L- $T"UC~’JF.£S 

"eecbacx control of vibrations in mechanically fle.-'.ible structures 
nas ICO li cation to many cur*^nt engineering oroolems. Control of flex'! ole 
structures, in general, was addressed in [30], [Cu] mo [38]. A new aoolica- 
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tion arta - control of large aerospace structures - Is very active at present; 
the Current control theory trends in this area have been surveyed In [12]. 

The results obtained In this paper have direct application to control 
of flexible structures, In general, and largo aerospace structures, In parti- 
cular; this application will be discussed In this section In order to Illustrate 
the use of the results of Sec. 6.0. This section will follow the format 
established In [38]; for clarity's sake, some of the format will be repeated 
here. 


We consider the class of flexible systems that can be described by a 
generalized wave equation; this represents the Idealization of many mechanically 
flexible structures. The generalized wave euqatlon Is given by: 


3^u 


♦ bl^ ♦ A„u 

at 0 


( 6 . 1 ) 


which relates the vector of displacements u(x,t) of a body a bounded open 
set with smooth boundary 3H In n-d1mens1onal Euclidean space s”, to the applied 
control forces ?(x,t). The operator is a ti Me- Invariant, symmetric d-f-'er- 
entia'i operator with compact resolvent and lower semi bounded soectrum. "he 
domain 0(4^) of is dense in the Hilbert space L^C.) with (•,•) benot^ng the 
usual Inner product and ' ' • denoting tne associated norm. ~ne natural camping 
in the structure Is modeled by tne term where o is a non-negative constant. 
The control forces 


F( 


X# w ' 


0 ' V ^ / 




!■ I 


b^(x) 



PAtrh' Ih 

OP POOR QUA^JTv (5.2 


are provided by M actuators with influence 
are measured by = a''eraging sensors 
y(t) ■ CgU(x.t) 

•where y,-(t) ■ c,(x;u(x,t)dx 'with j ■ 
sensor functions b^x), c^-x) are in L*(n) 

I _ 


funcficns b^(x). ’"ne dso*. acements 


1,2,...,“. "^ne actuator »rd 
and noma li zee to nave unit integral. 
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'Alhen the support of b^(x) is in a small neighborhood of a point , we say it 
is a point actuator and, similarly, we define a point sensor . The point 
actuator and point sensor situation is of special interest here. This class 
of distributed parameter systems includes interior and boundary control of 
vibrating strings, membranes, thin beams and thin plates. Although, only 
displacenrent sensors are considered in (6.3), this is not a restriction - more 
general sensors may be nradeled. 

It is well known ([25], p. 277) that the spectrum of A contains only 
isolated eigenvalues with corresponding eigenfunctions ■pj^ such that 
^ ^ ... 

and Ai,^ = We will assume that is positive. The resonant mode fre- 

guenci_es, of the structure are given by = (Xj^) and the corresponding 
eigenfunctions are the mode shapes . Thus A satisfies 

(A^u,u) > a| [u! 1^, a > 0 (5.4) 

1/2 2 

and has a square root A^ . Every vector usL (n) has a unique representation 


u(x) » : u.3.(x) 


k=1 


K < 


ra S'! 


R 


where u. = .Vuo.dx and we define the orthogonal orojections p.' , ?_' by 

N 

^ u. ?. 


a \ . 

”0 “ 1^;, “k'k 


( 6 . 5 ) 


p u » Z u^i, 
0 k*N+l ^ 


Let V be the domain of A, and W be the domain of A new ocerator A is 

0 0 

defined in H by 

a) D(A) * V X W i H, 


b) A 


i ■ [••- • v] 


( 6 . 7 ) 


for ueV, wcW 
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'he energy inner product (’.Og defined on by 


I ^ 


5 Ag'/?U2> + {w,. Wjl 


(5.8) 


for u^, Ug E V and w^, v <2 e W , and the Hilbert space H is defined as the 
closure of in this energy inner product. The associated energy norm is 
denoted by Ij-jlg and is a measure of the total potential and kinetic energy 
in (u,-— ) where u is a solution of (6.1). Let v » [u^, |^] be in H and 
write (6.1), (6.2) and (6.3) as 

'll- Av.Bf . v^eH^ 

(6.9) 

y » Cv 

where B • 5 and C • [C, 0]; titis 1s in the form (2.’). 

L'oJ ° 

It is easy to see that 


K I 

• I 


and I i V M 


7 • 2 

[\. u. - + u. ]. Also, A generates a C semigroup 'J(t) 

K K % 'J 


iSy [1-^3 P> il"59. The following result, due to Goldstein and Sosencrans 
[39], gives conditions under which this semigrouo U(t) satisfies (2.2) and 
hence (6.9) is a (strictly) dissipative hyperbolic system; 
if b < 2a then 


' I I 


< Y : 'v(O) 


, t > 0 


where y* [^(b/Ga) ][l-(b/2a)]“ ‘ . Conse-buently , U(t) satisfies (2.2) w’th 


Mq - (y)^'^^ and 


b/2. ^lote that the hyoothesis for this result means 


that £ < a, I.e. the damoinq must be smaller than the souare c* the lowest 
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mode frequency . Another bound of the type (6.11) Is available in [36] without 
the restriction £ < a; however, the damping coefficient b is not quite so 
directly related to z in (2.2). The bound used above will be very convenient 
to illustrate our results. 

Recalling (6.6), we define the projections Pf|, on H by 





u 


Po 

.W. 


r R 1 

'u'' 


Po ^ 

_w_ 


- Po’^'" - 






( 6 . 12 ) 


and note that they are orthogonal in the energy inner product and that Pj^^.4 and 
Hjj»PpH form reducing subspaces of H, as described in Sec. 4.0. In fact, is a modal 
subspace, and Vj^“Pj^v and Vj^'Pj^v satisfy (4.1), (4.4)-(4.5). Consequently, the 
system (6.9) partitions into the form (4.6)-(4.3)^and the corresponding modal 
ROM is given by (4.10) where A^, C^j may be identified with the matrices 


L’"N 


•bl 


M 



1 

o 

1 


n Q 

m mm 


[C^° 0] , respectively, with * diag [\.] 




3.^ = [(3^.3)^)]' . and = L(c,,o,^)]. Thus, "he ccntrol labi li ty and 
observability of the ROM can be determined from (.\j^, 3,^“', wnen the damping 
b is small ([38], Theorem 2.1). Of course, it is quite easy to assess the con- 
trollability and observability of this latter system and the conditions deoend 
on the interaction of the actuator-sensor influence funczions wicn the mode 
shaoes o, , . . . , o^. 

If the ROM is control lacle and ooservable, the usual ~inite I'T-ensicna’ 


state space controller can be designed based on the modal ROM (e.g. [26]]. 
has the form: 
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f(t) * V|^(t) (6.13) 

Vnj(t) • Aj^ Vj^(t) + f(t) + Kj^j(y(t) - y(t)) 

(6.14) 

y(t) • Cjj Vfj(t) , Vjj,(0) ■ 0 

This is a particular case of (5.6) with dimension S«2N (twice the number of 
modes in the ROM) and parameters: H.j.j«0, F * - K|y|Cj^i and 

Also, note that e,^(t) • Vjjj(t) - V|^(t), i.e. » l 2 jyj. Since ROM is 

controllable and observable, Gj^ and can be chosen to achieve (5.5) and 
(5.10) with any desired stability margins and "p. 

Thus, (6.13)-(6.14) represents an imolementable (finite dimensional) 
modal controller for the flexible structure (5. 1 )-(6.3), or equivalently (6.9). 
It is based on a modal RGH, i.e. it is designed to control the first N modes 
of the structure (of course, different sets of N modes could have been used - 
it is not necessary that they be the first N modes). The basic assumption 
made in synthesizing such a modal controller is that the ^ modes chosen for 
control will adeouatelv ’•eoresent the important vibrational behavior of the 
whole flexible structure ; this assumption can often be satis'ied in oractice, 
and it is the basis for many structure controllers. 

Having synthesized an implemen table modal controller to adequately 
tailor the vibrational response of M critical modes in the flexible structure 
we must not think we are done. The closed-loco benavior of our controller is 
affected bv the resicual modes in H, for which we have not designed, as well 
as the controlled modes in for which we have designed. These interactions 
occur tnrcugn the control and observaticn soil'over terms de'^ined in Sec. -.1 
and can drastically altar the ces'red oerformance of the controlled structure 
in fact, an example is presented in [30] where the closed-loop system become'', 
unstable. 
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Consequently, the analysis of closed-loop behavior discussed in 
Sec. 6.0 is especially pertinent. The principal results Theorems 6. 2-6. 3 
give an indication of how much spillover can be tolerated in the closed- 
loop system. In the context of flexible structures. Theorem 6.2 may be 
compared with th» results of [38] where observation spillover was assianed 
negligible (however, note that no damping was present in [38]). Combining 
the results of Theorems 6.2 and 6.3 for this application, we find that the 
closed-looped system, consisting of the flexible structure (6.1)-(5.3) and 
the finite dimensional modal controller (6.13)-(6.U)jWill remain ex- 
ponentially stable with stability margin t ^ given by (6,6), if the observa- 

/ 

tion spillover coefficient F is sufficiently small to make e >0, or 
equivalently: 

r < (6.15) 

V* w 

In this application, £ = j where b<2a in (6.1) and (6.4), is given by 


^0 = 


2a^ 

2a-b 


1/2 


(6.16) 


Also, from (5.4) since orthogonal projections are used: 


M* . (1 . ^ -n’> 


1/2 


(6.17' 


where and -£, and, from (6.4), 


ib'gj 

^ -Vr * 4-^) 


(6. IS) 


where 4 » t- - c., and i\, come from (5.5) and (5. 1C) resoectivelv. 
The ''single bar" norm denotes the matrix norm comcatible witn the euclidean 
norm. Finally, tne control and observation spillover coe""icient3 S arc 7 
are given by 

,\l/2 


3 • <( : : (5.,:,!-) 

^ ^ E \i-l loN*) ‘ / 


\C . j 


ORWENAL PAiJh' ; 
OP POOR QUALITY 
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iVoii s 

E ^ \j-l k«N+l ^ / 


( 6 . 20 ) 


It may seem like quite a complicated business to obtain all the data 
to check (6.15); however, most of the information can be obtained from the 
finite dimensional ROM and its corresponding controller, and the rest comes 
from the structure properties a,b. The spillover coefficients (6.19)-(6.20) 
depend on the controller gains Gj^, Kj^ and the actuator-sensor influence on 
the residual modes. Consequently, (6.15) gives an indication of how much 
spillover can be tolerated in the closed-loop system. 

Additionally, (6.15) gives the control system designer an aid to 
redesign the system for better overall stability; Quite simply, either the 
controller gains must be reduced or the actuator-sensor influence functions 
b., Cj must be chosen to minimize their influence on the residual modes. 

Of course, this latter could be accomplished if one had completely free 
choice of b. , c,. ; one would choose them to be in which is orthogonal 

to and both 7 and 3 would be zero. Please note, for a flexible structure 
this would mean the influence ^unctions were linear ccmoinations cf tne con- 
trolled modes only. Yet, in practice, most actuators and sensors nave a 
very localized influence (they are nearly opint devices), and this oractical 
constraint would make it very unlikely that zero soil lover could be acnieved. 
Nonetheless, adjustments in the location o* actuators anc sensors could be 
made in order to reduce the spillover. Further, one could reduce soillover by 
controlling TOre modes at the cost of increasing tne ccntrcller Comdex- ty 
and corraspono'i ng computational bu'"cen on one cn-'-ne comcuta'-. Itne- 
approaches to spillover reduction ana compensation (e.g. orefiltering) are 
surveyed in [12]. The aprior-. bound (5.151 can be used to assess the trade- 

<k 

cffs in these approacnes for -lexible structure control. 
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3.0 CONCLUSIONS 

In order to better understand the behavior of highly oscillatory dis- 
tributed parameter phenomena such as mechanically flexible structures, the 
class of dissipative hyperbolic distributed parameter systems (OPS) has been 
Introduced. The fundamental property of such OPS is that the open-loop 
system (2.1) must be exponentially stable, i.e., satisfy (2.2) for e > 0. 
Conditions have been presented in Section 3 for exponentially stable opera- 
tion of the infinite dinwnsional controller, (3.2) and (3.4) or (3.6), in 
closed-loop with the OPS; however, the principal concern of this paper has 
been with finite dimensional controllers for the DPS. Such controllers are 
practical in the sense that they can be imoleirented with a small number of 
actuators and sensors and one^or more, small on-line computers. 

The finite dimensional controllers (5.6) developed here in Section 
5 are based on reduced-order models (ROMs), i.e., finite dimensional approxi- 
mations, of the actual OPS. The class of ROMs used here has been derived 
in Section 4 from reducing subspaces; a large subclass of reducing subsoaces 
contains the modal subspaces which have been used extensively in flexible 
structure modeling and control. The controller synthesis is done entirely 
in terms of the ROM. 


Since the finite dimensional controllers developed in Section 5 
must operate in closed- loop with the actual DPS and not just the ROM, it is 
crucial to analyze the spillover , i.e., controller interaction with the 
residual (unmodelled) part of the DPS. One cannot expect to be able to :o 
this in utnx3St detail since the resicuals of the DPS are the part tnat is 
least <ncwn; however, spillover bounds can be prescribed for stable oceration 
of the closed- loo? system and this has been done in Section 5. The principal 
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results are Theorems 6.2 - 6.3 ; the value of these results is that the 
stability conditions may be checked using only knowledge of the ROM and 
bounds, 3 and r, on the spillover present in the open-loop DPS. The bounds 
3 and r may be calculated directly, or estimated where necessary, from the 
ROM projections of the actuator and sensor influence functions b^ and Cy 
When the spillover bounds are small enough to satisfy the stability conditions, 
the finite dimensional controller will function to enhance the behavior of 
the ROM subsystem of the OPS without causing any instability in the residual 
subsystem. When the stability conditions are not satisfied, the system may 
be stable or it may not; however, this should be a warning that some spill- 
over comoensation be added and the conditions rechecked. 

Several closely related OPS topics seem to require further in- 
vestigation: 

(1) sharper spillover bounds 

(2) bounds obtained from singular perturbations 

(3) nonreducing subspacss and model error 
(a) non! inear control 

(5) adaptive control 

Since the spillover bounds have been octained from regular perturba- 
tions, they may be sharpened by making tighter estimates, e.g., for in 
(6.4); connections need to be made between our results and ot.her stability 
results for large scale and DPS obtained by Siijak [31] and Michel and 
Miller [32]. In addition, bounds obtained 'or DPS via singular perturbations 
look promising for many applications, e.g., [33]. 

Many ROMs are orojections onto nonreducing sucspaces; for e.xamcle, 
when approximate modes are substituted for exact modes, the approximate modai 
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subspaces are not reducing subspaces. In the case of nonreducing subspaces, 
the controller interaction with the residual subsystem can occur through non- 
zero model error terms k 9^^ and Aj^j^ * ^ ^N’ ** spill- 

over terms. The bounds obtained here for spillover should be extended to 
model error bounds, as well. 

Finally, preliminary results in nonlinear and bilinear control of 
OPS have been obtained in C34]-[36] and the surface has barely been scratched 
in adaptive (or self-tuning) control of OPS, e.g., the survey [37]. Much 
more remains to be done in both of these areas. 
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Sinct (A, B) and (A , C ) art axponentlally stabiUzablt, thert 


♦ ** 


•xist bounded G and H such that A-^BGandA H gantrati U^(t) and 
V(t}» respectively; U-jCt) satisfies (3.5} and V(t) satisfies 


ll»(t)ll < »2 t *** , t > 0 


( 1 . 1 ) 


where M 2 > 1 and £2 > 0. Choose K ■ -H and from [21], p. 251, A - C gene- 

rates U 2 (t) where U 2 (t) ■ V*(t); consequently, !'U 2 lt)i! ■ MV (t)|| ■ !lV(t)!|. 


Therefore, 


-C2t 


il«(t)!l < M 2 e " Mvgll 


Also, from (3.3}-(3.4) and (2.1), 


■ (A + 3G) v(t) ♦ BG e(t) 


or, from (3.5) and (1.2), 


!!v(t)!( < M^ e ' (T 

'/«nere 1 • z, - Let M.. = M, M„ and : 2 min(£,, £,}. 

1 i t I 6 • 2.,. 

If .1 > 0 . iXt)’: < M. e‘-^ i:v (1 * I’bg;; ■ - t— ). 


(1.2) 


(1.3) 


(I.A) 


'/«nere 


If i < C. 


: V(t): ' < M, ] (1 ^ ' 'BG: 

^ m O' 

In either case, vor t ^ 0, 


1 - e-“« 


i!v(t)|i < M. e"-^ IjVg!! (1 ^ !'BG|i/lAl) 

From ( 1.2) , for t >_ 0, 

' :«(t}" < M, e‘-‘ ! 'v ' ’ 

— ' e ■ 


(1.5) 


( 1 . 5 ) 


"here’cre, the closec-loco system *3 exrcr.er,t*al'y stacie. 
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APPENDIX II; P«^QOF OF THEOREM A.1 


By th« definition of reducing subspaces (^.2c and d) , s P^ A Pj^ ■ 0 
and Aj^ 1 Pj^ A Pj^ ■ 0 and (2.1) becomes (4.6)-(4.8). Since the definition 
of reducl’ig subspaces Is equivalent to Pg,^ and P|^ cotimutlng with A (see 4.3b), 
then, by [25], p. 173, they also cosmute with the resolvent operator R(X, A) 
of A for any X In the resolvent set; because A generates a semigroup, the 
resolvent Is non-empty. But R(X, A) Is the (unique) Laplace transform of the 
transform of the semigroup U(t), e.g. [25], p. 482; therefore, the projections 
imjst cc.i«mite with U(t). 

Let Uj^(t) 2 Pj^ U(t) Pf^ then U^(t) Is a semigroup on because 
Uj^(t) Uj^(t') • P,^U(t) U(t') pJ^ - U(t) U(t') ’ ?p 'J(t * f) P-, • 

Uj^ft + t‘) for t, t' >0; this uses the fact that P^ commutes with U(t). Of 
course, Uj^(O) ■ P^ which Is the Identify on and the strong continuity follows 
from that of U(t). It Is also clear rh;jt A« generates U-(t) because, for v in 

It '' . . 


'r 


Ujj(t) V - V V - V - 11m U(t'' V - V _ 

• 1^15 1 t-0 1 ^ t-0 ^ 

Pg^ A V ■ by the continuity of Pp. "hp growth orooerty (-1.3) *ol'ows 

from (2.2) and : ;U (t)! 

U(t). 


■?„ L'(t)" < "U(ti 


because P^ cermutes wl tn 


ORIGINAL FACE IJj 
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APPENDIX III: PROOF OF THEOREM 5.1 


Since (Aj^, Bj^j) is controllable and finite dimensional, can be 
chosen so that > e. From (5.2)-(5.3) and (5.3), 

<% e ^ l|v„(0)|| (III.l) 

IIVt)|l < (llv^il * / e" llBjG^II ||v„(t) i| dr) (III. 2) 

because Uj^(t) V(^(0) = Pj^ U(t) from Theorem 4.1. Combining III.l and 
III. 2, obtain 

l!v„(t)ii < (1 (in. 3) 

Therefore, using (III.l) and (III. 3), 


-St 


liv(t)!l < !|vj^(t)|l + l|vj^(t)!| < e 

[M„ t M^(e + M,(i - .) '■ ij,) n^ii (i!i-«) 

Note that (III. 4) is the same as Triggiani's bound in [17], Theorem 6.1 




From (III.*!), obtain (5.4) for nonorthoconal projections because e 

m \ ^ • A ^ 

\ ’-N (' - s '''") i e i I " Pfl PPd 


-M • 


> 1 from 


^ = ( : •V|^(0) I ! + \ i |Vj^(0) i ;) and ; !v{-) ; ‘ ! v,j(t) 


><-, a 


Theorem 4.1. However, when the projections are orthogonal, then ''vj,(t)|' ^ 

7 ,, , , . 

< (M e””^)^ [■/il.. /rt\ M J. I ! .. /n\'lN2 i M.. tn\‘.'.2. 

<0 - :i.^h (Yn «''')^ [l!^(CI)il^ ' !!v,(0)“^’ 

Mv{t)il < (1 ‘ 3j, * 

2 ^ 7 2 2 

because (a a b) a" £ (1 + a ■*■ :i") (a r b") which isLemma 3.1 in [19j. 


Vj^(O) I ! + 

j ! iVj^(o) 

,«N [i 




MM,, e”-' i ’ 
0 11 

< 

o 
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APPENDIX IV: PROOF OF THEOREM 6.1 


With Cp - 0 in (6.1) 
q(t) - e 


From (5.5) and (5.10), it is easy to obtain 
l|eN(t)|| <Mp e"'"' ||T^ V|,{0)|| 


and 


“Cut t 

V'lii * 

1 ^ 12 ! ' ' ) i 


(IV. 1) 


(IV. 2) 


(IV. 3) 


but, using (IV. 2), (IV. 3) becomes 

<s^t i !Bu H,,! i I !T.! I 

® [1 ^ 3||v^(0)i| (IV. 4) 

0 

because F can be chosen so that A > 0. From (IV. 2) and (IV.4), the desired 

V# 

result follows because and are not less than unity and A. > 0. 
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APPENDIX V: PROOF QP THEOREM 6.2 


Let Cf^ ■ 0 in (6.1) and note that I .|2 
(6.1), obtain 


fq(t) » A.,^ q(t) 


3V|^(t) 

at" ' 


^21 ^ Ar V(^(t) 


0 and A 22 ■ Aj^ 


From 


(V.l) 


where liq(t)j! satisfies (6.3) in Theorem 6.1 by aopropriata choices of F 
and Gj^. From Theorem 4.1, Aj^ generates the semigroup Up^(t) with growth 
property (4.9). Using this in (V.l) produces 






.-CT I ,T 


A 2 ^" !!q(t)M dt ) (V.2) 


where j !l 2 il 1 • • 

Since (6.3) is similar to (III.l) and (V.2) is similar to (III. 2), 
the proof of this theorem - Theorem 6.2 - can be carried out in the same 
manner as the oroof of Theorem 5.1 in Apoendix III (where a(t) ta<es the 
place of Vj^(t) and replaces M^.) 
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APPENDIX VI: PROOF Q F THEOREM 6.3 

Let w(t) - Cq*'’(t) 


w(t) ■ w(t) + '^(t) 


(VI. 1) 


where ■ 


h\ ° 
^21 


and AAg » 


u m^2 
0 


From Theorem 6.2. generates a semigroup Sg(t) with growth property, from 
(6.5); 


ISgCt)!! < Mg e 


-£t 


(VI. 2) 


where Mg = Mj.MqM*. Also, from (6.1). 


1/2 




(VI. 3 ) 


Since is bounded. generates a semigroup S(t) (e.g.. [13] ?. 30 

or [6] p. 210); therefore 

w(t) = S(t) w(o) 


(VI. A) 


However, *rcm (VI.l), 

w(o) » S ft) w(o) r S,(t-t) aJ w(t) dT 

U I'h Mr- 

W 

"Taking bounds in (VI. 3) and using (VI. 3), obtain 

t .. 

"w(t) ' ' < Mp e"”" [| :w(o) i I ^ e’“ ' x(") ' ! dr J 

- t 0 

Let n(:) = e^^ ''w(t) and, *rcm (7!. 5), obtain 

-!( t! < M.. n(o ‘ * 7.' h(-) dt 




(VI. 6) 




From the Oronwall Ineouality (e.g., [6j p. 12A or [Id] p. II. -As) apolied .o 
(VI. 7). 
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h(t) < Mg e‘^ h(o) (VI. 8) 

or, equivalently, 

llw(t)i! < e‘^'^ |!w(o)|| (VI. 9) 

where e' is given by (6.6) and this is the desired result (6.7). Clearly, 
the semigroup S(t) is exponentially stable when t is positive. 
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ABSTRACT 



Large flexible aerospace structures, such as the solar power satel- 
lite, are distributed parameter systems with very complex continuum des- 
criptions. This paper investigates the use of finite element methods 
to produce reduced-order models and finite dimensional feedback con- 
trollers for these structures. Our main results (Section 5.0) give con- 
ditions under which stable control of the finite element model will 
produce stable control of the actual structure. 


r 




t 
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1 .0 Introduction 

Large flexible aerospace structures, such as the solar power satel- 
lite^^ are distributed parameter systems with very complex continuum 


descriptions. Often, their dynamics are 


modeled by large-scale 


systems of ordinary differential equations which are obtained via the 
finite element method for structures. To obtain implementable, i.e. 
finite dimensional, feedback controllers for these structures, we shall 
use this model as the basis for the control synthesis. 

The question of closed-loop operation of such controllers within 
the actual structure is certainly a fundamental one. In ‘this paper, 
we shall point out some connections between the finite element models 
and the closed-loop stability. These results will be used to give an 
answer to the following version of the above fundamental question: 

When does stable control of the finite element model ensure stable con- 
trol of the actual structure? 

"or distributed parameter systems, in general, the finite element 
method (FEK) for time-dependent problems is very well explained in [2] 
Chaoter 7, [3] Chapter 9, and the excellent survey [4]. The related 
approaches called "spectral methods" are described in [5]. The use of 
FEM for analysis of distributed parameter control problems involving 
interior and boundary controls has been considered, for example, in [6] 
[10]; all of these references are concerned with parabolic (i.e. diffu- 
sion-like) distributed systems and, with the exception of [10], they 
are concerned with open-loop (i.e. non-feedback) control. 

In contrast, our concern is feedback control of hyperbolic (i.e. 
highly oscillatory) distributed parameter systems and the use of FEM 
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to develop low-order finite dimensional controllers for such systems. 
Control approaches for structures have been surveyed in [11]-[13]. The 
most popular approach for structural control is the resolution of the 
structural vibration into normal modes, (semi) discretization of the 
distributed system by truncation of this modal expansion, and controller 
design based on the truncated modal model; this is called trodal control . 
Modal analysis of structures has been described in many references but 
an enjoyable survey is [14]. The use of modal analysis in feedback 
controller design depends on a complete knowledge of the modal data 
which is rarely available for something as complex as a Targe aerospace 
structure. Instead, this modal data is approximated by finite element 
structural analysis computer codes (e.g. NASTRAN), and it is this approx 
imate mooal model that is used for the controller design, e.g. [15]. 

The stability of closed-loop operation of the finite dimensional 
modal controller in the infinite dimensional distributed parameter 
structure is questionable even when the actual modes are known, as shown 
by a simple example in [16]. The mechanism for instability is the 
combination of control and observation spillover in the residual (un- 
modeled) modes of the structure [13]. However, closed-loop stability 
bounds have been obtained for highly oscillatory distributed parameter 
systems with low-order finite dimensional controllers [17]. These 
bounds indicate the tolerable levels of spillover for stable operation, 
but they are based on "reducing" subspaces. For large aerospace struc- 
tures, this corresponds to exact knowledge of the modal data used in 
the controller, although the residual modal data need not be so pre- 
cisely known. This paper addresses the closed-loop stability question 


» 


It 


# 


r 


for structure controllers based on finite element approximation of the 
modal data. The distributed parameter control problem for flexible 
structures is described in Section 2 and the finite element approximation 
in Section 3. The corresponding finite element modal controllers are 
presented in Section 4. and our main results on the closed-loop stability 
are Theorems 1 and 2 in Section 5. 


2.0 The Distributed Parameter Control Problem for Flexible Structures 

We consider the class of flexible structures that can be described 
by a generalized wave equation; this represents an idealization of many 
mechanically flexible structures. The generalized wave equation is 
given by: 




( 2 . 1 ) 


which relates the(vector of) displacements u(x,t) of a body Ci, a bounded 

open set with smooth boundary Sfi in n-dimensional Euclidean space r”, 

to the applied control forces F(x,t). The operator is a time-invariant, 

self-adjoint differential operator with compact resolvent and lower semi- 

bounded spectrum. The domain D(A^) of A^ is dense in the Hilbert space 

L^(n) with {’.Oq denoting the usual inner product and I1-!1q denoting 

the associated norm. The natural damping in the structure is modeled by 

the term b I 7 where b could be any bounded self-adjoint operator commuting 

with A^ but is taken here to be a non-negative constant, for convenience. 

0 

The control forces 


F(x.t) « B f(t) * I b.(x)f.(t) 

” i«l 

are provided by M actuators with influence functions b^(x). The dis- 
placements are measured by P averagit>g sensors 


( 2 . 2 ) 


r 


4 


y(t) - CpU(x.t) (2.3) 

where yj(t) ■ Cj(x)u(x,t)dx with j • 1, 2, . . . , P. The actuator 
and sensor functions b^(x), Cj(x) are in L (n) and normalized to have unit 
integral. When the support of b^(x) is in a small neighborhood of a point 
x^, we say it is a point actuator and, similarly, we define a point 
sensor . The point actuator and point sensor situation is our principal 
interest here. This class of distributed parameter systems includes 
interior and boundary control of vibrating strings, membranes, thin beams 
and thin plates. Although, only displacen«nt sensors are considered in 
(2.3), this is not a restriciton - more general sensors may be modeled. 

It is well known ([18], p. 277) that the spectrum of contains 
only isolated eigenvalues with corresponding eigenfunctions such 
that 


and A^C|^ * ^'k'k' positive. The resonant 

mode frequencies u;. of the structure are given by ^ 

corresponding eigenfunctions are the mode shapes . Thus A^ satisfies 

(AoU,u)o > a! iu! , a > 0 (2.4) 

1/2 2 

and has a square root A^ . Every vector ucL (^l) has a unique represen- 
tation 


u(x) « : 


k-1 


(2.5) 


N „ R 


where Uj^ ■ define the orthogonal projections , P^ by 


n 


N 


V " “ = “k-k 

0 k-1 


( 2 . 6 ) 




“k»k 


k-N+l 
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Let V be the domain of and W be the domain of 
A is defined in H by 

a) 0(A) • V X W s 

b) A^“j • ^-bw - A^uj for ucV, wcW 

The energy inner product (•,•) is defined on by 


1/2 


A new operator 


(2.7) 


( 



’TJ I '^2' 

for U.J, Ug c V and w^ , Wg c W, and the Hilbert space H, called the state 
space , is defined as the closure ov in this energy inner product. The 
associated energy no rm is denoted by t|-ji and is a measure of the total 
potential and kinetic energy in (u,!^) where u is a solution of (2.1). 


) 5 (A^ A^ ^^“2^0 * ^'^1* '^2^0 


( 2 . 8 ) 


T T 


St' 


Let V ■ [u"^ , ] be in H and write (2.1), (2.2) and (2.3) as 



(2.9) 


y^ere B ■ 


» m 

g and C • [Cq 0]. 


It is easy to see that 



( 2 . 10 ) 


and Ijvll^ 


I 

k-1 




Also, A generates a semigroup U(t) 


hy [19], p. 11-59. The following result, due to Goldstein and Rosencrans 
[20], gives conditions under which this semigroup U(t) is (slightly) 
exponentially stable and hence (2.9) is a (strictly) dissipative hyperbolic 
system: 
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if b < 2a then 


v(t)ll < ||v(0) 


(2.11) 


t > 0 


» ^ — 


v^re c ■ b/2, and y • [T‘*-(b/2a)][l-(b/2a)]‘^ . Note that 

the hypothesis for this result means that e < a» I.e. the derolno must 
be smaller than the square of the lowest mode frequency . Another bound 
of the type (2.11) Is available in [21] without the restriction c < a; 
however, the dancing coefficient b is not quite so directly related to 
c in (2.11). The bound used above will be very convenient to illustrate 
our results. 

Recalling (2.6), we define the projections Pj^, P^ on H by 


P, 


p ■ 

u 

' 

p 1 


2 

j 

, W . 


.Po'-J 


'Z) 


w 


'’o"“ 


P ^w 
L 0 


and note that they are orthogonal in the energy inner product and that 
^N* reducing subspaces H, as described in [17]. 

In fact, is a modal siA)space, and Vj^® P,^v and Vj^« P^^v satisfy: 

3Vj^ 

^ ''n ®N ^ 


3t 


^•R ''r " ®R ^ 


(2.13) 


''n * ^R ''r 




'‘4 


i'f C 

N 


it ^ 


r 


r 


r 


I 



» 


whert Aj^ “ P|^ A P|^, B|^ ■ P|^ B, Cj^ ■ C Pj^, and, corrtsoondlngly , for 
(Aj^, Sj^, Cjj). Note that Aj^, Bj^, Cj^ may be Identified with the matrices 

. [Cf^° 0]. respectively, with ■ dlag Cx^, 

. . . , Xj^], Bj^° ■ C\bj,o,^)]'’’ , and C,^® ■ [{c^.c,^)], and, If point sen- 
sors and actuators are used, then B,^® and Cj^® involve only mode shapes 
evaluated at the actuator and sensor locations. The terms Bp f and 
Cp Vp are called, respectively, control and observation spillover , and 
It Is through these terms that the residual modes Vp are affected 
by any feedback controller which accepts the sensor outputs y and 
generates control commands f. 

Modal controller? can be based on the reduced-order model (ROM) 
whose states are Vj^ and whose parameters are (Aj^, Bj^, Cj^); such a model 
is obtained by Ignoring (for the purpose of controller synthesis) the 
residual modes, i.e. setting 3p “ 0 end C- • 0. The conf'ol labil ity 
and observability of the ROM can be determined from (Aj^, Cj^®, Cj^®) when 
the damping b is small ([22], Theorem 2.1). Of course, it is quite 
easy te assess the controllability and observability of this latter 
system, and the conditions depend on the interaction of the actuator- 
sensor Influence functions with the mode shapes c.| , . . . , $,j. 

If the ROM is controllable and observable, the usual f inite dime n- 
sional . state space controller can be designed based on the mooal ROM 
(e.g. [23] or [24]). It has the form: 




0 


In 

-bl 


N 


f'(t) » Vj^(t) 

r V,^(t) « Aj^ Vj^(t) ♦ Bj^ f(t) ♦ Kj^(y(t) - y(t)l 
(. y(t) ■ Vj^(t) , Vjj(O) ' 0 


(2.12) 


(2.14) 
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This controller has dimension 2N (twice the number of modes in the ROM). 
Since the ROM is controllable and observable, and K,, can be chosen to 
achieve any desired stability margins Oj^ and Op for 6|^ and 

" Si *'N >^spectively; optimal control, as well as pole placement 
techniques, may be used successfully for the controller design. 

Thus, (2.13)-(2.14) represents an implementable (finite dimensional) 
modal controller for the flexible structure (2.1)-(2.3), or equivalently 
(2.9). It is based on a modal ROM, i.e. it is designed to control the 
first N modes of the structure (of course, different sets of N modes 
could have been used - it is not necessary that they be the first N 
modes). The basic assumption made in synthesizing such a modal con- 
troller is that the M modes chosen for control will adequately represent 
tha important vibrational behavior of the whole flexible structure ; this 
assumption can often be satisfied in practice, and it is the basis for 
most structure control lars. 

Having synthesized an implementaole modal controller to adequately 
tailor the vibrational response of N critical modes in the flexible 
structure, we must not think w'.* are done. The closed-loop behavior of 
our controller is affected by the residual modes in Hp for which we have 
not designed, as well as the controlled modes in for which we have 
(. ’Signed. These interactions occur throurh the control and observation 
spillover terms defined in C2.13) and can drastically alter the desired 
performance of the controlled structure; in fact, an example is presented 
in [16] where the closed-loop system becomes unstable. Consequently, 
the closed-loop stability analysis presented in [17] is extremely impor- 
tant; these results give an indication of the amount of spillover that 
can be tolerated in the closed-loop system without loss of stability. 
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3.0 The Finite Element Approximation and Open-Loop Convergence 

Unfortunately, all of the controller synthesis and closed-loop 
analysis presented so far is based on reducing subspaces and and, 
in particular, on exact knowledge of the system modes 4>i , . . . , <l>|^ 
used in the ROM; this modal data must be approximated in most cases 
due to the complexity of the actual partial differential equations and 
associated boundary conditions which describe the structure (2.1). Such 
approximations do not lead to reducing subspaces, in general; therefore, 
further discussion is necessary. 

The most natural approach (and the one we shall use)' seems to be 
to approximate the modal data required for the modal ROM (Aj^,, Cj^) 
via the finite element method (FEM) for time-dependent distributed para- 
meter systems (2.9); the data required will be the ROM mode frequencies 
and the corresponding mode shapes evaluated in neighborhoods 

of the actuator and sensor locations when point devices are used. Once 
this FEM approximated modal data is available, the controllability and 
observability of the ROM may be tested, and the FEM based modal controller 
may be synthesized using the approximate modal data in (2.13)-(2.14) as 
though it were exact. Although the FEM approach requires a great deal of 
"number-crunching" to obtain the approximate modal data, this data is 
needed for the open-loop structural analysis and is calculated anyway, 
regardless of control considerations. Also, the large scale computer 
codes necessary for such FEM calculations are available and in general 
use (e.g. NASTRAN and others). 

This makes the FEM approach to structure controller synthesis 
appea"!lng; however, the closed-loop stability analysis is much more diffi 


cult. The fundamental question to be answered is - when does stable 
control of the FEM model produce stable control of the actual structure? 

To begin to answer this question, we must give a more precise descrip- 
tion of the FEM approach. 

A precise description of the FEM approximation of time-dependent 
distributed parameter systems like (2.9) is given, for example, in [2] 
Chapter 7, [3] Chapter 9, or the survey article [4]. We point out 
that, historically, the FEM was developed specifically to numerically 
approximate the behavior of a given distributed parameter system or set 
of partial differential equations; in control terms, it was intended 
to give ? close approximation of the dynamic behavior of open- loop sys- 
tems like (2.9). However, we will use it in a very different way, 
namely to develop reduced-order models of the system (2.9) for synthesis 
of finite dimensional, feedback controllers to be used in closed-loop 
with the system (2.9). In the usual applications, the stability of the 
ROM plays a big role, and hyperbolic systems like (2.9) are notorious for 
their lack, or very low margin, of stability when the FEM is used. In 
our application, the stability of the ROM is less important because the 
feedback controller will be designed to improve its stability; more 
important to our application is the stability of the residual (unmodeled) 
pa-t of the system over which we have no direct control. 

We describe the FEM approach to finite-dimensional controller syn- 
thesis for the distributed parameter system (2.9). Let be an increasing 
sequence of finite dimensional subspaces of the state space H for (2.9). 
Each subspace has dimension N. To make life easier, we assume that 
each is a subspace of D(A) so that its elements satisfy 


OKttHNAL PACK ^ 
OF POOR QUALm 


n 


the boundary conditions for A; however, so called non-conforming elements 
may be used in the more general case. In the FEM, each subspace 
consists of splines (i.e. piecewise-polynomial functions) of fixed degree 
defined over a mesh (usually, of triangles) laid out to approximately 
cover the structure 15 (see [3] Chapter 6). No matter how irregular the 
shape of the boundary of f5 such meshes can be fitted very closely; this 
is one of the principal assets of the FEM. To each mesh, a normal i zed 
mesh parameter h (where 0 < h 1) is assigned so that the mesh is refined 
as h -► 0 and the dimension N of the subspaces increases indefinitely. 


3.1 The Galerkin Approximation 

Let be the orthogonal projection from H into this is called 
the Galerkin projection . The corresponding orthogonal projection onto 
hIJ" is called (i.e. P^^ « I - Pj^). The "rate of convergence" of 
to H is said to be of order q when 


! !Pr i < K h^ (3.1) 

for V in H (or D(A)); this rate is related to the ability of splines in 
to interpolate functions in H. We shall not be concerned with the 
rate of convergence ; consequently we write (3.1) as 


lim ijPj, v|j * 0 for v in H (3.2) 

N-k*> 


and suppress the dependence on n for our discussion. 

Let 'i^^(x), . . . , '<*f^(x) form a basis in (i.e. they are linearly 
independent). These functions are called patch functions or assumed mode 
shapes . An appro,<imation of the solution v(x,t) of (2.9) can be formed 
in by 


Vf^(x,t) 


N 

I v.(t) \{x) 
k=l ^ 


S 


(3.3) 
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i.e. assume separation of time and space variables with all spatial 
variation lumped into the patch functions \^|^(x). The choice of the 
coefficients V|^(t) remains; these are obtained by substitution of (3.3) 
into (2.9): 

3v. 


'N 


3t 


A V|^ + Bf + E 


(3.4) 


where E is the equation error , and the V|^(t)'s are chosen so that 
= 0 


(3.5) 


This is called the Galerkin approximation ; when it is carried out with 
the spline subspaces described above, it produces the FEM approxima - 
tion (3.3) where the coefficients Vj^(t) are given by the entries of the 
solution vector v^^(t) = [vT(t), . . . , Vfj(t)]' for the system of 


ordinary differential equations: 

»N = \ ^ Bn ^ 

where = [(4^^,, ^,^)] , 

» 4 4j^ ) 1 , 

and Bf^f » [(4^^, Bf)]. 

The matrix is symmetric and positive definite because {4|^(x)}|^..| are 
linearly independent: 

^ . N N N N 

z ‘ z * I : z,(4p,4jz^ = ( I lO and 


(3.6) 


t«l k=l 




£*1 


rv 


k=l 


0*1 % 1 = 0 which leads to the Z|^ = 0 

due to the linear independence of the 4|^. Therefore, (3.6) can be solved 
uniquely for_Vj^(t) whenever _Vf^(0) is specified, and hence the FEM approxi- 
mation (3.3) is obtained. It is assumed that Vj,^(0) is given by the 


vector of coefficients of 
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“ Pn ''o 


(3.7) 


expanded In the basis • Note that 


"n * Pn ''n 


(3.8) 


The approximation (3.3) is called a semi discretization of (2.9) because 
time t remains continuous in (3.6). 

It should be noted that to obtain the most analytical benefit from 
the FEM, the approximation (3.3) should be obtained from the “weak" form 
of (2.9); however, we omit discussion of this technicality and refer the 
interested reader to [3]. 

3.2 Convergence of the Approximation 

Let e(x,t) = v(x,t)-V|^,(x,t) , the FEM solution error , and consider 
from (2.9), (3.4), (3.5) and (3.8): 

''n ^'O (3.9a) 


de 


St * ^RN ''n ■*' 


e(0) 




3.9b) 


''n * 


(3.9c) 


where « F^ A P^., B, « CP^, B^^ = Pj^ B and ^ 

Since Pj^ is orthogonal projection onto the finite dimensional subspace 
Hj^, (A|^, Bjj, Cj^) may be identified with appropriate matrices easily re- 
lated to the ones in (3.6) by choosing the orthonormal i zed basis for 

N 

Obtained from addition, is a bounded . linear operator. 

In general, Aj^j^ is not zero, and hence and are not reducing sub- 
spaces for A; furthermore, ''n ~ ® are not orthogonal 

to each other. Consequently, except for th'* special case where the exact 


modal data is known and is chosen to be a modal subspace, the FEM 
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does not partition (2.9) into the previous form (2.13) analyzed in [17]. 
The FEM partition (3.9) must be analyzed separately. 

From (3.9b) and the fact that A generates the Cq semigroup U(t), 
obtain: 


e(t) « U(t) Vq + U(t-T)[Ap^ v^,(t) + Bf^f(t)] dr (3.10) 

Take norms in (3.10) and use (2.11) to obtain: 

||e(t)l| < M|3 e-='(MP„ vjll . v(t)|! 

♦ l|B„ll !|f{T)|| + |!e(T)m Ht) (3.11) 

Now use Gronwall’s Inequality ([25], p. 124) on (3.11) to obtain: 
ne(t)|! e‘""’(||P„ v„!i . ;* »(t)il 

!!f(T)j!] dt) (3.12) 


Where = e - 1 1 Aj^,^ j 1 . 

convergence result: 


This can be used to prove the following 


Theorem 3.1 : The following two conditions guarantee convergence, over any 

finite interval of time, of the FEM approximation (3.3) to the actual 
open-loop solution of (2.9): 


(a) !|A 


RN’ 


< a where a is independent of N 


(b) The approximation is consistent , i.e. lim MA^j^ v|| » 0 for 


N^oo 


all V in D(A) . 

Furthermore, conditions (a) and (b) are equivalent to 

lim MArn'I'O (3.13) 

N-*-* 


Proof: From [18], p. 151 lemma 3.5, we obtain that (a) and (b) imply 

(3.13) because D(A) is dense in H. Conversely, it is easy to 
see that (3.13) inrolies (a) and (b); so 
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these conditions are equivalent. Use (3.12) with replaced by e-MqO 

enc let 0 <. t <_ T. From (3.2) and (3.13) , obtain lim |je(t)l| ■ 0. I 

N** 


This is only a slight modification of the usual open-loop convergence 
results (e.g. [5] Sec. 4), but it will be convenient for our needs. 

A priori estimates of the convergence rate of e(t) as a function of the 
resh parameter h can often be made. 


3.3 An Alternate Form of The FEM Approximation for Flexible Structures 
We want to carry the development of the FEM approximation (3.3) and 
(3.6) a bit further. We can rewrite (3.3) as 


v..(x,t) = [u^(x.t), Wj^(x,t)3 


T 


v.mere 



Z V, (t) (x) 

k»l ^ 

V (t) t 2(x) 

k“ i 


(3.14a) 

(3.14b) 


•j 2 T 

where ] is in Now the entries in the matrices Mj,j, 

A,^, in (3.6) are given by: 


■ ('ft'. *0 *k’>j ^ 

(C..A»k) . Aj tk'! ■ *0 

U CO 

(c.,Bf) « I (o/, b.) f. 

^ i«l ^ 0 ^ 


(3.15a) 

(3.15b) 

(3.15c) 


From '•his data (3.6) may be solved for the coefficients V(^(t) of the FEM 
approximation (3.3). 

An alternate form of the FEM approximation is obtained by choosing 

0 0 
spline subspaces S^ in Hq with 6^(x), .... Sj^(x) as a basis '^or . 
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Instead of making the approximation in the state space formulation (2.9), 
we do it directly on (2.1): 

(3.16) 


N 


(3.17) 


* I u.(t) e.(x) 

where U|^(t) satisfy the system of ordinary differential equations: 

u + b u u^ « f 

where “ [u^ . . . , Uj^]''' , « [(5^,, * C(©;^. Aq 

and * [(6^^, b^. )]• It is easy to see that Mj^ and are positive defi- 
nite matrices because is a self-adjoint, positive differential operator. 

From the extended eigendata problem : 


we obtain the approximate mode frequencies 


“k * 


1/2 


(3.19a) 


and approximate mode shapes : 


s'*’ ■ j", 'kt 'k‘'<> 


wnere 


t^ki 
Vn ■ Ci, . ■ 
we find that 




(3.19b) 

When we form the nonsingular matrix 

(3.20) 


'N 




V A ^ V » A ^ 

''n ''n 


diag [X 


1 


Xn] 


(3.21a) 

(3.21b) 


The modal coordinates q are defined by 


H.H 


''N -N 


( 3 . 22 ) 


and (from (3.17) we obtain 

£m ^ in S.N “ '^N ^ 


(3.23) 
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[ 

[I 


f 


r 


r 


t 


i 



I 


t 


9 


9 


t 




This system of ordinary differential equations is decoupled into N second 
order systems and, hence, is quite easy to solve for . From (3.22), 
we can obtain the coefficients for the alternate approximation (3.16). 

The alternate formulation (3.16) of the FEM for flexible structures 
is the standard one used in structural analysis; it is related to the 
one obtained in (3.3) or, equivalently, (3.14). Appropriate versions 
of (3.2) and (3.9) and a convergence result similar to Theo. 3.1 can be 
produced; fcr example, see [2] Sec. 7.3. For convenience we shall con- 
tinue to use (3.3); however, alternate versions of our later results 
hold for (3.16). 

4.Q Feedback Controllers Associated with the FEW Approximation 


The FEM reduced-order model associated with (2.9) is defined on 
and given by 
, 3v, 


'N 


n 




v„(0) • Vj 


(a.l) 


^ y ' 'n 

where (Aj^, Bj^, Cj^) are defined in (3.9) and are determined by approximate 
modal data. Since is a finite dimensional subspace, (Aj^, Bj^, C^^) may 
be identified with their matrices in an appropriate basis of and 
(4.1) is equivalent to a lumped parameter, state variable system for 
which a well developed feedback control theory exists, e.g. [73] or 
The controllability and observability of (Aj^, Bj^, ere -icsily 
checked via simplified conditions obtained in [16]. These conditions 
are the same even though approximate, rather than exact, modal data are 
used; the conditions apply because of the low level of natural structural 
damping present in these structures (i.e. b is small in (2.1)). 


1 

1 
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The FEM Feedback Controller is based on the ROM (4.1) and defined 
by 




*N ■ \ * ®N' * '^N 




N' 


(4.2) 


where, due to the controllability and observability of the ROM, we can 
adjust the controller gains Gj^ and sc tn: and Aj^ - Kj^ Cj^ 

have any desired stability tnargin. The controller (4.2) has finite dimen- 
sion N (where N * dim H,.); so we can abuse notation by using v^ instead 

3v ^ ^ 

® N 

of . This controller consists of a linear feedback control law and 
full order state estimator (full order in the sense that it is matched 
to the full order ROM). Much lower order controllers than (4.2) may 
be developed, but we leave that issue for another time. 

We define the estimator error ej^ = v,, - Vj^ and, from (3.9c) and 
(4.2), obtain: 


■ h ®N h 


Also, from (3.9a) and (4.2), we have 


(^. 3 ) 


Vn * (Af^ + Gj^) Vfj + Gf^ e^ (4.4) 

If there were no FEM solution error (i.e. e « 0), then (*.3) and (4.4) 
would be designed with a desired stability margin. Consequently, the 
controller (4.2) would stabilize the model (4.1) by design ; however, 
our principal interest is the closed-loop stability of the actual struc- 
ture (2.9) with the controller (4.2). 
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5.0 Closed-Loop Stability; Main Results 


The closed-looD system consists of the actual structure (2.9) with 
the FEM controller (4.2). The closed-loop system state is [v where 
V • Vjij ♦ e and e|^ ■ - Vj^; equivalently, this state may be written 

“ • t''N «K 


Where w satisfies (from (3.9), (4.3) and (4.4)): 


^ • A 

st C 


where A^ ■ 


^1 1 ^1 2 

L ^21 ^22 


and 


'll 


(5.1) 



A * 

1 

o 

0 A|, - Cjj 

‘ • 

. 



^21 * h ^22 * ^ 


Sines A.j.| represents i finite dimensional system, it is bounded and gene- 
rates a Cq semigroup Uj^(t) with the following growth property: 

-Cf.t 

. t > 0 (5.2) 

where i 1 and > 0; c^, is the designed stability margin obtained 
in Sec. 4.0 by the choice of controller gains Gj^, and Kj^. We have the 
following closed-loop stability result : 


Theorem 5.1 : If the FEM model (A^^, Bj^, Cj^) is controllable and observable 

for evc'*y N and if conditions (a) and (b) of Theo. 3.1 are satisfied, then 
controller gains Gj^ and can be chosen so that A^ generates a Cq semi- 
group UQ(t) such that 
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w(t) ■ U^(t) u(0) 


(5.3) 


and 




l|Ut(t)|l < « . t > 0 

«here 1 1 (-| I i (.1 |v„| |^ * I I«nI 1 ^ ♦ I |t| . 


(5.4) 


Me - % Mo (1 ♦ 


£0 • E - Mj S„ 




2 1/2 


S ■ ^ 


"n ■ " 


) 


« • MArnM ^2MBr! 


'N ■ ''"RN" ' ‘■N“RN ^)* ^ *'** (2.11). 

The closed-loop system is stable, if 6j^ is bounded for all N and, for N 
sufficiently large, 

Mj. Bfj < c (5.5) 


Proof: Choose Gj^j and Kj^ so that tfj > c since (A^^,, Bj^, C|^) controllable 

and observable. Let + :.A where 


All ^12 



1 

O 

O 

— — J 


and 

AA * 


CM 

CM 

O 



^21 ° J 


Note that Aq generates a Cq semigroup L'g{t) such that 
!|Uo(t) 1 1 1 Mg e ° . t > 0 

where Cg = m-’n (cj^, c) ■ c; this follows from arguments similar 
to those used :o prove Theo. 5.1 in [17]. Since !1^A|[ • 
llA^iM < 6j^, AA is a bounded perturbation of the semigroup gene 
rator Ag; hence, by [25], Theo. 10.9, ?. 210, Ag generates a Cg 
semigroup Ug(t) which satisfies (5. A). The conditions (a) and 
(b) of Theo. 3.1 imply that 
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11m - 0 

and this, together with (3. 2) and the boundedness of G|^, In^lles 
11m Bm * 0. If N can be chosen large enough to satisfy (5.5). then 
> 0 and the closed-loop system (5.3) Is exponentially stable. 
Note that may Increase as N » so It Is not always certain 
that (5.5) can be satisfied even though 6^^ - 0. * 

Theorem S.2 : In Theo. 5.1, if there exist positive constants a, and 

0-2 such that 

a,-l 

(5 

then there is a finite integer N* such that, for all N ^ N*, (5.3) is 
exponentially stable. 


Proof: This is a corollary to Thee. 5.1. From (6.6), we outain 

K. < a, “ •* 0 as N - Therefore, choose N* such 

w N “ I U 

that a.] eJJI < c and, for all N ^ N*, (5.5) is satisfied; 
the result follows from Theo. 5.1. • 

The above re.»;ults provide conditions under which ti -,2 clcsec-loop system 
(5.3) is stable and the FEM solution error e(t' satisfies 

-Cft 

!ie{t)M 1 !!-3(t;M l 2 Mj- e ^ \\y^\\ , t > 0 (5.7) 

because 


:..( 0)''2 « 2 " P , v .'|2 

Pi W 


' ' V ! 


1 2Mv,' 


.2 


This is an a priori convergence estimate for the closed-loop system. 




i 

22 ' 


Since e « v - - (v - v) + (Pj^ v - and v - v «_ Pj^ v -► 0 

as N « by (3.2), we may be interested in the behavior of the term 
Pj^ V - as N -*• «. However, from the properties of the Galerkin pro- 


jection P|^, 


r?^ e 


"“r * ■ "n* • '’n " - *K 


^ e . Pp P;, V + Pr(P„ V - v„) . P„ V 


(5.8a) 

(5.8b) 


i .e. , Pj^jV - v^i * Pj^ e and v - P^^ v * Pj^ e form the orthogonal components 
of e in and , resoectively. From (3.9b) and (4.2), obtain: 


3 Pn e 
3t 


= An Pn e + An,^ Vf^ 


(5.9a) 


3v 


“ Anil Pm 6 ■*■ An Vn + (Anij ■*■ Bn G.,) V^ + Bp Gm S. 


(5.9b) 


3t ^ "RN "N ' ' '"R 'R ■ '"RN ' ‘'R ir ''N “R '*N "N 
where Vj^ = P^ Vj^, An^ = Pj^ A Pj^, etc. Also, (4.3)becomes 

- S " S S Pn ® " S S ''r 

Let (j * tvN.fif^iPj^ e , be the closed-loop system state, and, from 
(4.4), (5.9) and (5.10), obtain 


3lO n 

t "> C “ 


where * 


ni ^12 


L^21 ^22 


and 


^11 


71 


An + Bn Gn Bn Gn 

“1 


^ - 

0 

0 An - <N Cn 

'■nI 

1 , Qi2 * ^ 

S ^R 

o 

o 

% . 

1 

%R 

m m 

^RN ®R S ®R ^RN^’ 

and 

Q 22 ' Aj^ 



(5.11) 


c 
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The following alternate closed-loop stability result may prove useful in 
those cases where conditions (a) and (b) of Theo. 5.1 are not satisfied 
by the FEM approximation: 

Theorem 5.3 : Assume, for every N, 

(1) the ROM (Ajj, is controllable and observable 

(2) A^j is (exponentially), stable with margin 

(3) Ap generates an exponentially stable Cq semigroup Up(t) with 

|!Uj(t)|| e . t > 0 

«) Hi" I|a„rII ■ 0 
00 

Then controller gains Gj^ and can be chosen so that in (5.11) gene- 
rates a Cq semigroup UQ(t) with 

-^rt 

MUQ(t)|| iMq e ^ . t > 0 (5.12) 

where Mq = max (M.|i, M^) , 

N'^r! 1 ^ ' 

and Q.J 1 generates the Cq semigroup U.|i(t) with 
!IUn(t)i! < e , t > 0. 

If K|^ is bounded for all N and, for N sufficiently large. 

Me !•„ < ^0 <5-'3) 

the closed-loop system (5.11) is exponentially stable. 

The proof of this result is sufficiently similar to that of Theo. 5.1 
that we shall omit it. In some cases, the conditions (l)-(4) of Theo. 
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5.3 may be easier to verify than those of Theo. 5.1. Again, the stability 
of (5.12) leads to an a priori convergence bound on the FEM solution 
error: 

|le(t)ll < 2Mj. e ^ IIvqII . t > 0 (5.14) 

6.0 Conclusions 

The use of the finite element method (FEM) for numerical simulation 
of complex distributed parameter systems is well known. In this paper 
the use of the FEM to produce reduced-order models and finite dimensional 
feedback controllers for distributed parameter systems ha-s been investi- 
gated with emphasis on application to control of large aerospace struc- 
tures. Our principal concern has been to determine the stability of 
tne controller in closed-loop with the actual distributed parameter 
system even though the controller is synthesized using a FEM approxi- 
mation. Our main results. Theorems 5. 1-5. 3, give conditions under which 
stable closed-loop operation can be expected. Such conditions are 
especially desireable for large aerospace structure control because the 
FEM model is the tool most widely accepted for structural analysis and 
large scale computer codes already exist to produce such models for very 
complex structures. 
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ABSTRACT 


Implementable feedback control of distributed parameter systems must 
often be based on reduced-order models due to the infinite dimensional 
nature of the actual open-loop system. In this paper, we analyze the be- 
havior of controllers designed via reduced-order models obtained with singu- 
lar perturbations techniques. When such controllers are used in the actual 
distributed parameter system, the closed-loop stibility is in question. The 
results presented here provide bounds on the smallness of the singular per- 
turbations parameter to ensure stable operation; such a priori bounds may be 
used to evaluate potential reduced-order controllers for distributed para- 
meter systems. 


1. INTRODUCTION 


Distributed Parameter Systems (DPS) are described by partial differen- 
tial equations. Examples of DPS include heat diffusion and chemical pro- 
cesses. wave propagation, and mechanically flexible structures, such as 
large flexible spacecraft and satellites, high-speed aircraft and surface 
vehicles, and large civil engineering structures, e.g., bridges and tall 
buildings. 

The state spaces for these systems have infinite dimension. Thus, it 
is impractical or impossible to implement feedback controllers based on 
confilete models of these systems; hence, reduced-order models must be used 
for the controller design. A great variety of reduced-order modeling tech- 
niques exist for general systems, e.g., the bibliography of [1], and new 
techniques are presently being developed in the specialized areas lik» large 
aerospace structures, e.g., [2], [3]. Of particular interest are the model 
reduction techniques based on asymptotic methods, such as multiple time 
scales and singular perturbations, e.g., [4]-[7] and the excellent survey [8]. 
These papers address the large scale, but finite dimensional or lumped para- 
meter, systems; to our knowledge, very little has been done with asymptotic 
methods for DPS with the exceptions: [9] Chapt. 5, [10], and the survey [11]. 

Our interest here is less in the area of derivation of reduced-order 
models for DPS (although that is an in^ortant topic in its own right); 
rather, we are concerned with the successful operation of the controller, 
designed on a reduced-order model, in closed-loop with the full system. Other 
work on this basic topic may be found, for example, in [2], [12]-[15], and, 
for singular perturbations methods, [16]-[24]; none of these singular pertur- 
bations papers deals with DPS, with the exception of [24]. This paper 
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extends the preliminary results on singular perturbations and stability for 
DPS obtained In [24] and obtains an upper bound on the singular parameter 
for stable, closed-loop operation. Our results for DPS are In the same spirit 
as those of [21]-[23] for lunged parameter systems; the method of proof 
Is different due to the Infinite dimensional nature of the DPS problem. 

We consider linear , time-invariant distributed parameter systems with 
the following form: 

. * “ 0 - 1 ) 
y « Cv Df 

where the state v Is In a Hilbert space H with Inner product {.,.) and 
associated norm H-H and the control vector f and observation vector y 
have dimensions M and P, respectively, which denote the number of (Inde- 
pendent) actuators and sensors. When the system (1.1) Is a distributed para- 
meter system (DPS), the dim H ■ <» and the operator A Is (usually)an un- 
bounded, differential operator with domain D(A), containing all states which 
satisfy the boundary conditions of the problem Srde ns e in H, while the in- 
put and output operators B and C have finite rank M and P, respectively, and 
D is a PxM matrix. 

The operator A generates a semigroup U(t) on H. This semigroup 
U(t) Is usually exponentially stable, I.e., it has the growth property: 

-6nt 

l|U(t)|| ^ , t > 0 (1.2) 

where > 1 and 6 q > 0. In physical systems of "hyperbol1c"-typejenergy dissi- 
pative mechanisms make (1.2) true, even for DPS; however, 6 q may be quite 
small . as it often is In the case of large aerospace structures [2]. 

In the following sections, we develop the general idea of reduced-order 
models for the systems of type (1.1) and obtain results on the operation 
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of controllers, based on these reduced>order models, In closed-loop with 
the actual system (1.1). 

2. REDUCED-ORDER MODELING <\ND CONTROLLER DESIGN 

Let and be subspaces of the total state space H with dim > 

N < • and H ■ Define the projection operators Pj^ and Pj^ (not 

necessarily orthogonal) and let v and Vj^ • Pj^ v. This decomposes 


V into V 


+ Vp and the system 1 

[1.1) into 



Vn- 

■ 'Vi ' 

'n ^ ®N 

f 

* Vj^(O) - 

’’n *0 

(2.1) 



''n ^R ''r ®R 

f. 

. Vp(0) - 

'’r *0 

(2.2) 

y • 

■ ' 

'n ■*“ ‘'R ''r ^ ^ 




(2.3) 

Note that all 

parameters, Aj^, 


etc. , with 

the exception of A^^, are 


bounded operators because they Involve projection onto the finite dimen- 
sional subspace H|^. We abuse notation slightly by writing v^ for 3Vg/3t. 

Henceforth, we assume that A|^ generates a semigroup Uj^(t) on which 
satisfies the growth condition: 

lIUfj(t)I! < Kp e , t>0 (2.4) 

where Kj^ > 1 and > 0. In the special case of reducing subsoaces (i.e., 
both and are A-invariant) , ■ 6 q and Kp ■ llPpH Kq where 6g, Kg 

are as given in (1.2). In general, K|^ and 6^^ depend on the choice of 

subspaces. The growth condition (2.4) will hold for "parabolic", as well 
as "hyperbolic" systems even though the more stringent (1.2) does not hold. 
The reduced -order model (ROM) for the system is (2.1) and (2.3) with 
and Af^p assumed to be zero: 

I ’n ■ 'n * B), f ■ Pn "o 


(2.5) 
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Thus, the ROM depends on what choice of subspace Is made and what type 
of projection is used (o»* alternatively, what Is). The subspace 
Is called the ROM subsoace and the subspace H|^, the residuals subsoace . 

The terms and A|^ V|^ are called model error and B|^ f and C|^ V|^ are 

called control and observation spillover , respectively. 

3. REDUCED-ORDER MODELING AND CONTROLLER DESIGN VIA SINGULAR PERTURBATIONS 

For certain choices of subspaces In Sec. 2, It Is possible to 

produce ROM's such that (2.1)-(2.3) become: 

'■'n * ^1 ''n ^1R ''r * ®N ^ 

® ''r “ ^RN * ^R ''r ®R ^ 

y ■ Cj^ Vj^ ♦ Cj^ Vj^ + D f (3.3) 

where the singular perturbation parameter e > 0 represents some small para- 
meter dependence, such as electrical networks with parasitics or dynamical 
systems with small masses and time constants (see [25]}; alternatively, e 
may represent a ratio of time-scales In the system (see [7]). This approach 
Is especially valuable when the ROM Is based on vibration modes for a mech- 
anically flexible structure and the frequency separation of slow and fs:;t 
modes can be used to provide the deconposition of (3.1)-(3.3). A survey of 
singular perturbation techniques for model reduction Is presented In [8]. 

It should be noted that one of the most difficult tasks may be to formulate 
a particular DPS pn^lem Into the singular perturbations format (3.1)-(3.3). 

The parameters such as A^^ may Involve regular perturbations In c, as 
well, I.e., 

A,^ • aJ ♦ e aJ (3.4) 

for all e and ^ * 0- Also, It Is possible to have multi 


v4)ere ^ Is bounded 


I 
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parameter singular perturbations where several c parameters are present; see 


[26]>[27]. To slR^lIfy the presentation, we shall not consider either 


of the above situations here. 


The singular perturbations R(^ Is oL*«1ned when e ■ 0 and Is given by 


^ ^ ''O 


(3.5) 


I y “ ^ ^ ^ 

where ^ A,^ ^ - Aj^j^ A^ ^ ^R ^ ^RN’ 

and U • 0 - Cj^ Aj^’^ Bp. This Is somewhat different from the usual RW (2.5) 
due : the parameter correction terms Involving Ap’^. In the special case of 
reducing subspaces Hp, Hp (as In [24]), the model error terms Aj^p and App are 
zero, and ^ * Ap, Bp ■ Bp, and Tp ■ Cp; however, 0 D. Therefore, even In 


this special case the (3.5) Is not quite the same as the usual ROM (2.5). 


Note that the ROM (3.5) Is finite dimensional (dim Vp ■ N) and the para> 


meters ()Tp, fp, 7p, U) may be 1<tent1f1ed with their corresponding matrices In 


an appropriate basis for Hp. Henceforth, we assume the ROM (3.5) Is cen- 
trollable and observable ; for any ROM, this Is easy to verify with the the 
usual ranL tests for finite dimensional systems ([28] Chapt. 11). 


Also, note that, although Ap Is usually not bounded for a DPS, Ap' 
1$ bounded due to (2.4). In fact, from [28], Theo. 8.9, 


(3.6) 


The reduced-order controller based on the R(W (3.5) Is given by 


f 

z ■ )Tp z + ffp f l^(y-y) , z(0) • 0 

y ■ ^ z ♦ U f 


(3.7) 


A 
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where dim z * N. This finite dimensional controller can be implemented for 
the DPS, and the gains G|^, can be designed for closed-loop stability when 
c > 0; however, successful closed-loop operation is in question when e > 0. 


4. MAIN RESULTS: CLOSED-LOOP STABILITY 


The fundamental question is whether a stable closed-loop reduced-order 
system (e = 0) will remain stable when the same controller (3.9) is used 
witp the actual system (3.1)-(3.3) when e > 0 (although small). The answer, 
as given belcw, is that it will when e is small enough ; bounds on the small- 
ness of e are obtained. 

The closed-loop behavior of the total system (3.1)-(3.3) with the 
controller (3.7) is governed by 

w. 


f 


■"11 

Hi2 


S 



1^2 - 


.^21 

“22- 




l'^2 J 


where w, = [vl el]^ , e. = z - v^,, w, = Vp, and 


'll 


'N 'N- 






-1 




^RN 


(4.1) 


'12 


'NR 


L h ^R " ^NR 


'21 


-^RN h h 


and H 22 ® A|^. Note that ^^21 bounded operators, and H 22 

is bounded by (3.6); therefore let 


-1 
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1 1^22 1 1 £. ^4 = * 1 1^22 ^ 21 11—^5 

When e = 0, (3.1) reduces to 
»1 * »1 

which represents the reduced-order design consisting of the ROM (3.5) 
the reduced-order controller (3.7) in closed- loop. The reduced-order 
loop operator is given by: 


(4.2) 


(4.3) 

and 

closed- 


Hi * - H ^2 ^22 ^ ^21 


^ h h h S 


0 J 

which can be made stable by the choice of gains 6|^, K^, since the ROM 
is assumed controllable and observable in Sec. 3. Therefore, let the 
tion matrix U^(t) associated with satisfy: 

-6,t 

||U^(t)|l < e ' . t > 0 


Furthermore, is bounded: 


where 


IlHill iMg 

Mg = Ml + M2 Mg < K^/6 t 

When e > 0, we can write (4.1) as 


(3.5) 

transi- 

(4.4) 

(4.5) 


w = H w = 


"ll 

“l2 


"22 


w 


(4.6) 


T .T-.T 


where w = [w^ w^] • The primary question to be answered in this section is 
whether the semigroup U (t) generated by H in (4.6) remains stable for small 
positive e, i.e.. 


Uc(t)l| 1 e 




t > 0 


(4.7) 
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The answer is that it does remain stable and the following theorem — our 
main result gives a bound on the smallness required of c: 

THEOREM 1 ; Let the controller gains Kj^ be chosen so that in (4.3) 
has stability margin 6-| as in (4.4). There exists §n c^ > 0 such 
that, for all 0 < e < Cq, the closed-loop system (4.6) — consisting 
of the full-order DPS (1.1) and the reduced-order controller (3.7) 
based on the ROM (3.5) — is (exponentially) stable and the controller 
state z converges (exponentially) to the reduced state v^; this means 
(4.7) is satisfied. An upper bound for Cq is given by: 

. ^1 

^0 (1 + Mg) K,^ Mg (M2 + Mg) ‘ ^ 

where y is defined later in Lenrnia 1, and K , 6 , in (4.7), are given 
by: 

6^. = min (6^ , ^Cq) (4.9) 

2 2 

Kc = Ki Kj^ (1 + a + a^) (3 + 3 Mg + Mg"^) (4.10) 

where 

6^ = 6^ - £q M 2 M^ (1 + Mg) (M 2 + Mg) 

^2 - ^R " ^0 “^R ^2 * ^5^ 

Gt = M2/IAI , A = 6-j 62/^0 , 

and the constants are obtained from (2.4), (4.2), (4.4), and (4.5). 

The proof of Theo. 1 requires the following two lemmae: 

LEMMA 1: There exists e.| > 0 such that, for all 0 < e £ the nonlinear 

mapping h(L) = H 22 '^ H 21 + e H 22 *^ L (H-ji - H-j 2 L) defined on 
n * {L I ||L - H 22 ~^ H 21 II < 1} has a unique fixed-point L* = L*(e) 
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in n (hence, ||L*|| < 1 + Mg). Furthermore, Z 2 = Wg + L* trans- 
forms (4.1) into 


*1 


cz. 


- e 


H 


12 


H22 + e L* H^2 


w 


1 


(4.11) 


-1 


where = H -|2 H 22 L* (H^^ - H ^2 •-*)• upper bound for e-j is 
[M^ (1 + Mg) ( M 2 + Mg)]"^ and we write * yCM^ (1 + Mg) (M 2 + Mg)]"^ 
where 0 < y < !• 

LEW1A 2 : There exists Cq such that, for all 0 < e < Cg , 


1/2 


-6.t 


(4.12) 


[|!wi(t)ir + l)z2(t)ir] iK^e 

[!1w^(0)1!2 + ||z2(0)!12]^^^ 

2 

where = K-| Kj^ (1 + a + a ) and a, 6 ^ the same as in Theo. 1. 

An upper bound forCgis given by (4.8). 

The oroof of Lemma 1 involves a contraction mapping argument ([29], 
Theo. II. 1.1, p. 24); however, the results of [30] might be used to obtain 
a different (possibly, larger) bound for e-j. It is well known ([29], p. 24) 
that 

L* » lim L|^ 

where Lj^ is obtained via the algorithm: 

^k+1 " * k = 0, 1, 2, . . . , (4.13) 

and Lg any menter of ft. This algorithm may be implemented to calculate the 
desired L*. The following gives an indication of the convergence rate: 


1 14 " '■*1 1 - I l*-o ’ ^*1 1 


(4.14) 


where y is defined in Lemma 1. 


OHIGINAL 
OF POOR QUAliTY 


The proof of Lernna 2 , and hence Theo. 1, makes use of Lemma 1 and 
results from [28]. The proofs of Lenma 1 and 2 and Theo. 1 appear in 
Appendices I -I II. 

5. CONCLUSIONS 

Our main result (Theo. 1) provides an upper bound (4.8) on the small- 
ness of the singular perturbations parameter (e) which ensures stable closed- 
loop operation of a finite dimensional controller with the full -order distri- 
buted parameter system. Since the term H ^2 ^nd H 21 In {4*6) appear because 
of spillover and model error, the bounds (and hence, the stability) improve 
when these terms can be reduced (i.e., when Mg and Mg can be made smaller). 
Lenina 11^ taken together form an infinite dimensional version of the 
Klimushchev-Krasovskii result [8]. 

When the distributed parameter system can be put into the singular 
perturbations format, (_3.1)-(3.3) , the stability condition presented here 
can be checked with only a limited knowledge of the unmodeled residuals pre- 
sent in the full -order system. This makes it possible to synthesize low- 
order controllers for distributed parameter systems via general reduced- 
order modeling techniques and analyze their operation in closed-loop with 
the actual system. Such a result appears to be particularly applicable to 
distributed parameter systems with multiple time-scales or high-low fre- 
quency separation, as is often the case in large aerospace structures [2], 
[31]. The general modeling issue, i.e. obtaining a singular perturbations 
format, for large-scale or distributed parameter systems is quite complex; 
see [32] for further discussion. 
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Appendix I: Proof of Lemma 1 

We want to apply the contraction mapping theorem to h to obtain L*. 
Consider any L in n and note that 

lll-ll < ML - Hj,!! + IlHjj’’ Hjill 1 1 + Mj (1.1) 

Also, ||h(L) - Hjj'’ Hj,!! - ellHjj’’ L (H,, - L)|| 

< e M, ||L|| (M, + Mj ||L||) 

£0 «4 (1 + Mj) (Mj + Mj) (1.2) 

because Mg » M,| + Mg Mg in (4.5). Therefore h(L) is in n when 0 < e £ 
where 

04 ' [H 4 (1 ♦ Mg) (M^ + Mj)]"’ (1.3) 

because ||h(L) - H 2 .||| <1. 

Consider L and L' in 12 and, from (I.l), 

||h(L) - h(L')|l • c IlHjj-’ [L(H„ - L) - L'(H,, - L')]|| 

5. ^ M^ (||L(H ^1 - H ,|2 L) - L(H,j,| - H^jg L')[| + 

- H^g L*) - - H^g L*)|I) 

< e M^ (1|L H^gll + I|H^^ - H^g L'||) 1|L - L' 1 1 

< e M^ (2(1 + Mg) Mg + M^) ML - L'M < ||L - L'|! (1.4) 


when 0 < e < Eg where 


-1 


[M^(M^ + 2(1 + Mg)Mg] 


(1.5) 




I 
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But, from (1.3), 

- [(1 + Hj) M, (M, t Mj (1 + Mj))]'’ 

< [(1 + Hj) (H, ♦ 2 Mj (1 + Mj))]'’ 

■ ‘b'O 

Therefore, < (1 Mg) * £g; hence, we choose and, from (1.4), 

h is a contraction on J2. By the contraction mapping theorem ([29] Theo. 

II. 1.1, p. 24), h has a unique fixed point L* in for any 0 < e < , i.e. 

L* » h(L*) » H22"’ H2^ + e H22' ' L*(H^^ - H,2 L*) (1.6) 

and, from (I.l), 

1!L*11<1+M5 (1.7) 

Note that - H^2 L* * - e and, by substitution of Z2 into 

(4.1), obtain + H^2 ^^2 - L*w.|) » (H^ - e W^) + H^2 ^2 

9 0 

e Z2 = e(w2 + L* w-|) 

» H21 + H22 W2 + e L* ((H^ - e W^) w^ + H^2 ^ 2 ^ 

= H22 [h(L*) - L*] + (H22 + e L* H^2^ h 

» (H22 + e L* H^ 2 ^ h 

This is the result desired. # 
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Appendix II; Proof of Lemma 2 

Assume 0 < e < Cq < e.| and, from Leimia 1, obtain 
- (H^ - e W^) + H^2 ^2 


eZ2 • (H22 + e L* H^g) Zg 


Choose T ■ t/e and ZgCt) ■ ZgCt). Therefore, from (II. 2) 

A 

^ ■ ^2 3? ■ ' 4 • <”22 * “■* "l2> *2 • 


(n.i) 


(II. 2) 


Hence, from H 22 * Aj^ and (2.4) and [27], p. 215, 

A 

I|22<^)II <i^r ‘ ^ l|J2<»)ll (”•3) 

where ^2 ■ - Cq (1 + Mg) M 2 Kj^ because e ||L*H^ 2 M ^5^ ^2’ 


Choose 


-1 


so that &2 ^ ( 11 . 3 ) and ]|z2(t)ll » 11 z2(t)| 1, obtain 

4 ‘ 

l|Z2<‘>ll IUrOIII 


(II.4) 


<K„e ^0 ‘ ||Z2(0)|| 


because e < c 


0 * 


,-l 


(II. 5) 


Choose Eq < 6^ [(1 + Mg) M2 M^ (M2 + Mg) K,] ' (II. 6) 

where 5.j and K.j are given in (4.4); assume £g satisfies (II.4), as well. Let 
5.| » 5.J - Eq K.J Mg M^ (1 + Mg) (M2 + Mg) which is positive by (II. 6). Since 
I |W-| 1 1 < ^2 M^ (1 + Mg) (M2 + Mg), we have, from (4.4), (II. 1), (II. 5), and 
[28], p. 215, 
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b I 

^ f 


n 


|w,(t)ll <K^ [|h(0)|l *n^fl I |ds3 


-6,t 


< e ^ Cl|w^(0)|l + Mg K,^ /J e^* ds ||Z2(0)||] 


•dit 


< K, e'°’‘ [ 1 |»,( 0 )|| K„ ] 


(11.7) 


A 

where A ■ 6, - — - . 


’ ‘o .“St 

Eq 

If d > 0, then (II.7) implies Ilw,(t)|| < e ^ (ilwi(0)|i+ 


«2 


^ ||Z2(0)|1). If A < 0, then (II. 7) Implies |Iw^(t)l| < e 


-«lt 


MgKR 


Wl(0)|| + -jjp l|z 2 ( 0 ) 11 ). In either case, if 6^ ■ min { 6 ^, ~), then 


-6.t 


■S! 


Il*<i(t)ll < ' '''(lIxidDII +-jsf lUidDII ) 

and, from (1 1. 5) and (1 1. 8) and K.| and Kr ^ 1, 

-V 2 


(II. 8) 


||w^(t)|l2 + i ® (IU 2 {o)M^ + (llw^(o)ll 

Ma a -6 t 2 A 

+ f|^ ||Z2(0)|I)^) < (K^ Kr e ^) (l+a + a-^) 


{llwi( 0 )ll 2 + | 1 Z 2 ( 0 ) 11 ^) 

Ma 

where o ■ and the easily verified inequality: 

+ (a + b)^ < (1 + o + a^) (a^ + b^) 
is used with a ■ ||w^(0)|| and b ■ |Iz2(0)|l. Therefore, 


1/2 


- 6 ^t 


w^(t)ir + 1122^^)11 ) i*^c ® + I1Z2^0)II ) 


1/2 


(II.9) 


Wt imist choose Cq to satisfy (II.4), (11.6) and Cq < ■ Y * 

{1 ♦ Mg) (M 2 + where 0 < y < 1. Note that, frow (4.2), M 

Kj^/6p. Therefore choose Cq to satisfy: 


and this will meet a11 the requirements. However, we can refine (11.10) 
further because 0 < y < T Md M 2 CM 2 + Mg]"^ < 1; consequently, (11.10) 
becomes: 




which Is the desired result. I 
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Appendix III: Proof of Theo. 1 

Assum 0 < e £ Cq« where Cq given by (4.8). 


■ w 

♦ L* w^ , 


V ^ mi 



' 

- Q 

’ 

'I 

CM 

AM 

1 


CM 


Consider thet, from 


(III.l) 


where 




’i o' 



■ I o' 


Q • 


and 


1 


.L* I, 



-•«-* I. 


From Lemma 2, (II. 9) and (III.l), obtain 

l|w(t)l| < I IQ*’ 1 1 1 151 1 IWOlll 


(III.2) 


But IIQwll . 4lw,||2 ♦ l|uj ♦ L* w,||^ 

? 1/2 2 

> 1(1 MiL*ii HiL*ir) iNir (HI. 3) 

by use of the easily verified inequality: a^ + (b ♦ oa)^ < (1 + a + a^) • 

(a^ ♦ b^) with a ■ 1 Iw^ } | , b ■ | jw^l | , and a » | |L*| | . Therefore [ jQ] | < 

* (1 + ||L*|I + IlL^ll^)^^^ and. similarly. 1|Q‘^|| < (1 + ||L*M ♦ \ 

Note that, from Lemma 1 (1.7), ||L*|I < 1 ♦ Mg. When these inequalities are 
substituted into (III. 2 ;, it is clear that the desired result is obtained. # 
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A3STS/.CT 

Using s siagulsr psrcurbscioQ fosaulsclon of che liaser cias-lsvsrisnc distributed 
psrsacter systea, we develop e aethod to design finite-diaeasionel feedback cospes- 
setert cf any fixsd order which will stabilize the infinite>-dlBeQsio&al distributed 
paraaeter syatea. The synthesis conditions are given entirely in cervs of a finice- 
diaensimal reduced-order aodel; che scabiUcy result depends on an infinice-diaen- 
sional version of che Kliaushchev-Rrasovskii Icasna developed in tlOl. 


1.0 imoDuenoN 

Many engineering systns exhibit a distributed paraaeter nature and auac br des- 
cribee by partial differencial equations. Exeaples of such distributed paraaeter 
syscecs (DPS) include heat diffusion 4 chcaical processes, wave propagation, 4 
aechanically flexible structures. Various aspects of the cracrol of DPS have been 
ccnsiderec in, for exaaple, [I] - 15]; our experience in DPS has been shaped by 
applications in large aerospace structures [6]. 

7r.e state spaces for DPS have infinite diaensien; sc, at best, reduced-order models 
3 :st be used in controller synthesis. However, che closed-loop stability of the 
itiinice dlaensiooal I^S with a finite diacnsicnal feedback concroller becomes a 
fundacental issue. The synthesis of finite dimensional controllers for DPS and the 
analysis cf their closed-loop stability by singular and regular perturbation techni- 
ques have been our main areas of emphasis [7]; this theory has been developed with 
flexible structures and ocher highly oscillatory DPS applications in aind. 

Sven in large-scale, lusped paraaeter syscesm, such as electric power distribution 
networks, it is necessary to perfora aodel rciuccion ai^ reduced-order concroller 
synthesis and to analyze closed-loop stability. The use of asyaptocic methods, 
especially singular perturbatiMis, has been vary successful in this regard (e.g. 
[S]). Vs have extended certain of these singular perturbations methods for DPS 
to provide estimates of stability in an infinite dimensional setting [9) - [10] and 
applied th» to Mchanicall'* flexible structures 111). In this paper we will use 
these singular perturbation results to analyze che general finite dimensional com- 
pensator for llnaar DPS. 
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Wc consider linesr, tlsc-inverient DPS with the form: 
• Av + Bf, v(0) • Vq 

y - Cv 


( 1 . 1 ) 


where the state v(t) is in a Hilbert space H (called the state space) for all 
taO; the state space has inner product (*,•) and associated norm ||*||. The 
control vector f and the observation vector y have dimensions M & F, respectively, 
which denote the number of (independent) actuators and sensors in use; thus, rank 
B«M and rank C«P. the operator A is an unbounded differential operator with 
domain D(A) dense in E and A generates a C. semigroup U(t) on H. Hany practical 
DPS can be stated in this form [7]. 


2.0 REDDCED-ORDER MODELING OF DPS: A SINGULAR PERrjRBATIONS FORMULATION 


Since the state space B of the DPS in (1.1) is infinite dimensional, we must obtain 
" a reduced-order model (RDM) upon which to base the finite dimensional controller 
design. In general, this is done by selecting a finite dimensional subspace. 

(with dim. Hj; • N < «) contained in D(A) . Thia subspece Hn is the ROM subspace ; 
its complement Hp is the Residual subspace, and together they decompose the state 
space: 

t ( 2 . 1 ) 

We define the (not necessarily orthogonal) projection operators P^j and Pr and let 
v^ ■ and vj^ ■ Pj^v. These decompose the state v ■ vjj + and the DPS (1.1). 


In this paper, we will make two basic assumptions on model reduction processes 

(1) the subspaces Hj;, Hr are reducing subsnaces (i.e. they are both A- 
invariant) ; 

(2) the subspaces Hr, Hr^ can be chosen to give a singular perturbations formula- 
tion with small parameter eiO and an eyponentiallv-stable residual subsystem. 


These assumptions yield the decomposed DPS (1.1) in the following form: 


— ■ S "s V 

(2.2a) 


(2.2b) 

^ ■ S ’r 

(2.2c) 

where A»^ ■ PrAPr, Br ■ PrB, Cr ■ CPr, etc. and Ar generates the Cq semigroup 
UR(t) • PpU(t)Pp with exponential growth condition: 

II Uj^(t)|| < e“^R^ t i 0 

(2.3) 

where Kj^ i 1 and > 0. All operators except Ar are 

Cp VP are called, respectively, control & observation 

bounded. The terms BRf & 
spillover; they represent 


the interconnections through which the controller can affect the residual sub- 
system. 

The ROM of the DPS, in this formation, is obtained by setting £-0: 


t 




(2.4) 


[ ^ “N *N “N* 

where Di^ > - ^R'^R ^R 4 ^ is a bounded operator due to (2.3). The ROM is a 

finite disensional system and the parameters (A>;,B^,Ci},%) identified with 

matrices by choosing a basis in %; this will be done wherever necessary. Ue 
shall assume , henceforth, that the ROM is controllable and observable (easily 
verifiable conditions exist for finite dimensional systems) . The ROM is completely 
determined by the choice of the subspaces H}] & Eg. 


3.0 FliJITE DIMENSIONAL CCSIPENSATOES FOR DPS 


The form of the finite dimensional feedback compensators used here will be the 
following: 


f - L^^y + 


(3.1) 


where the conpensator state z has dim. z ■ a < N < ® and Ln, 112,121*^22 
matrices of appropriate sizes. We say the compensator is output feedback when 
ct-0. The order a of the compensator is assumed to be fixed at some acceptable 
value which reflects the available capacity of the on-line computer being used to 
implement (3.1). lo [10], the conpensator order was a « N. 

The compensator design is synthesized as though the ROM (2.4) were the full DPS 

(1.1), i.e. as though £*0. Let 


where 




Fs, “ 


(3.2) 


Hi “ ■ ^11 ^11 

H2 “ " ^11 V ^12 

Hi ■ Hi^H ^ °N ^11^ 

H2 “ H2 Hi °N ^12; • 


Let the composite state qjj ■ 
(for £-0) : ^ * 


in H» X R ; from (2.4) & (3.1), this satisfies 


ar-^% (2.3) 

The following theorem gives the conditions under which a stable design can be 
synthesized: 


Theorem 3.1 : If the ROM (A>j,3i5,Cjj,Dj5) is controllable and observable and 

M+P+aiN+1 (3.4) 

then L may be chosen so that F given by (3.2) has any desired pole locations in 


I 


the complex plane. 

The proof follows easily from results in [12] or [13] and depends on_the finite 
dimensionality of the closed-loop system when e«0. Since the gains L can be 
synthesized, the compensator gains Lll,Li2>L2i.I-22 found; however, we must 

assume that |j || is sufficiently sxaall to make both (I^ Dn)“^ & Up+DN^il)”^ 

axis t . 


4.0 aoSED-LOOP STABILITY FOR TEE DPS WITH A FINITE-DIMENSIOm 
COMPENSATOR 

Let e be such that 0 < e < eo» consider the closed-loop system consisting of 
the infinite-dimensional DPS (2.2) and the finite-dimensional compensator (3.1) . 
Let the closed-loop state w in H ■ Hm x R<^ x Hp be given by 


w 


w. 


where wjj 


■['^slis in X R®^ & W|^ ■ v^ in The norm on fi is defined by 

1/2 


! w 


( V. 


I! 2 




N" ■ """ ■ "'R' 

From (2.2) & (3.1), the clos.:d-loop state satisfies: 


3w, 


N 


at 


®11 ”n ^12 ''r 


3t “ ^21 "n ®22 ”r 


where 


®11 


h ®N^iiS 


H. 


12 


^2 lS 


N 11 R 


.^21^ . 


\42 


(4.1) 

(4.2a) 

(4.2b) 


®21 “ ^®rHiS ®R^12^ 

^22 “ ^ * ®R^llS 

But, (4.2) is in the form of the infinite dimensional Kliaushchev-Krasovskii le™<a 
whose proof appears in the revised version of [10]; estimates on the upper bound 
eo are presented there, as well. From the infinite-dimensional K-K lemma, we 
obtain the following closed-loop stability result: 


Theorem 4.1 : If H 22 generates an exponentially stable Cq semigroup and the finite 

rank operator 


s ■ - “12 “22'' S21 


(4.3) 


has all its poles in the open left half -plane, then there exists positive eq (whose 


I 



» 


value can be estimated from the model data) such that, for all 0 < e < cq, (4.2) 
is stable. 


This theorem says that under certain conditions the finite dimensional controller 
based on the ROM (e>0) will continue to stabilize the infinite dimensional DPS 
when e is "small" but non>zero. The "smallness" required of c can be estimated 
from bounds on the control and observation spillover terms, the constants 
in (2.3X and the ROM data. 


The following lemma greatly simplifies Theo. 4.1: 

Lemma 4.2 : H - F « (4.4) 

The proof appears in the Appendix. Consequently, Theo. 4.1 may be restated as 

our main result: j 


Theorem 4.3: If the stable compensator synthesis conditions of Theo. 3.1 are 

satisfied for the ROM (2.4) and if the control and observation spillover coeffi- 
cients £ ■ IjBjJi & r ■ are sufficiently small that 


~ > 6 r 




' 11 ' 


(4.5) 


where Kr & Or appear in (2.3), then there exists positive Cq such that for all 
0 < e tQ the same finite-dimensional compensator (3.1) stabilizes the infinite- 
dimensional DPS (1.1). 


This follows from Lemma 4.2 and use of the semigroup- generator perturbation result 
(1143, Theo. 10.9) applied to 


^22 * ^ ®R^llS 


because is a bounded operator. 


) 5.0 CONCLUSIONS 

Our main result (Theo. 4.3) gives conditions under which a general finite-dimen - 
sional compensator , based on a singular perturbations reduced-order model, will 
stabilize an infinite dimensional (linear) distributed parameter svstem . This is 
an extremely useful result since all practical feedback compensators must be 
I finite dimensional in order to be implementable with on-line digital computers. j 

The result is valid for large-scale lumped-parameter systems, as well. 

The most difficult assumption to satisfy is the choice of reduced-order subspace i 

to achieve a singularly perturbed formulation (2.2) of the distributed para- 
meter system (1.1); the other assumptions are reasonably easy to satisfy. Modal 
I methods have worked for mechanically flexible structures, e.g. [11]. A general 

discussion of this modeling difficulty is given in [15]; it is a fundamental 
problem in large-scale or distributed parameter systems and should not be over- 
looked. j 
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APPENDIX: PROOF OF LEMMA 4.2 


« - ^11 “ ^12 «22 


-1 


‘21 


Ajj + Z L^^Cj, 

Liidp - 

where Z - I^j - 

^ hi’ 


h * s 

S VllS ®N^12 

l4i S hi J 




^22 ■ hi ^R®22 ^R^12 


I 






I 


r 


I 


» 


I 


I 




I 


where ■ 

^hi ~ ’•lA' Hi 

hi " 

A - HA)’" H2 

hi ■ 

HiA ”Ai) 

T m 

i-22 

)"22 * H1V12 

& D • ■ 
N 

■ HH’" H 


Therefore, H ■ F if end only if 
(.) 2 - Oj, - 

'»> V 2 i\ ■ - - hlV 


However, (a) (b) : 

■ =sV'(»22^2'' - ®r'-1iV22'">»S 

" '•R'Sl ^®22 ■ ®R^ll’"il"*22 \ 

■ VV'<V®22‘' 'r ■ ^ 22 "' =R 


Thus it reaains to prove (a) which we recast as: 
id) (I^, + QV)'^ . - Q(WQ 

are both bounded (finite rank) operators. 

Since (WQ + I^) (1^ ^ ^ ^ " ^2'’”^ ^ ^22 ^ assuaed to exist 

for snail enough spillover bounds, therefore 


where 




2^, - QC.Q + Ij)-2 V 


where the required inverse exists. Thus, to prove (a) we need only prove (a'); 
however, this follows directly: 

(I^, ^ QW) - Q(WQ + 

- + Q« - Q(WQ + 1^)‘‘ W 

- Q[WQ(Vp + 

- + QW - Q(WQ + I^)"^ W 

- Q[I^ - (WQ + 

“ V. ^ - Q(WQ + I^)"^ W 

- QW * Q(WQ * 1 )“^ W - 

•• (I^ + QW)‘^ - Q(WQ + I^)‘^ W 


) 


and (a’) Is proved. 
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Abstract 

Largi Spact Structurts art trap tad as a spacial 
clast of highly oscillatory distributed paranatar syt* 
tans. Practical controllars will naad to ba finite di- 
mensional ; however, this calls the closed-loop stabili- 
ty into question due to spillover Interactions. Recent 
results, obtained using regular and singular perturba- 
tions theory, give stability bounds for closed-loop 
control of this type of distributed parameter system; 
this paper surveys these results and their Implications 
for large space structures. 


Ignore these external disturbances to simplify this 
discussion and consider the control forces as generated 
by H Independent actuators: 

M 

Ffx.t) • B-f • I b,{x) f,(t) (1.3) 

0 0 VI ’ ’ 

Ohservatlon Is done with P Independent sensors whose 
output vector Is given by 

y(t) ■ CjU(x,t! * EgU^(x,t) (1.4) 


1.0 Introduction 

Large aerospace structures, such as satellites and 
spacecraft, are necnanically flexible structures which 
win usually recuire active (feedback) control systems 
to suopress' vibrations and satisfy stringent require- 
ments for pointing accuracy, orientation, snape, and 
stat1on-*eeo1r.g. The problems and current directions 
in large aerospace structure control nave been surveyed, 
for exaraole, in II] -[3]. 


1.2 SMte Space Description 


The OPS (1.1), and (1.3)-(1.4) can be rewritten in 
many state space formulations (not all of them are equlv- 
alenth however, the most natural one Is tne state 
v(x,t) * [u(x,t), u^(x,t)] whitn satisfies: 


f v^ ■ Av + Bf 

1 y • Cv 


(1.5) 


Suen structures are hignly oscillatory, distributed 
oaramete*~ s /stems (DPS). The viewpoint we have oeen 
pushing for’ IPS "n general 1$ a practical control the- 
ory, i.e. one t.hat incluOes such practical constraints 
as'a finite number of (localized) control devices and 
finite-dimensional (implementable) controllers; this 
nas been discussed in [4], Other helpful references on 
aspects of DPS t.neory are ^S]-[8]. In this paper, 
large aerospace structure control will be viewed in the 
franeworx estacUsned in [4], and recent stability re- 
sults using -ecu'ar ano singular perturbations tech- 
nicues [9]-]lc] -'ill be examined in the large space 
structure context. 

.aroe 4e*escace Structure Description 


Large aerospace structures may often be modeled by 
the followme DPS description: 






* (x 


( 1 . 1 ) 


wnere F.lx.t) is the displacement of the structure 
< n'’ off Its eouilibrium position in space. The 
damping and stiffness differential operators D© and A^ 
are oefinec on a domain 3(An) which Is a subspace of 
the Hilbert Space Hg • L‘(n] with inner product (•,■)«; 
D(A.) contains all smooth functions In Hg which satisfy 
the structural boundary conditions. Furthermore, the 
ooerators satisfy: 


with the state space being the Hilbert space 

e I W * 

H ■ D(Ag '‘) x Hq witn the energy norm 

llv!!^ • * iiAg’^^ull^ (1.6) 

and the operators (A, B, C) determined from (1.1) and 
(1.3)-(1.4). The input and output operators 3 and C 
have finite rank, and the system operator A has domain 

D(A) • D(A,''^^) X D(A.) and generates tne C. semigroup 

9 3 - 

'J(t). 

The oamoing present in the structure gives tne sys- 
tem a small amount of overall oissioation. It is tne 
job of the finite-dimensional feedback controller to re- 
distribute this dissipation to critical structual mooes 
to improve their performance without creating instabil- 
ities In any residual modes. 

E.D Mocel -Reduction and Controller Svntnesis 

The syntnesis of a finite-dimensional controller is 
based on a reduced -order (finite-dimensional) model 
(ROM) of the DPS. In general this Is obtained by de- 
composing the original state space H into the ROM sub- 
soace Hjj and the residual subsoace Hg: 

H • Hj^ J Hg (2.1) 


^ -'k * 


-0 •< ‘'<-'k'k 


(1.2a) where dim < «. Then projections Vj^ • P,.v and 

Vg > PgV are defined for these suospacis, and the DPS 
(1.5) is decomposed: 


■wnere ^ (x) are tne mode shapes . tne mode *reeuen- 
cies . ano the •ledal tamenq {which i$ generally 
thought tb oe sr^ll for large aerospace structurts). 

Alt.hougn the force distributicn -g(x,t) in (1.1) 
tne sum of :ist.rbance and control forces, we will 



f^it • S ' M - Sr ''r * 

(2.2a) 


'Sai, • ♦ A. Vg * B.‘ 

(2.2b) 

is 

^ ’ Si *N * ^R ''R 

(2.2c) 

* 
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whtrt . P^AP^. A..^ . P^APp, 8„ • P^. ttc. Tht 
ttrms Ajj^ Vj^ tnd Vj^ art called model error , and 
B(jf, C(j Vjj are called control and observation teillover . 

respectively. The parameter c > 0 appears in some 
cases for appropriate choices oT the subspaces and 
Hp; when it appears and is small, we art otaling with a 
singular perturbations fomwlation of the OPS. In many 
eases c ■ i tor the choice o^ subspaces and this is a 
regular perturbations formulation . Both of these form- 
ulaiions wii^ be dealt with Here. A basic, and reason- 
able, assumption made in both formulations is that Af| 

generates a stable C„ semigroup . 

2 . ' Regular Perturbations ROW and Controller Synthesis 

When t • 1, the regular perturbations ROW is ob- 
tained from (2.2) by ■ignor'ing the residua' state v^: 

(v^)t * Afi ''n * 

y • Cj, (2.3b) 

This is a finite-dimensional subsystem (since 

dim Hj, < ■); conseQuently, the following finite-dlraen- 


controller may oe synthesized; 



(2.4a) 

(v^)t * * ®S' ^ 

(2.4b) 

A A 

^ S "n 

(2.4c) 


wi th ■ diag • > * * ^ diag ,. .. 

and bJ, have terms involving the inner products in 

Ho of ectuator or sensor influence functions with mode 

shapes Oi Corresponding parameters appear for 

the residual subsystem (A|,, 6^, Cp) in (2.2). 

The modal controller associated with this medal ROM 
can be obu(ned directly using the parameters (2.7) in 
(2.4). Modal controllability and observability condi- 
tions are available and easily verified, e.g. [3]; con- 
sequently, the controller synthesis may be carried out 
by standard finite-dimensional techniques. Note that 
any subset of critical modes may be used in this ap- 
proach. 

2.2 Singular Perturbations ROM and Controller 
Synthesis 

The singular perturoations formulation of (2.2) 
occurs wnen e > o ano small. )~he singular perturbations 
ROM is obtained from (2.2) by setting c • o and is 
given by 

^''N^t ’ \ * (2.8a) 

y • ^ f (2.8b) 

where - A^p Ap Ap^^, 5'j, • Bp - A^p Ap Bp, 

■ Cp - Cp Ap App , and 7 • - Cp Ap Bp. This is 

clearly di'^ferent ^rwr tne ROM obtaineo in (2.3). How- 
ever, it is a finite-dimensional ROM ano the controller 
associated with it is given by: 


The controller gains 5^ and are designed to 
staoilize and improve tne performance of the ROM (2.3) 
in closed-loop with the controller (2.4). Of course 
the closed loop stability of the actual OPS (1.5) with 
the controller (2.4) will be a principal guestion; reg- 
ular perturoation theory is used to answer this Ques- 
tion in [9], 

For large space structures, tne most popular model 
reduction technioue nas Seen tne modal one, i.t. 


(2.5a) 


(2.5b' 


''N * * 


- ! 


V- • PsV • 


■"r- 


f * (2.9a) 

'^N^t * \ 'n * *N* ' (2.9b) 

y • Tj, Vp - Uf (2.9c) 

This controller differs from that in (2.4); however, 
the gams Ip anc c are sesignec Py standard finite- 
dimensional techniiues to stabilize and improve tne per- 
formance cf the RCM ;2.8l in closed-loop with the con- 
trcller (2.9). Again, tne closed-looo stability of tne 
reouceo-oroer controller (2.9) with tne actual DPS wnen 
c is nonzero is in Question. To attempt to answer this 
funoamental Question a singular perturbations stability 
theory »as developed in [lOj for OPS of the form (2.2) 
with finite-dimensional controllers (2.9); this theory 
gives estimates of tne tolerable size of e. The corre- 
sponding singula' oerturoations theory for finite dimen- 
sional systems ’s surveyed in [Hj-[12]. 


where P? u • : u.Ol, ?S u » : UySL. and u^ 

N k a R K a a 

satisfies (by using (1.2) in (1.1)): 


“k * 


• F^ 

k k c 


'o’*k^ 


. 6 ) 


In this fornjlation. the projections Ps and Pg 
are orthogonal n H and tney ’'reduce" A, i.e. tne 
model error te-'is Ap. • 0 and Aj^ ■ C; of course, this 
depends on kr.,.ing tne true structure modes rather 
than approximations O’* tnem. Sc tne 'egular pertur- 
bations RCM .1.3) O’* the structure nas peraneters: 
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For large space structures, as in most applications, 
there is no preferred singular perturbations formula- 
tion. The general approach nas been applied to space 
structures in [12]. However, as pointed out in [14], 
the nardest task in t.nt application of the singular 
perturbations approacn is casting the problem in an 
appropriate singular formulation (2.2) wnere tne small 
parameter t will nave a connection with pnysical param- 
eters and it will be possible to verify wnetner c is 
sufficiently smell or not; this is true for applica- 
tions of the approacn to DPS, as well as lumped parair.- 
eter systems. 

Here «e will mention two directions • time scaling 
anc 'reouency scai'rg - wmeh lead to singularly oer- 
turoed formulations of tne mocalmodel Qiscussec *or 
structures at the enc of Sec. 2.1. 

In the time-scaling apereecn (see [12] for a general 
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discussion), w« tssuiM that ttw structurt modas of 
(2.6) saparata Into vloritlons on two time scales • a 
slow one and a fast one - whose ratio Is t. I.a. (2.6) 
becomts; 

“k * 2Vk“k * -k -k • ^ 

t^Uj^ * c2C|jU|jU|^ ♦ uj[u^ • k ^ !)♦) (2. 'Ob) 


Some rtlaballng of modes might be rtquirtd to obtain 
(2.10). Our assoclatad singular perturbations formula* 
tion Is obtained by Uking v,. as In (2.6a) and, In 
place of V||, using: 



whart V. Is as In (2.5b). This yields a singular form- 
ulatlon"(2.2) of tna modal modal of the structure with 
c being the ratio of the time-scalts; the mooal param- 
eters are the. same as those given in (2.7) with the ex- 
ception that C^« [C^ cEg] replaces in the residual 

subsystam. The ROK (2.8) is obtained by setting t ■ 0. 
and the corresponding controller (2.9) can be synthe- 
sized; this controller will differ from the original 
modal controller tz.aj. 


In the freouency scaling aboroach (see [15] for a 
general discussion), we assume tna t c ■ for N 

sufficiently large. I.a. the rtsidual structural modes 
nave suffic'lent frequency separation from the ones con- 
tained In the ROm. In this case tna associated singu- 
lar perturbations formulation (2.2) 1$ obtained by 
taking v^ as In (2.5a) eng, in plada of v^, using; 


b 



ol 


l/c 


( 2 . 12 ) 


where V. Is as in (2.5b); the modal parameters arc tna 
same as’'those ln (2.7) with the exceptions; 


A. • ! . * i c£jl replace 

n . TO *0 V I R H j 

/"R ■‘'R- 

Aji anc C-, wnere • dlag •••]< 

% • 2 dlag ...]. and for 

k > N*l. Something interesting happens here when we 
saT c • 0 to obUin the SOM ’2.8) ano corresponding 
controller ;2.S): since ig becomes zero wnen c * 0 
(and A,.. * 0, Agj. • 0 because of tne use of structure 
modes ).^tne singular perturbations ROM (2.8) Is exactly 
tne same as the regular perturoationse one (2.3), i.e. 

• Aj,. • Bj, • C^. and 7 • 0. Consequently, 

for the freouencv scaling accroach the controller (2.9) 
i s exactly tne same as tne orioinal mooal controlTer ] 
when sufficiently many mooes are controlled, kowever, 
the singular perturbations stability results of [1C] 
would seem to yield different conditions on the size of 
c • l/i-^*T than the regular perturbation results of [9]. 


controller synthesis tKhnIques for OPS have appeared 
In [9]-[10] with corresponding stability results ob- 
Ulned from both singular and regular perturbations 
formulations of tha problem. These have been surveyed 
here by considering their Implications for modal con- 
trol of large aerospace structures. Both singular and 
regular perturbations formulations yield useful results, 
but the singular formulation seems extremely versatile 
due to the variety of singular formulations available 
In large space structure applications; two of these 
singular formulations - time and frequency scaling * 
were discussed here. 
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0n«11ne adaptive control 1s essential for Large Space Structures 
(LSS). where the modal parameters are poorly known, due to modeling 
error, or changing, due to variable configurations. It Is especially 
important that such adaptive controllers produce stabilizing controls 
during adaptation due to the small damping present In LSS. However, 
any such controller must be based on a reduced-order model of the LS!i. 
The spillover from the unmodelled residuals, as well as the modeling 
error, can deteriorate the performance of the adaptive controller and, 
if uncomp jnsa wed, this spillover can defsat the whole purpose of the 
adaptive control. 

This paper investigates adaptive control for LSS using direct and 
indirect schemes and points out the mechanisms whereby observation 
spillover can upset the stability of the controller. The framework for 
nonadaptive control of LSS is reviewed and many of the generic problems 
of adaptive LSS control are pointed out within this framework. These 
generic problems must be overcome for successful operation of adaptive 
LSS control . 


1.0 INTRODUCTION 


This paper deals with the basic problems Inherent in adaptive con- 
trol of large space structures (LSS), such as satellites and spacecraft, 
where the structural parameters are poorly known or slowly time-varying. 

With the advent of the Space Shuttle Transportation System, it has 
become possible to conceive very large spacecraft and satellites which 
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could be carried Into space and deployed, assembled, or manufactured 
there. Such LSS would serve a variety of civilian and military needs 
Cl]i [2]. Including electrical energy generation from the solar power 
satellite - a structure nearly the size of Manhattan Island - to be 
constructed In space and operated In earth geosynchronous orbit [3]. The 
control technology needs for such LSS have been discussed In a variety 
of articles, e.g., [4], [5], and the developing LSS control theory and 
technology has been surveyed, for example. In C6]-[10]. 

The size of these structures, their low rigidity, and the small 
damping available In lightweight construction materials combine to make ■ 
LSS extremely mechanically flexible. In theory, LSS are distributed 
parameter systsims whose dimension Is Infinite; however. In practice, 
their dynamics are usually modeled by large scale systems based on 
approximate elastic mode data. Active control schemes for LSS are often 
required to meet stringent requirements for their shape, orientation, 
alignment, and pointing accuracy. Such active control Is limited by the 
capacity of the on-board control computer, the modeling Inaccuracy In cur- 
rent finite element computer codes for analyzing structural dynamics, and 
available control devices (actuators and sensors); therefore, the con- 
troller must be based on some reduced-order model (ROM) of the LSS. 
Fundamental problems of LSS control Include: 

(1) selection of appropriate modes to control for desired system per- 
formance; 

(2) development of ROM for analysis and controller design; 

(3) computation of system model and control parameters; 

(4) controller design with multiple distributed actuators and sensors; 

(5) the number and location of sensors and actuators for efficient 
control ; 

(6) the effect o', and compensation for, residual (unmodelled) modes 
and modeling error on the closed-loop system performance; 

(7) adaptive and self-tuning controllers for LSS with poorly known or 
changing parameters snd configurations. 

Item (7) Is the basic topic of this paper but It must be considered In 
the context of the other Items with which It Is completely Intertwined. 

The need for adaptive control In LSS arises because of Ignorance of 
the system and changing control regimes . The former is due to 


(a) Ignorance of the systMi structure and order, and (b) Ignorance of the 
systttR parameters; the latter occurs because of changing configuration 
of the LSS. Chances In configuration way be due to construction In-space, 
thermal distortion, or reorientation of subsystems, e.g., rotating solar 
panels or sunshlelds; these changes usually produce slowly time-varying 
parameters. Ignorance of the LSS system structure at^ order Is due to the 
fundamental problem of modeling a distributed parameter system, e.g., 
faulty physics, reduced-order models, and Ignored nonl Inearltles; this 
means that the order of the ROM Is lower than that of the actual LSS. 
Ignorance of the system parameters , while directly related to the system ' 
structure, Is due to the Inherent modeling error present In even the best 
structural analysis computer codes and to the limitation of testing huge, 
lightweight LSS on earth; this produces constant but poorly known system 
parameters. There Is a very clear need for an adaptive LSS control metho- 
dology that can begin with the best available computed parameters and 
self-tune Its way toward the correct parameters while stably controlling 
the LSS and, possibly, adapting to variable configurations. 

A great variety of adaptive control schemes exists for lumped param- 
eter. small scale systems [11]; In particular, model reference adaptive 
methods have achieved a great amount of success 1n producing stable, 
convergent adaptive controllers, and adaptive observers for systems whose 
structure Is known and whose parameters are constant but poorly known 
or slowly time-varying, e.g.; [12]-[25j. Adaptive schemes may be direct , 
i.e., the available control parameters are directly adjusted (adapted) to 
improve the overall system performance, e.g., l 25"-'26], or Indirect , i.e., 
the system parameters are identified (based on the assumed system struc- 
ture) and the control commands are generated from these parameter esti- 
jnates as though they were the actual values, e.g., [20], [24], [27]. 

The abundance of adaptive control methods is overwhelming and an 
understanding of the Interrelationships and structural commonality of these 
methods is desperately needed; see, e.g., [28], [29], for some beginnings 
in this direction. Furthermore, the use of such methods on distributed 
parameter or large scale systems, like LSS, is greatly limited by the ROM 
problem - the adaptive scheme must be based on a ROM of the actual system 
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and, hence, the order of the model Is, and must remain, substantially 
lower than the controlled system. In addition. It seems essential that the 
LSS adaptive controller provide a stabllltino control In such highly 
oscillatory systems as LSS . 

Th.1s paper deyelops a framework for LSS adaptlye control problems and 
points out generic problems In the use of the most natural direct and 
Indirect adaptive approaches. In other forms, these problems will haunt 
every use of adaptive control on LSS and must be solved before the valu- 
able benefits of adaptive control can meet the needs of this new applica- 
tion area . A few preliminary attempts at adaptive control for specific 
distributed parameter systens or LSS have been made In [30]*[36]i also, 
for the corresporuling parameter s'* anti fleet Ion problem for distributed 
parameter systmns, see [37]. 

2.0 NONADAPTIVE LSS CONTROL 

Following [6], the LSS may be described by the partial differential 
equation; 

ra(x) u^^(x.t) + Oq u^(x,t) + Aq u(x,t) ■ F(x,t) (2.1) 

where u{x,t) represents (possibly, a vector of generalized) displacements 
of the structure 0 off Its equilibrium position due to transient dis- 
turbances and the applied force distribution F(x,t). The mass distribu- 
tion m(x) Is positive and bounded on :i. The internal restoring forces of 
the structure are represented by u where is an appropriate differen- 
tial operator with domain 0(Ag) defined in a Hilbert space with 
inner product (•••)q« In most cases, A^ has discrete spectrum, 

i.e. , 

*0 *k ■ “k \ 

where are the mode frequencies of vibration and ®|j(x) sre the mode 
shapes . The damping term Dq u^ is generated by an appropriate A^-bounded 
differential operator and may represent gyroscopic damping as well as the 
very small (^%S critical) natural damping expected in the LSS. 

The applied force distribution is given by 

F(x,t) ■ Fj.(x,t) + F^(x,t) 


(2.3) 


s 


I 




. t 


» 


I 


I 


I 


I 


I 


*1 


whtrt Fq rtprtstnts ixttrntl disturbancis and F^ rtprasants tht control 
forcts dut to H actuators: 

• M 

Fj.(x,t) ■ Bq f ■ 5^ b^(x) f^(t) (2.4) 

in 


whtrt b| art tht actuator Influtnct functions (usually point dtvicts) and 
art tht control ccnroands. Observations art produced by P sensors: 

y ■ Cq u ♦ Cg u^ (2.5) 

whtrt yj(t) ■ (Cj,u)g ♦ (Cj,u^)g for 1 < j < P with c^ being tht position 
stnsor Influence functions and cj the velocity sensor ones (usually point 
devices). 

Tht state variable form of (2.1) and (2.3)-(2.5) 1$ obtained by 
taking 

v(x,t) • Cu(x.t), u^(x.t)]' 

In H 2 0(aJ) X Hp with energy norm: 

11 vli^ ■ (mu^, u^) ♦ (aJ u, A^ u) (2.6) 


This produces 


v^ • Av ♦ Bf; v(0) ■ Vg 
y ■ Cv 


(2.7) 


C = :Cg Cg] and A 


where we consider the disturbance- free case (Fg = 0) and define B = [0 Bg] , 

^ 0 I 

This distributed parameter system Is 

very oscillatory In the sense that the semigroup U(t) generated by A has 
very little damping: 

!1 U(t)il < Mg e‘®^ for t 1 0 (2.8) 


where c ^ 0 and small and Mg >, 1. 

The desired performance of the actively controlled LSS greatly effects 
the design of the controller. Many desirable properties of the active 
structure can be obtained with constant feedback gains applied to the sys- 
tem state v(x,t); such solutions arise for regulator problems and sta- 
bilization (pole placement) problems for LSS. However, the full (infinite 


dimensional) state v is never available from a distributed parameter sys- 
tem; only the P sensor outputs y are available. 

Implementable controllers for LSS (and most distributed parameter 
systems) must be based on finite dimensional on-board control computers 
which process the sensor outputs y and produce control commands f; thus, 
a reduced-order model (ROM) of the system (2.7) must be used for the 
controller design. A ROM can be obtained by projecting the system (2.7) 
in H onto an appropriate finite dimensional subspace the projections 
P (onto H, / and Q (onto the residual subspace) are usually, but not 
always, orthogonal. Let v^^ • Pv and Vj^ * Qv and, from (2.7), we obtain: 


h “ ''n ^NR ''r ^ 

^R * ''n ^R ''r ®R ^ 

^ "n ^ ''r 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


where » PAP, Aj^j^ * PAQ, = PB, etc. The terms Bj^ f and Cj^ Vj^ are 
called control and observation spillover ; the terms Aj^j^ Vj^ are called 
model error . The ROM for this system is given by (2.9) and (2,11) with 
* 0 and * 0: 


'N 


Afi v^ + B^ f 


( 2 . 12 ) 


y « C.^ Vj^; v^(0) . Pvg 

The ROM state Vj^ and the residual state Vj^ form the true system state v 

given by: 


2 

with total energy 1| v {j 


I'lWll ^ ^ ^ (if projection is orthogonal) 

All implementable controller designs based on any ROM must be evaluated 
in closed-loop with the actual LSS (2. 7), and it is in this evaluation that 
the effects of model error and spillover due to the residuals become 
apparent. 

If the actual mode shapes are known, the modal ROM is a sensible 
choice: 


<t»i ••• 

and the model error terms Aj^j^ and become zero. Of course, any collec- 
tion of modes could be used; usually, the most easily excited or critical 
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ones win be chosen. However, 1n many cases, the partial differential 
operator A Is too complex to provide closed>form mode shapes. Instead 
finite element approximations of the mode shapes are computed (e.g., 
via NASTRAN),and these approximate mode shapes can be used to form the 
ROM for controller design; note that some model error Is present when 
these approximations are used. Henceforth, we will assume the actual 
mode shapes are available to simplify the discussion but much of our 
analysis remains valid for approximate mode shapes and other types of 
ROM. 


Modern modal control (MMC) for LSS, as developed In [38], uses the 
modal (or approximate modal) ROM and develops a controller consisting 
of a state estimator based on the ROM and a constant gain control law: 


and 


'n ' *N 'n * 

y * V|(^, * 0 


(2.13) 

(2.14) 


This controller design requires the ROM (A|^, Bj^, Cj^j) to be controllable 
and observable for the calculation of control and estimator gains Gj^j, 
Kj^. These conditions, in modal terms, provide insight into the number 
and location of actuators and sensors. From [38], (Aj^, Cj^) Is con- 

trollable and observable, when position sensors are used (Cg *0), If 
and only If 


(1 ) min(P.M) ^ max mode frequency multiplicity in the ROM 

(2) each sub-block of and associated with a mode frequency of 
multiplicity must have rank at least equal to 

Similar results hold for other types of sensors, e.g., velocity, acceler- 
ation, or mixtures of types [38], [391. These results are easy to in- 
terpret in terms of the mode shapes, e.g., if no repeated frequencies 
exist, then the above result says that a single actuator and sensor, 
not necessarily collocated, will do the job as long as neither is located 
at any of the ROM mode shape zeros. Since LSS have many symmetries and 
rigid body modes, It Is not often that a LSS control problem will have a 
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controllable observable ROM with only one pair of devices; this has con- 
sequences for the adaptive control problem to be discussed later. 

Let the estimator error - Vj^ be defined and, from eqs. ( 2 . 9 ) 

-( 2 . 11 ) and ( 2 . 13 )-( 2 . 14 ), obtain 

* ®N ''n ®N S ®N 

■ S *N ^ 'Si ‘'R ''r 

'^R " ®R S ''n ®R S ®N * ^R ''r (2.17) 

This shows the effect of spillover on the closed-loop system: even 
though Aj^ + Gj^, ~ ^ and Aj^ are stable, the closed-loop system 

need not be stable. When either control or observation spillover Is 
absent (Bj^ • 0 or ■ 0), then stability Is assured; otherwise, spill- 
over causes pole-shiftinq and can Induce Instabilities [38], [40]. 

Bounds on the destabilizing effect of observation and control spillover 
were produced In [38] and can be extended to the case where some model 
error and small nonllnearltles are present [41]. Such bounds give an 
Indication of how much spillover the closed-loop system can tolerate. 

A variety of methods have been suggested to reduce spillover [6]. 

One obvious way would be to prefilter the sensor outputs, with a bandpass 
filter, to substantially reduce observation spillover. This alleviates 
the worst pola-sl»1ft1ng problem; bounds cn system performance with con- 
trol spillover alone can be found In [42]. Note that the post-filter 
of the controller outputs could do the same job by reducing control 
spillover; this Interchange of filter and controller is possible due to 
linearity and time-1nvariance. The trade-off with this means of reducing 
spillover is that the prefilter Introduces phase distortion which can 
have a destabilizing effect of Its own. Therefore, a very high order 
filter may be required to keep the phase distortion acceptable; phase- 
locked-loop quadrature filters may be another solution [43]. E ven in 
nonadaptive LSS control, the spillover and model error problem is a funda 
mental one . 

Finally, we should note in this section that digitally implemented 
controls would be based on discrete-time versions of the distributed 
parameter system (2.7). One such version Is obtained by using a uniform 
time step At: 
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V 


v(k + 1) « «v(k) + E f,^ 
y(k) - Cv(k) 


(2.18) 


where ^ = U(At) and E s B 


.At 

/ U(t) dr B and the control command Is 
*' 0 


a constant fj^ over the interval (k 1) At ^ t < k At. Other versions of 
this could be obtained with nonuniform time steps. When the ROM procedure 
of projecting onto the subspace is used, we obtain the discrete-time 
versions of (2.9)-(2.11) ; 


V|,(k 

^ 1) 

■ A 

N 

* *NR 

, Vj^(k) + Ej^ f,^ 

(2.19) 

v„(k 

+ 1) 

■ ♦rn 

1 * «R 

[ ^k 

(2.20) 


y(k) 


v„(k) t C, 


(2.21) 

«N 

> the 

ntodal 

subspace. 

the above become: 


v^(k 

+ 1) 

3 $ 

N 

v^(k) ♦ E|j 

'k 

(2.22) 

v,(k 

+ 1) 

• ‘^R 

* 'o 

8r \ 

(2.23) 


y(k) 


V|ij( k) + C|^ 

v,(k) 

(2.24) 


where (2.22) is the same as that obtained by airectly discretizing the 
ROM in (2.12); if the exact mode shapes are not available, these two 
oiscreti rations may yield different results. In aodition, the sampling 
process can alias residual modes and increase observation spillover and 
the zero-order hold process can spread-out the control command signal 
spectrum and, hence, increase control spillover by increasing the energy 
in the residual mode spectrum; this has been observed and investigated 
in [44]. Therefore, the time discretization is a very important factor 
in the design of implementable LSS controllers. 


3.0 TOWARD ADAPTIVE CONTROL OF LSS 


In order to design MMC, or other controllers, for LSS as proposed in 
the previous section, we must have knowledge of the ROM parameters (A^, 
®N’ parameters are obtained from modal data; they are the 

mode frequencies for and the mode shapes ac actuator and sensor loca- 
tions for and Cj^, respectively. This data is required for three 


(1) to determine controllability and observability of the ROM and, 

hence, to help locate control devices effectively; 

(2) to design control and estimator gains; 

(3} to use in the state estimator's Internal model. 

However, we have noted in Sec. 1.0 the sources of error for this data; 
consequently, a need arises for an adaptive version of the MMC of Sec. 

2 . 0 . 

The most logical and reasonable procedure to obtain adaptive con- 
trollers for a LSS seans to be the following: 

Procedure for Adaptive LSS Control 

(a) choose a “nice" reduced-order model (ROM); a modal ROM 
would be the obvious choice; 

(b) use your "favorite" lumped parameter adaptive control 
scheme; 

(c) design the adaptive controller as though the ROM were 
the actual LSS to be controlled, i.e., ignore the un- 
modeled residual part of the structure; 

(d) use this adaptive controller in closed-loop with the 
actual LSS and hope for the best. 

There is nothing wrong with following this procedure as a best first guess 
- in a way, there is little else that one can do to produce an implement- 
able adaptive LSS controller. 

In some cases, spillover is sufficiently small or enough other mathe- 
matical structure is present in the system, e.g., a high level of damping 
in the distributed parameter system [31] or low level of performance re- 
quired from the controller (increased damping via direct velocity feedback) 
[32], to allow the adaptive controller to operate. However, these situa- 
tions are rare with LSS and one would not like to count on the "generosity 
of nature" or the temporary suspension of Murphy's Law as part of the 
above design orocedure. Consequently, we would add the following items to 
that procedure: 

(e) analyze computer simulations of higher-order models of the 
LSS in closed-loop with the adaptive controller base-' on 
the lower-order ROM (e.g,, simulate more modes than you 
plan to control ) ; 


(f) investigate the specific mathematical mechanisms whereby 
the residual (unmodeled) part of the LSS couples into a 
given adaptive control scheme (e.g.. find out where and how 
spillover affects the adaptive controller); 

(g) obtain mathematical results on the amount of spillover 
and/or model error that can be tolerated in the closed- 

loop system and still provide adequate adaptation and control; 

(h) develop spillover and model error compensation schemes to 
augment the adaptive controller when the residuals cannot 
be tolerated (e.g., when the conditions of (g) are not sat- 
isfi ed) ; 

(i) recheck (g) with this compensation in the closed-loop system. 

We believe that, within the basic framework of LSS control as de- 
scribed in Sec. 2.0, this Augmented Procedure (a)-(i) will go a long way 
toward revealing the problems of adaptive LSS control (and indeed, most 
other adaptive control situations where control must be based on a ROM) 
and will help to focus needed attention on these crucial issues. For 
example, although Step (e) would be done most likely at some point in the 
system development phase of any project as complex as the construction and 
operation of a LSS, often it is done much too late and the design is 
"set in concrete (or in th-'s case, graphite-epoxy}"; the other steps 
(f)-(i) may not be done at all. Yet, ignoring the effects of the resid- 
ual unmodeled LSS can produce some disasterous behavior in the adaptive 
controller; this was pointed out quite clearly in the LSS example in [36]. 

When an adaptive scheme is applied to such a lightly damped, oscil- 
latory system as a LSS, the stability of the closed-loop system during 
adaptation is a necessity ; therefore, we do not view convergence and 
(global) stability results as luxuries for adaptive LSS control and shall 
only consider appropriate those lumped parameter adaptive schemes for 
which such results are available. However, even a globally stable adaptive 
scheme may prove to be unstable when it is used in closed-loop with the 
actual LSS instead of the ROM on which it was based . This is no^ a failure 
of the adaptive scheme; it is a failure to satisfy the mathematical hypoth- 
esis of the stability result associated with the scheme. 
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In the rest of this section, we shall study the use of two very 
general adaptive schemes which seem to illustrate the problems and poten- 
tial of adaptive LSS control : the indirect schemes of C231-[241, which 
use an adaptive observer and operate in continuous time , and the discrete- 
time, direct or indirect, schemes which are based on an Autoregressive 
Moving Average (ARMA) model of the controlled plant, e.g.. C251-C271. 

These approaches represent a good cross-section of available lumped 
parameter adaptive control schemes which have been shown to possess the 
desired stability properties. We emphasize that the point of this sec- 
tion is not to criticize or slander these schemes; rather, we mean to 
point out where the hypotheses of their stability results are violated, 
and must be modified, when we attempt to use them on LSS. We feel that 
consideration of these approaches within the context of the Augmented 
Procedure for Adaptive LSS Control (a)-(i) will illustrate the generic 
difficulties in the application of existing, well-behaved, lumped param- 
eter adaptive control schemes to LSS. 

3.1 Multivariable Systems Converted to Scalar Systems 

The results of many stable adaptive schemes, e.g., [13]-[14], [23]-[24], 
are limited to a single actuator and/or a single sensor; yet, we have seen 
in Sec. 2.0 that most LSS co.'.trcl problems will involve multiple actuators 
and sensors. One way to deal with this (although, admittedly it has its 
drawbacks) is to convert the controllable observable LSS problem via out- 
put feedback into one that is controllable and observable from a single 
actuator and/or single sensor; this can be done with almost any output 
feedback gains [45]-[47]. These gains would have to be based on the best 
available calculated ROM data and the designer must hope that they will 
continue to do their job during adaptation. 

The output feedback modifies the original system (2.9)-(2.11) to be- 
come: 

s' “n * S S' ''r * S ^ 

S ■ (Sn * S S S' " S S S' S * S ' <3-2' 

y ■ 'n ’'n * 'r 'r 

where is the output feedback gain matrix, b^, bj^, Cj^, Cj^ are vectors, 
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and f,y have been renamed. Let be etc. and we have that 

the new ROM (K^, b^, cj) is a controllable, observable single-input, 
single-output system and (3.1) and (3.2) become: 


* Nr N N ^ 
"^R “ Nn ''n N ''r N ^ 


3.2 Indirect Adaptive Controller Design 


We apply the design of the adaptive controller in [24] directly to 
the ROM consisting of (3.4) and (3.3) with the assumption, for now, that 

*NR ■ ° 'r * 

The control law is given by 

f(t) ■ gj ^^{t) + f^(t) (3.6) 

where g|^ is a constant gain vector, f^ is a "sufficiently rich” external 
signal (more about this later), and v^^ is derived from the following 
adaptive observer (or state estimator): 


» F V|,^ + gy + hf; V|^j(0) - Vq 


(3.7) 


where F is in arbitrary, stable matrix and g, h are unknown parameter 
vectors. The appropriate matching conditions are: 


' f + 5*cJ • \ 


(3.8) 


where g*, h* are constant 
-et pj = Cg*"^ h*"'' 

(i* » F p*; p*(0) » p 


Let p^ = Cg*"^ h*"^ Vj^j(O)'*']'^ and note that 


0 


(3.9) 


where F = diag [0 0 F] and, when g » g*, h * h*, we have from (3.8) and 
(3.7): 


v^,^(t) » [M(t) p*(t) (3.10) 

where M(t) 5 e^^^ ' [I|^ y(-r) f(T)] dr and we have used the fact 

that h, f are scalars. 
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Mow (3.7) can be rewritten as 

v^(t) - CM(t) I^] p(t) (3.11) 

Mt) - F M(t) + CI^ y(t) f(t)] (3.12) 

M(0) - 0 (3.13) 

where M(t) Is as defined In (3.10)^and we have yet to produce the adaptive 
law to generate p(t). This adaptive law is given by the following: 

Pit) - F p(t) - a(t) CR(t) p(t) + r(t)] (3.14) 


where p(0) ■ Pq Is arbitrary and 


R(t) 

R(0) 


•q R(t)- F R(t) 


R(t) F + CM(t) I^]‘ [M(t) I^] 


(3.15) 


and 

f(t) » -q r(t) - f'*’ r(t) - [M(t) y(t) 

' (3.16) 

r(0) » 0 

where a(t) 1s the adaptive gain and the constant q is chosen to exceed 
twice the absolute value of the real parts of the eigenvalues of F. The 
adaptive gain Is chosen so that 


a(t) « Y + 1 ju(t)I I 


(3.17) 


where 

u(t) - + (N/2)^ (;y(t)| + if(t)i) (3.18) 

with X positive and F + F^ ^ 2X 


3.3 Convergence Results: What Goes Wrong? 

All of the above is exactly as stated in [24] where it is also shown 
in Appendix I and II that 


y(t) > MM(t)|i (3.15) 

and, with f(t) sufficiently rich in frequencies, there is a t^ such that 
R(t) > 0 > 0 for all t 1 (3.20) 
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In addition, It Is shown in Appendix III that 

R(t) p*(t) + r(t) s 0 (3.21) 

This Is very crucial to the stability results of [24] and It Is here that 
observation spillover (I.e., the fact that we are using a ROM) appears - 
(3.21) Is not valid when C|^ f 0; however, 

R(t) p*(t) + r(t) • 4^(t) (3.22) 

t ""T 

where Aj^(t) s J e’^*^ ^ ^ (t - t) 


Let enj(t) 5 Vj^(t) - Vj^(t) and Ap(t) 2 p(t) - p*(t); then 


e^(t) - [M(t) I^] Ap(t) 

(3.23) 

and 


Ap(t) • [F - a(t) R(t)] Ap(t) - a(t) A^(t) 

(3.24) 

Consider V(t) 2 Ap(t)^ Ap(t) and we obtain: 


V(t) < -Y p V(t) - 2 y Aj^(t) Ap(t) 

(3.25) 


where V(0) - ApJ APq and Ap^ “ Pq • Pq* 

This follows [24] except for the additional term In (3.25); also 
note that 

Let AVfjCt) = v^j(t) - v*(t) and AVj^(t) = Vj^(t) - v*(t) where vj^ and vg 
represent the ideal states of (3.4)-(3.5) when the parameters are 
exactly known: 

"N ■ * »ll sl> ’'N * ’fuR * "n ^ 

'r ■ <’^RN * ‘>N 3n) ’'n \ 'R * *>R % 

with v*(0) ■ v^(0) and v*{0) « ^^^(0). 

When the implementable adaptive control law (3.6) Is used, we obtain: 

'N^ •^''n ^NR ■''r ®N (3.29) 

^''r “ ^^RN ^''r ^R ‘^''r ‘’R ®N 

and we have the following result: 


THEOREM 3.1; Assum* 


(1) Cj^ ■ 0 

(2} X|f + b|^ gj stable 

(3) stable 

( 4 ) ■ 0 

(5) fg(t) sufficiently rich (I.e., It has at least 3N/2 
distinct frequencies) 

2 

Then there is a 5 > 0 such that for all IUpqH <, 6: 

(a) Is bounded and (eventually) vanishes with an 
arbitrary exponential rate 

(b) lira AVyCt) ■ 0 
t— ^ 

(c) 11m AV(,(t) ■ 0. 

t— 

Therefore, even though the closed-loop system with the adaptive 
controller Is highly nonlinear, it Is stable while the adaptation 
Is taking place. In particular, (1), (3), and (4) are satisfied 
If there Is no observation spillover (Cj^ ■ 0) and no irodel error 
■ 0) and some damping in the residuals (i.e., A^^ is stable). 

Proof ; This result follows from the results of [24] because Aj^(t) = 0 In 
(3.22) when c^ ■ 0. The stability of (3.29)-(3.30) is determined by 

S * ^NR " ^R " Nr “ Nr’ N “ N’ 

and cj ■ 0 since it is a row vector of Cj^. # 

Note that the stability of the system with adaptive control is 
determined by that of (3.27)-(3.28) - the ideal case where the parameters 
are known and the external signal f^ is apol led . This is natural, since 
adaptation cannot take place without f^ present; however, after adapta- 
tion, we would most likely want to turn off f^. In addition, we could 
choose the f^ signal so as not to excite the residual frequencies whenever 
sufficient spectral separation is present. Still, we would need to turn 
on f^ now and then, in order to "tune-up” the controller. 

3.4 Spillover Compensation for the Indirect Adaptive Controller 

The above result is merely a slight extension of the results of [24] 
to a special case of the adaptive controller based on an ROM instead of 
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the full system. However, It does suggest that some form of compensation 
should be used to eliminate the observation spillover. Such compensation 
must be essentially Independent of the parameters of the ROM; yet, most 
methods of spillover reduction require knowledge of the R(M (and some 
residual) parameters. 

One approach to spillover compensation already suggested In Sec. 2.0 
Is preflltering the sensor outputs to remove or greatly reduce the term 
C|^ preflltering can be achieved with low*pass or band*pass 

filters when the ROM frequencies are known and separated from the residual 
frequencies. However, the modal frequency data Is part of the poorly 
known parameter data. 

In an attempt to resolve this predicament, we could try using phase- 
locked loop (PLL)-based filters with the center frequency of each loop 
tuned to the best approximation available of the corresponding ROM data. 

The PLL will adapt Itself until It tracks the actual mode frequency, and 
after "lock-on," It will behave as a narrow band-pass, linear filter which 
tunes out the observation spillover from the other frequencies C^3]. Of 
course, sufficient spectral separation must be present, the calculated 
values of the ROM modal frequencies must be sufficiently good, the dis- 
tortion Introduced by the filter must be sufficiently small , and the 
adaptive controller must not shift the poles around too much. Thus, the 
PLL prefllter Is' not a panacea,^ But, It might work to reduce spillover 
and. If It does. It has the added advantage that Its output could also 
reveal better estimates of the modal frequencies; this would take some 
of the load off the adaptive observer. If the adaptation mechanism causes 
too much pole shifting, the ROM frequencies could be excited via f^ and 
Identified In open-loop by the PLL filters before the adaptive controller 
Is turned on. Note that modal frequency data Is usually better known, 
via computer approximation, than modal shape data; hence, this approach 
might not be unreasonable. Another possible, but untried, approach to 
spillover compensation might be an adaptive version of the orthogonal 
filter In [48]. 

Note that some preflltering (and postflltering) always takes place 
due to the bandwidth limitations of the sensors (and actuators). Whether 
this can be used to advantage In adaptive LSS control Is still a matter 
for speculation, e.g., [33]. 

ORIGINAL PAGE Is 

OF POOR QUALHY 


A probitm that ar1s«is with tht ust of tha rtsuUs In [24] and thair 

iMdIfIcatlon to LSS Is that tha constant faadback gains gl must ba ca1> 

^ T ^ 

culatad. In advanca» to stablllza A|^ ♦ b|^ g^^. It would ba battar If 
thasa gains wara adaptad along with tha paramatars In tha obsarvar. Of 
coursa, aftar adaptation has takan placa* thay could ba racalculatad 
from tha ”tunad-up" paramatars, but. In soma casas, tha adaptation phasa 
Is navar ovar, a.g., s1ow1y*vary1ng paramatars. Othar approachas could 
ba usad for adaptiva pola-placamant, a.g., [20], [21], but thasa also 
hava thair limitations. 

3.5 ARMA-6att1n' 

Many discrata-tima adaptiva control schamas dapend on an Auto* 
Ragrassiva Moving Avaraga (ARMA) raprasantatlon of tha plant In dlscrata- 
tlme, a.g. , [25]-[27]: 

N N 

y(k + N) - ^ y(k + r - 1) ♦ ^ f(k ♦ r - 1) (3.31) 

r«l r«1 

for soma H and appropriate matrices o«, S«. What the ARMA says Is that, 

r r 

after N time steps, the present output Is only related to the past N out- 
puts and inputs. Existence of an ARMA Is directly related to the finite 
dimensionality of the plant (N Is usually that dimension) and Is obtained 
using the Cayley-Harallton theorem for matrices. For LSS, only a quasl- 
ARMA can exist; these were considered in detail In [4S]. From tha 
Appendix of [49], we obtain the quasl-ARMA for the LSS (2.18) or (2.19)- 
( 2 . 21 ): 

N N 

y(k + N) • ^ y(k + r - 1) + ^ E^ f(k + r - 1) ♦ R(k) (3.32) 
r-1 r«l 

N 

a(k) • Vj^(k + N) + ^ Vj^(k ♦ r - 1; 

r-1 


where 


and Is dafincd In tht Appendix of [49]. Since R(k) > 0 when ” 0 
and ■ 0» we have the following result: 

THEOREM 3.2: When the observation spillover (Cf^) and the 
model error term are both zero* the quas1*ARMA (3.32) 

Is a true ARMA for the LSS (2. 18), and any stable adaptive 
scheme based on this ARMA will be globally stable when used 
In closed-loop with the actual LSS (2.18). 

When the rather stringent hypothesis of Theo. 3.2 Is not satisfied 
(as It may not be In practice), any adaptive LSS control scheme based on 
the quasl-ARMA (3.31) must Ignore R(k) In order to be Implemcntable. How- 
ever, R(k) Is the term where the residual effects - spillover and model 
error - enter the scheme and can cause Instability. Again, as In Sec. 

3.4, preflltering or other compensation might be tried In an attempt to 
reduce or eliminate this term. 


4.0 CONCLUSIONS 

In an attempt to point out the crucial Issues and generic problems 
associated with adaptive control of large aerospace structures (LSS), we 
have reviewed the framework for nonadaptive LSS control (Sec. 2.0) and, 
withiii this framework, hav* proposed a general procedure, based on 
reduced-order models (ROM) of the LSS, for obtaining and assessing the 
problems of adaptive LSS controllers (3.0). In addition, we have taken 
a closer look at the use of certain well-known, lumped parameter, stable 
adaptive control schemes In chls procedure. Taking these schemes as 
representative of the basic Ideas present In all lumped parameter adap- 
tive control approaches, we have obtained corresponding LSS adaptive 
controllers and found the following generic problems associated with 
adaptive LSS control : 

(1) LSS are distributed parameter or large scale systems; 
therefore, the plant dimension Is always larger than 
the dimension of the adaptive controller, which must bs 
based on a AOM; 
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(2) LSS control must often be done with more than on: actua- 
tor and sensor; conversion of multivariable to scalar 
systems via output feedback introduces problems of its 
own (e.g., stability of the residuals); 

(3) LSS control problems are often non-minimum phase due to 
noncoil ocated actuators and sensors; 

(4) Interaction of the residuals with the adaptive control- 
ler may negate the stabilizing properties of the con- 
troller due to observation spillover; this interaction 
is much worse due to the nonlinear nature of adaptive 
control ; 

(5) Methods of spillover compensation for LSS often require 
knowledge of the ROM (and some residual) parameters - 
the very data that are poorly known; 

(6) The adaptation mechanism may shift the closed-loop fre- 
quencies around; this counteracts the benefits of any 
prefiltering unless sufficient spectral separation is 
maintained; 

(7) Indirect adaptive controllers need sufficient excitation 
from an external signal f^; however, this signal may 
substantially excite the residuals. 

(8) Discrete-time adaptive controllers can only be based 
on quasi-ARMA rather than strict ARMA representations 
of the LSS; this may negate the stability properties of 
such a controller. 

Stable adaptation is essential for such highly oscillatory systems 
as LSS, yet our preliminary stability results, Theos. 3.1 and 3.2, both 
require that observation spillover be somehow completely eliminated 
before it reaches the adaptive control logic; certainly, this is not 
an easy thing to do in general! Perhaps, global stability is too much 
to ask for LSS adaptive control because observation spillover will 
always be present to seme degree in LSS control. However, it seems rea- 
sonable to hope for the development of spillover bounds to give some 
idea of regions of stability for the successful operation of LSS adaptive 
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control . Some comparison should be made between stable adaptive con- 
trollers based on ROM and stable, robust control schemes, e.g., [50], 
[ 51 ]. 

In closing, m would like to say that it is not our intent to present 
a gloomy picture for the application of adaptive control to LSS . In 
fact, the need for adaptive control in LSS is already becoming quite 
clear, and recognition of this need comes, for a change, at an appropri- 
ate time - before any LSS have been built and put into space. However, 
the development of adaptive control for LSS will not take place overnight 
and will not be done by one or two people. Consequently, what we have 
tried to stress here for the interested researcher are some of the funda- 
mental problems that arise and the basic steps which need to be taken 
toward the goal of successful adaptive control of LSS. In the long run, 
we have high hopes for the success of this endeavor and we expect that 
adaptive control theory will profit by its association with large aero- 
space structures, as well. 
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A 3 S-nUCT 

Taa quanticaciva affacca of splllovar on raducad 
crdar adasclva concrol ara txaminad via a :a$c 
axasslt. 

: z.t?odvct;ok 

A aajor aaauspcion in i.ia various ractnc prcofs 
'!i; - ='* luap«^-?ara=atar systan (LPS) adapcive 

::r.:roi algorichos is tha adaquacy of tha ordar a: 
i.'.t iODarallar. Tha niooai of cha plant usac in 
$;acii:.-iag c.ha orcar of the controller is recuirad 
:: rate.-, or axcaad that a: t.ha plant. For axanpla 
, ftr tha aodal rafaranca adaptive control of the 
• »ir. z 

v(i) « :■ a.vC'K-i) - 'a.u(k-j) (1) 

i-1 *’ :•! ^ 

in the absence nf accurate a priori specification a: 
tht platt paraaatars a^ and Si an adjustable concrol 
lav of tha fora 

uU-i; • c,(k-l)r(<-l) * _ c . (k-l)u(k-j) 

- .ii? 1 

n . 

- ' d^ (k-l)y(k-i; C) 

is Tacuired to insure the asyaotocic convergence of y 
ta tha output s of tha aodal 

n 

Sv'k, • gr(k-l) -► j t^. s(k-i). (3) 

1«1 ‘ 

;<ota mat in (2) tha susation lisits a and n must 
correspoed exactly to tha aoving-avaraga and auto- 
regression lengths in (1), raspactivaly . As shown in 
(^) , tha problem of (1) • (3) is solvable via a 
variety of indirect and direct adaptive concrol tech- 
niques chat prove quite siailar in this exact ordar 
case. However, in the raducad order case where the 
ccntrcl input is provided by 

P - 

u(k-l) - t,(k-l)r(k-l) * I c (k-l)u(k-j) 

j-2 J 
c 

- ' a,(k-l)y(k-i) (4) 

i "l 
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with n > p and n > q these slailarlclas disappear 
along with the ability for the output error s - y to 
converge to zero. In this case it is expected that 
the results of indirect anc direct adaptive control 
will be quite eiscinguishaclc. which counters the 
current trend of exciaisiing their equivalence (')-tlO). 

The use of a reduced-order adapcive controller, 
such as (4), can arlsa from a dtsira to rccuca control 
coBplexity with a coosansurate racuctioe in controller 
sensitivity (_^) or from necessity in an attaepc to 
use a finite diaansional ccntrollar on an infinite 
dinensional syscen { 13) . A recent interest (_^) - (15) 
in develoDing adaptive controllers for flexible struc- 
tures by asplication of existing IPS adaptive contrsl- 
le- strategies ro truncated expar.sicn descriptions c: 
the distributed paraaeter systea (S?S) beha\ior of 
flexible structures has led to rwo qualitative descrip- 
tions of the eisbehavior of reduced-order adapcive 
controllers. These interpretations of the additional 
cifficulties facing reduced-orcar adapcive concrollers, 
which are bypassed by axact-ordtr aeaptive controllers, 
are susmarized in the next section. The last section 
formulates a test problem, which initializes attemets 
to quantify these qualitative insign ts. The aost 
current results of this investigation will be compiled 
in the conference presentation. 

II iliDUCc>OrOr?. ZTTiCZS IN .aIAPTI’.I CONTSCl 

Though nuaercus strategies heve been espoused 
( lb ) for extracting a low ordar aodel froa a too com- 
pltx system description, none currently seea fully 
aoclicable to the real-tiae, recursive requirements of 
on-line edaptivc control algorithas. These reduced 
order moccllng techniques seea to fall into rwo broad 
categories: (i) Ixtract tnosc "nodes" (or component- 

subsystems) froa the full systea description that arc 
aost influential in tht perforaance of the aodel in 
its subsequent use. This strategy is followed, for 
example, in (^) and (^) . (li) Parameterize the 
reduced structure to provide the best prediction of the 
desired output. This latter strategy will produce a 
aodal Che "modes" of which need not correspond to any 
of Chose of the full system as ooced in (14) . this 
latter strategy includes the aodel reference approach 
prominent in adaptive systems (19) . 

The presentation in ( 14 ) esphasizes the difficulty 
of reduced-order adaptive control by cootraating these 
two aodel reductiOD strategics. Consider the approxi- 
aacion of 

d(9,k) - I M.(l4)9,(9) , (5) 

j-1 J J 

ORIGINAL PAGE IS 
OP POOR QUALITY 


WP2-0 


9 vh«r« d can tarv« as chc dtflaction of a floxiblt 
•crueeurt, u< cho nodal asplltudos, and dj tha nod* 
shapaa for dI^S or with «< s 1. as cha output of a 
dacouplad (partial fraction axpansion) daacrlpciao of 
a LPS with Uj at tKt individual "nodal” outputs (20) . 
Aasualnt t’bat tha nodal asplltudat obay a diffarenca 
aquation such as 

i 

n 

W,(h) - I [a,,W (k-i) b,.F,(k-i)]. (6) 

J 1,^ ji 1 ji j 

vhara tha Tj ara tha nodal inputs and tha Sj4 and bj^ 
ara tha unknown paranatart of intarast, tha first 
strataty for adaptive Idantification of (6) raquiras 
• axtractioo of tha appropriata Vlj fron (5> foraing 

Mof [l.S] 

d^(e,k) - I Wj(k)Sj(6). (7) 


SHAXF: y(k) • 

a(k)y(k-l) ♦ b(k)u(k-l) 

(18) 

a(kt-l) 

■ a(k) ♦ uv(k)y(k-l) 

(19) 

bOc^i) 

• b(k) pv(k)u(k-l) 

(20) 

v(k) ■ 

y(k) - y(k) ♦ X[y(k-1) - y(k-l)]. 

(21) 


If tha y in (9) ware ganaratad by a first ordar nodal 
as ara (10), (12) and (IS), than aach of tha thraa 
schasMs (if u, s, and 1 wara ehosan to assura sttbi* 
llty) convarga to tha uniqua eorract valuas givan 
tufficiantly rich u. Howavar for tha raducad>ordar 
idantification tasks of (9) - (21) aach schana ylaldt 
dlffaranc naan convarganca points as shown in Pigura 
1. Tha input u in this caaa is whita tad xaro-naan. 
Tha contours in Figura 1 raprasant constant nomslitad 
naan square error 

t{a*(k)} - £{[y(k) - iy(k-l) - bu(k-l)]^) (22) 
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which as noted in (^) , (^) and (_2£^ tanains an open 
problem without knowledge o; N and the ty tha 
second stratagy would astisata the Wj such that 

Mof[l,.S-]. 

d(e,k) • I w_(k)»j5) (8) 

j 

closely approMaated d in (5) in, c.g., a least* 
squares sense. The astinatad nodal amplitudes 
need not provide useful values for the cstiucion o: 
the a<j and b-i in (6) with W, replaced by W, as noted 
in . In IS) tna chosen at each time instant 
to closely match c to c may not even obey an under- 
lying difference eeuation such as (6). The problem 
with closing tha aoaptive control loop is tnen inten- 
sified by tha difficulty in idancifying (6) and then 
in feeding back the w< instead of the U, for modal 
control. 

Tha difference in various identification tech- 
niques in the reduced-order L?S case can be amphasizad 
by a simple example froc (Zl) . Consider identifying 
tha stable, second-order plant 

y(k) • 1.131-y(k-i; - 0. If y ; -.-I : * 0.05u(k-l) - 0.-u(,k-l) 

(9) 

uaing a first ordar nodal and three cooson adaptive 
identification schemes, tna equation error fonsulatcq 
UiS algorit.hm (22), the gradient based output error 
formulated Steam's algorithm (SA) (21) > and the 
stability based output error formulated SHARF (2c) . 

The algorit.nms are 


UiS: 

y(k) • 

a(k)yCx-l) b(k)u(k-l) 

(10) 


i(k+l) 

■ a(k) + u(y(k) - y(k) ]y(k-l) 

(11) 


b(k+l) 

“ b(k) * cly(k) - y(k)]u(k-l) 

(12) 

Sa: 

y(k) - 

i(k)y(k-l) t. o(k)u(k-l) 

(13) 


i(k*l) 

- a(k) + u(y('x) - y(k)]X(k) 

(14) 


Kk) • 

y(k-l) + a(k)X(k-l) 

(15) 


bCcel) 

• b(k) 0[y(k) - y(k) ]y(k) 

(16) 


TOO - 

u(k-l) T- a(k)T(k-l) 

(17) 


for various choicts of a and b. Note that tha 
gradient-based SA is susceptible to local convarganca 
and for this problem tha SA mean convergence point is 
dapandanc on the initial guesses for s and b. Though 
if any of these schemas is used for simultaneous 
idtntif icsticn and control the input would no longer 
be white, clearly aach would have a quits eiffartne 
control barfomance. I: becomes unclasr as to which 
wculd adabtivtly provide iha bast rcauced-ordar modal. 

An alternate interpretation of the "dual" ill- 
effects of reduced-order adaptive control has aaargad, 
also iron the fla.xibla structure control preblaa ( 12) 

' If) . "or this viewpoint consider saparatisn o: a 
single- input . single-output system state-space de- 


scripcion of 

:h« 

raducad-oroar 

modal and tha ramaindcr 

(or risidual) as 






+ i • lu(k) (23) 

■ 


1 x^(k+l) ! 

i^N 


ibRi 

y(k) - c,j 


r^^.(x)“ 

{24) 


^(k)j. 

From (23) - (24) the "spillover" of the residual modes 
into the observation of y(k) via c^xsik) and the spill- 
over of the control designed for the raducec modal 
into rhe residual modes x^(k-*l) via bqu(k) ara clearly 
displayed. Also the possible coupling of the reduced 
and model modes via A^p and Aq^' is lamcdistely 
separent. Shown in (12^ ( If ) cs tne predictable 
fact that if Asu - 2 ‘R “ £ ihen, atsuming the 
residual modes art not unstable, any full-ocdar-stable 
adaptive controller Identifying A^, bs, and c^ 
axplicltly or implicitly from y and u alone would 
retain its stability. If A«q ■ 2 cr and bj ara 
nonzero then tha degradation of reduced-order adaptive 
control has two sources. Unmodalad componancs in / 
via cqxq will ganaraca an error that ia indistiaguish- 
abla from paraaetcr error cauaing adaptation. The 
application of control u to xr via bqu will causa 
further axeitaclon of tha raaidual modaa further per- 
turbing tha reduced modeling of y. .Sot only will tha 
paraaattr estimates become incorrect but also the 
state axcimatas from an adaptive observer will be 
incorrect leading to an uapradictablt "controlled" 
casponsa. 
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(30) 


# Thu« two problaw of Inopproprlae* p«ru«c*t 
«nd itttt ••ciaatlon art eht^taea ones netad in cha 
aarXiar intarpracacion with Vj rathar chan Wj 
axtraction. Tha ability co obtain dp corraaponda to 
an affactiva caroing of cp. Sinca auch ability doaa 
not appaar forthcoeing in ganaral, tha affacta of 
apillovar on raducad-ordar adaptiva concTol via tha 

^ taesne ttratagy notad at tha baginning of thia 

* aaccion muat ba avaluacad. a taac axaapla follo%ra 
vnich if daaignad to provida inaight into thia 
pcoblaa. 

Ill TZST ZXAHnZ 


Conaidar controlling tha aacond-erdar plant with 
f tranafar function 



(25) 


and tharafora diffaranca aquation daacclpcion 



(a^*a,)y('n-i) - (aj^*,)y(k-2) 

- “ (bj^a,-cb,a, )u(k-2) 


(26) 


2'“' b(kfl) 




and 


d(k) - f^(k)r(k) + fj(k)y(k). (31) 

Tha coBBon adaptiva conctollac froa 

f (k) - f,(k-l) I- or(k-l)[dr(k-l) cy(k-l) - y(k)] 

/ , (32) 

f,(k) • f,(k-l) + uy(k-l) [dr(k-l) ♦ ey(k-l) - y(k) ) 

^ (33) 

with 'j and p approprlacaly conacrainad ui2J. para- 
aatarixa (31) for tha diract caaa. 

Pralialnary raaulca arc aabigueua. Mona of tha 
aiaulacad caaca claarly aupport tha unlvaraal uaaful- 
naaa of naglacting high fraquancy aodaa, tha aupari- 
orlcy of aodal rajaccion or actaaptad "full" adaptiva 
cracking, or tha claar banafit of indiracc or diract 
adaptiva tachniquaa. In othar worda, ona approach ia 
battar in ona axaapla than anothar, which raquiraa 
furthat invaatigaticn of raducad-ordar adaptiva 
control. Tha aixad avicanca will ba intarpracad ia 
:.na confaranca praaantation. 


jndar tha aaauapcion of a firat order awcal. I'aing 
tha aaoaration taehniqua of (23) - (2i), (26) can oa 
rawrittan at 


9 


• X. ( ii*l) 

* 

x.(k-i) ^ 


;‘l 2 : |X.(k1 , 

I , i * i 

3 a,':x, (k'' 



(2T) 


v(k) • J. r:x. (k) 


(28) 
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